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I.  RESEARCH  HIGHLIGHTS 


In  August  of  1981,  we  recieved  approval  from  the  Army  Research  Office  (ARO)  J 

i 

to  conduct  research  on  the  propagation  of  millimeter  waves  through  the  | 

atmospheric  boundary  layer.  Georgia  Institute  of  Technology  (GIT),  whose  exper-  | 

tise  in  millimeter  wave  sources  and  receivers  is  reknown,  were  co-investigators.  i 

The  GIT  responsibility  was  to  develop  transmitters  and  receivers  in  accordance 

with  specifica  tions  in  our  jointly  derived  experiment  design.  Investigators 

the  National  Oceanic  and  Atmospheric  Administration  designed  and  implemented  the 

meteorological  support  stations  and  would  perform  the  bulk  of  the  data  analyses 

and  the  theoretical  support.  The  overall  mission  of  the  experiment  was  to 

establish  a  data  base  for  propa  gation  of  millimeter  waves  at  a  variety  of  fre-  | 

quencles  in  a  well-documented,  wide-variety  of  atmospheric  conditions. 

Before  embarking  on  such  an  ambitious  research  program  as  the  NOAA/GIT 
experiment,  we  performed  more  than  one  year  of  preparatory  research  and  site  | 

selection.  This  latter  task  was  accomplished  after  an  examination  of  more  than 
twenty  sites  (see  Table  I)  over  a  period  of  eighteen  months.  This  was  necessi¬ 
tated  because  of  the  stringent  requirements  necessary  for  efficient,  economic 
operation  and  to  ensure  the  correct  conditions  for  the  applicability  of  the 
Honln-Obukov  similarity  theory,  that  allows  us  to  accurately  characterize  the 
boundary  layer  processes  attendent  at  such  a  site  (Table  II).  We  settled  on  an 
exceptional  site  at  Flatvllle,  Illinois  (see  Fig.  1)  and  signed  contracts  with 
the  affected  farmers  in  early  1983.  At  that  time,  the  equipment  was  deployed  on 
the  1.4  km  Flatvllle  site  as  depicted  in  Fig.  2. 

Before  selecting  instruments  for  final  deployment,  we  performed  preliminary 
tests  on  both  the  millimeter-wave ^  ^  and  atmospheric  instuments.^’^'^  A  major 
result  of  the  NOAA  effort  was  to  redesign  the  Lyman-alpha  humidiometer  to 
Improve  its  frequency  response.  Our  design  changes  were  adopted  by  the  manufac¬ 
turers.  During  this  period  GIT  successfully  designed  and  constructed  the 
transmitter  and  receiver  vans  and  antennas.  This  was  no  simple  task  because  of 
the  stringent  receiver  stability  requirements  of  maintaining  less  than  one  wave¬ 
length  of  random  oration  of  the  receiver  array  over  the  approximately  twenty-six 
thousand  wavelength  long  array.  The  final  design  incorporated  a  13-m  I-beam, 
the  supports  of  which  were  bolted  to  a  concrete  pad  beneath  the  receiver  van 
(see  Appendix  I). 
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TABLE  I:  SITE  SURVEY 


FIELD  SITE  CANDIDATES  VISITED  - 

Illinois  River  Valley  near  Havana,  Illinois 

Kilbourne,  Illinois 

Chanute  AFB,  Rantoul,  Illinois 

Flatvllle,  Illinois 

Champaign-Urbana ,  Illinois 

-  Sandusky ,  Ohio 
Wallops  Island,  Virginia 

RedsCone  Arsenal,  Huntsville,  Alabama 
Great  Bend,  Kansas 

-  Sikeston,  Missouri 

-  Otis  Air  Force  Base,  Massachusetts 

Kennedy  Space  Flight  Center,  Cape  Canaveral,  Florida 

FIELD  SITE  CANDIDATES  RULED  OUT  BY  TELEPHONE  - 

-  O'Neill,  Nebraska 

North  Air  Field  (USAF),  North,  South  Carolina 

Aberdeen  Proving  Ground,  Maryland 

Moody  AFB,  Valdosta,  Georgia 

Robins  AFB,  Warnei -Robins ,  Georgia 

Fort  Rucker,  Alabama 

Fort  Stewart,  Georgia 

-  Fort  Benning,  Georgia 

-  Morris  AAF,  Fort  Gillen,  Atlanta,  Georgia 

-  Whiteman  AFB,  Missouri 
Eglin  AFB,  Florida 


TABLE  II:  SITE  EVALUATION  FACTORS 


Topography  -  flatness 

Fetch  -  lack  of  trees,  direction  of  prevailing  wind 
Annual  precipitation  and  its  predictability 
Occurrence  of  snow  and  fog 
Availability  of  electrical  power 
Access  roads 

Land  use:  crops  and  their  uniformity  over  section 
Number  of  owners  involved 
Cost  to  use 
Security 

Fences  and  other  interference  objects 
Telephone  lines 

Accessibility  via  common  carrier 
Temporary  housing,  e.g.,  distance  to  motels 
Storage  facilities 

Proximity  to  mm-tech  facilities,  NWS  observers 
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Optical  Cj!  and  gf  Optical  Roin 

Crosswind  8  Gauge 

(Paths  1  Km  each)  (Paths  50  m  each) 


After  the  redesign  of  the  Lyman-alpha  humid iometer,  a  preliminary  test  was 
run  to  choose  the  optimum  available  humidity  sensor  from  a  group  of  candidates.7 
Another  equally  important  aspect  of  this  work  was  to  gain  familiarity  with  the 
instruments  in  a  boundary  layer  experiment  and  to  see  if  we  could  eliminate  the 
need  for  the  labor-intensive  radio  ref ractometer  needed  to  measure  the  refrac¬ 
tive-index  fluctuations  in  the  near  millimeter  wave  range.  It  was  concluded 

that  it  was  possible  to  derive  the  refractive-index-fluctuation  statistics  in 
2  2 

terms  of  and  from  our  theoretical  calculations  (Appendix  II  and 

ref.  8),  using  a  combination  of  Lyman-alpha  humidity  and  fine-wire  temperature 
measurements,  with  sufficient  accuracy  for  our  purposes.  A  noteworthy  result  of 
the  preliminary  experiment7  was  £he  strong  degree  of  absolute  correlation  bet¬ 
ween  the  humidity  and  temperature  fluctuations.  The  correlation  nearly  +1 
during  the  daytime  when  both  temperature  and  humidity  gradients  were  negative 
and  nearly  -1  when  the  temperature  gradient  reversed  a  night  due  to  radlational 
cooling  of  the  ground.  The  transition  periods  were  quite  brief  on  the  order  of 
15  minutes. 


Also  during  the  preparation  phase,  we  made  progress  on  theoretical  studies 

9 

of  propagation  of  millimeter  beam  waves  in  complex  refractive  media  and  the  use 
of  optical  and  millimeter  wave  devices  to  measure  surface  layer  fluxes.^  This 
latter  work  is  significant  because  we  believe  that  a  combination  visible,  infra¬ 
red  and  millimeter  wave  system  can  give  path-averaged  values  of  heat,  momentum 
and  moisture  fluxes  that  may  be  more  representative  and  certainly  more  quickly 
measured  than  results  using  conventional  in  situ  devices.  This  technique  would 
be  extremely  advantageous  for  use  at  inaccessible  sites  such  as  irregular 
terrain  or  over-water  paths.  To  test  the  ideas  proposed  in  Ref.  10,  was  one  of 
the  NOAA  groups'  major  goals. 

At  the  conclusion  of  the  above  research,  we  decided  on  a  set  of  meteoro¬ 
logical  instruments  to  be  deployed  on  site  at  Flatvllle  at  the  two  locations 
along  our  1.4  km  propagation  path  shown  in  Fig.  2.  Table  III  lists  the 
available  instruments  and  their  associated  data  recording  channel  numbers.  The 
data  flowed  to  recorders  in  the  receiver  trailer  via  a  buried  fiber-optic  cable 
from  both  Met.  Station  A  and  Met.  Station  B.  Each  station  consisted  of  the 
standard  in  situ  instrumentation  listed  as  channels  0-12  and  deployed  on  the 
meteorological  tower  as  shown  in  Fig.  3.  In  addition,  three  optical  remote  sen¬ 
sors,  channels  13-15,  gave  path-averaged  values  of  path-crosswind,  optical 
refractive-turbulence  strength  and  rain  rate. 


Propeller-vane  anemometer 

Three-axis  sonic  anemometer  with  wire  temperature  sensor  in  center 

Lyman-a  hygrometer  with  wire  temperature  sensor  in  center 

Rotatable'boom 

Vaisala  humidity  probe 

Aspirated  psychrometer 

Electronics  cabinet 


TABLE  HI.  METEOROLOGICAL  INSTRUMENTS 


Channel  // 

Station  A 

Station  B 

0 

Lyman-a  humidiometer 

Lyman-a  humidiometer 

l 

T'  (temperature  sensor)  for  Ly-a 

T*  (temperature  sensor)  for  Ly-a 

2 

T'  for  sonic  anemometer 

T'  for  sonic  anemometer 

3 

psychrometer 

psychrometer 

4 

disdrometer 

Vaisala  temperature  probe 

5 

Vaisala  relative  humidity  probe 

Vaisala  relative  humidity  probe 

6 

sonic  anemometer  U 

sonic  anemometer  U 

7 

v 

v 

8 

W 

W 

9 

angle  encoder 

angle  encoder 

10 

barometer 

pyranometer 

11 

propellor 

propellor 

12 

vane 

vane 

13 

optical  cross-wind 

optical  cross-wind 

14 

rvnb  i  o  1  /"*  2 

v,n 

optical  Cn2 

15 

optical  rain  gauge 

optical  rain  gauge 

One  of  che  optical  rain  gauges  gave  both  rain  rate  and  drop-size  distribu¬ 
tion  (channel  4,  station  A).  The  rain  sensor  data  were  compared  with  a  set  of 
six  weighing-bucket  rain  gauges  deployed  as  shown  in  Fig.  2.  The  optical  rain 
sensor  is  discussed  in  detail  in  Appendix  III.  On  site  we  used  a  50-m  laser 
path.  By  the  spring  of  1983  we  were  recording  both  millimeter  wave  and  meteoro¬ 
logical  data.^  Once  on  site  we  carried  out  five  experimental  phases.  These 
data-taking  sessions  were  in  June-July  1983,  Nov. -Dec.  1983,  Feb. -March  1984, 
June-July  1984,  and  Feb. -March  1985.  These  varied  sessions  allowed  us  the  maxi¬ 
mum  opportunity  to  take  multi-wavelength  data  during  clear-air,  rain,  fog,  sleet 
1 2”  1 6 

and  snow.  More  accurate  dates  and  precise  millimeter  frequencies  are  shown 

in  Table  IV. 

Fundamental  theoretical  work  on  propagation  theory^  was  undertaken  to 
verify  that  the  mutual  coherence  function  (MCF),  derived  for  optical  waves  pro¬ 
pagating  through  the  atmosphere,  could  be  applied  to  the  millimeter-wave  case. 
Optical  waves  readily  satisfy  the  condition  X  «  £q,  where  X  is  the  wavelength 
and  lQ  the  smallest  size  refractive-index  irregularity  in  the  medium,  necessary 
for  the  validity  of  the  parabolic  (narrow-angle  scatter)  approximation,  milli¬ 
meter  waves  do  not.  Reference  (17)  proves  that  there  is  a  broader  range  of 
applicability  to  the  optical  result  and  the  optically  derived  MCF  applies  in  the 


TABLE  IV.  SCHEDULE  OF  FLATVILLE  OBSERVATION  PHASES 


Session  I:  4  June  -  21  July  1983 

Tapes  3-15:  f  -  116.300  GHz  11  -  16  July 

Tapes  16-29:  f  -  172.910  GHz  17  -  20  July 

Session  II:  31  Oct.  -  10  Dec.  1983  (Break:  19-29  Nov.) 

Tapes  2-32:  f  -  116.300  GHz  11  Nov  -  2  Dec 

Tapes  33-45:  f  *  173.001  GHz  3-6  Dec 

Session  III:  5  Feb.  -  9  March  1984 

Tapes  1-33:  f  =■  142.000  GHz  11  -  20  Feb 

Tapes  34-75:  f  -  116.300  GHz  22  Feb  -  2  Mar 

Tapes  76-93:  f  -  142.000  GHz  2-6  Mar 

Session  IV:  31  May  -  4  July  1984 

Tapes  1  -  20:  f  =  142.004  GHz  8-18  June 

Tapes  21  -  22:  f  -  118.754  GHz  19  June 

Tapes  23-24:  f  -  120.750  GHz  20  June 

Tapes  25-26:  f  -  118.754  GHz  20  June 

Tapes  27-65:  f  -  142.004  GHz  21  -  30  June 

Tapes  66  -  70:  f  -  173.005  GHz  1-2  July 

Tapes  71-76:  f  -  172.620  GHz  3  July 

Session  V:  31  Jan.  -  24  Mar  1985  (Break:  22  Feb.  -  4  Mar.) 

Tapes  1-13:  f  -  229.250  GHz  8-11  Feb 

Tapes  14  -  78:  f  *  230.441  GHz  11  Feb  -  18  Mar 

millimeter  wave  region.  This  was  confirmed  in  our  subsequent  data.  Further 

theoretical  work  on  absorption  sharpened  our  understanding  and  the  accuracy  of 

lg_20 

absorption  spectrum  in  general,  and  the  refractive  and  absorptive  behavior 

of  the  atmosphere  and  its  constitutent  gases  in  the  millimeter  wave  region  in 

,  18,19 

particular. 

After  collecting  291 , ..thirty-nine-minute  magnetic  tapes,  each  containing  60 
channels  of  either  meteorological  or  millimeter  wave  data  at  more  than  10  dif¬ 
ferent  frequencies,  a  massive  data  collection,  handling  and  processing  problem 

21 

presented  itself.  The  processing  architecture  is  documented  in  detail  in 

Appendix  IV  and  in  Reference  22.  Procedures  used  to  pre-process  data  to  reduce 

23  24  23 

noise  problems  ’  and  algorithms  to  process  rain  data  were  developed.  The 

26 

final  result  is  a  database  of  answer  files  where  time  series,  spectra,  varlan- 


m 
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ces  and  higher  order  moments,  probability  density  functions  etc.  of  both 
millimeter-wave  and  meteorological  data  are  stored  for  ready  access  by  the 
NOAA-GIT  team  or  other  interested  users.  This  so-called  users  guide  is  the  sub¬ 
ject  of  Section  11  of  this  report. 

Although  the  database  has  been  tapped  only  modestly,  we  have  performed 
several  theoretical  studies  and  tested  them  using  our  Flatville  data. 

References  27  and  28  compare  the  observed  clear-air,  millimeter-wave,  amplitude 
and  phase  difference  spectra  with  a  theoretical  model  derived  by  the  authors 
from  work  based  on  single-scatter  theory.  The  results  show  a  very  close 
agreement  between  theory  and  experiment  and  also  indicate  that  the  millimeter 
wave  link  can  act  as  an  accurate  monitor  of  the  path  crosswind  in  a  fashion  ana¬ 
logous  to  our  on-site  optical  remote  sensors  shown  in  Fig.  1.  Other  work  ana¬ 
lyzed  phase  difference  power  spectra  of  infrared  waves,  to  determine  their 

29 

sensitivity  to  rain  rate  and  eventual  use  as  a  long-path  rain-rate  monitor. 

30 

Another  study  reported  a  detailed  analysis  of  an  oscillatory  rain  rate  fluc¬ 
tuation  and  showed  it  correlated  nearly  perfectly  with  intensity  fluctuations 
observed  simultaneously  on  the  millimeter  wave  link.  Other  work,  conducted  pri¬ 
marily  by  the  GIT  staff,  analyzed  the  Flatville  observations  for  effects  of 

31  32  33 

severe  weather,  fog  and  snow.  *  A  recent  paper  compares  the  behavior  of 
variances,  structure  functions,  and  probability  distribution  functions  of  ampli¬ 
tude  and  phase  differences,  and  mutual  coherence  functions,  with  the  flrst- 
order-slngle-scatter  theory  prediction  using  the  observed  meteorological  data, 
and  finds  excellent  agreement  between  theoretical  prediction  and  experiment  in 
nearly  all  respects.  Finally,  Ref.  34  describes  the  first  demonstration  of  a 
successful  derivation  of  the  fluxes  of  heat,  momentum  and  moisture  from  multi- 
wavelength  scintillation  data.  These  results  were  compared  with  directly 
measured  fluxes  using  the  meteorological-station  data,  and  good  agreement  was 
found.  This  was  one  of  the  key  results  of  the  experiment  and  further  work  in 
verifying  the  results  under  a  wide  variety  of  weather  conditions  is  under  way. 

Although  we  have  just  begun  to  tap  the  data,  many  theoretical  predictions 
have  been  confirmed  and  many  tantalizing  new  results  observed.  The  ultimate 
value  of  the  database  will  be  determined  by  Interested  researchers  both  inside 
NOAA-GIT  and  outside  agencies.  We  have  taken  great  pains  to  make  the  database 
accurate,  accessible  and  available.  If  Army  and  other  workers  will  use  the 
results  of  our  efforts  and  expand  the  research  activity,  we  are  convinced  the 
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Flatville  data  set  will  be  more  than  worth  the  extraordinary  effort  it  took  to 
produce  it. 


II.  DATA  BASE  USER  GUIDE 


II. 1  General  Descriptions  of  the  Data  Base 

The  Data  Base  described  herein  is  a  collection  of  statistical  quantities, 
probability  distributions,  and  spectra  of  data  taken  by  the  NOAA-GIT  mm  wave 
experiment  held  in  Flatville,  Illinois  between  October  1983,  and  March  1985. 

The  purpose  of  the  study  is  to  provide  a  general  purpose  scientific  data  base 
which  contains  quantities  essential  to  the  investigation  of  millimeter  wave 
propagation  in  the  turbulent  and  turbid  atmosphere.  Thus,  the  data  base 
necessarily  contains  meteorological  data  designed  to  complement  the  millimeter 
wave  data.  A  detailed  description  of  the  experiment  and  data  analysis  is 
described  in  Ref.  21  (Appendix  IV).  The  meteorological  data  contained  in  the 
database  is  weather  dependent  because  of  the  extended  use  of  state-of-the  art 
Instrumentation.  For  example;  temperature  fluctuations  were  measured  using  2 
micron  diameter  platinum  wires,  and  humidity  fluctuations  were  measured  using 
Lyman-alpha  hygrometers.  During  precipitation  events,  these  instruments  were 
covered  because  the  rain  drops  would  break  the  fine  wires,  and  the  excess 
moisture  would  saturate  the  Lyman-alpha  gap  rendering  these  data  meaningless. 

In  such  circumstances,  structure  parameters  are  not  calculated  and  only 
macroscopic  meteorological  data  is  entered  in  the  database.  These  data  are 
derived  from  instruments  which  respond  too  slowly  to  detect  turbulence  scales. 

The  database  resides  on  the  CYBER  180/855  at  the  National  Bureau  of 
Standards  in  Gaithersburg,  Maryland.  It  is  not  portable  because  we  have  taken 
advantage  of  vendor  software.  We  decided  to  do  so  after  learning  of  NOAA's  com¬ 
mitment  to  the  purchase  of  CDC  hardware  and  believe  that  this  data  will  be 
available  for  a  long  time  to  come. 

The  system  design  of  the  data  base  is  "user  friendly";  which  in  this  case 
has  literal,  meaning. . .all  items  are  accessed  by  their  name.  Furthermore,  one 
may  select  data  which  falls,  into  a  predefined  category.  For  example,  one  may 
wish  to  collect  all  the  mutual  coherence  functions  (a  complex  quantity)  for  all 
blowing  snow  data  at  a  specific  mm  wave  frequency.  This  is  done  by  setting  keys 
which  selectively  gather  records  from  which  data  is  extracted  by  name.  This 
report  explains  how  to  do  this. 
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The  data  base  is  a  specially  structured  random  access  file  in  the  sense  of 
a  FORTRAN  '77  definition.  The  special  parts  of  the  file  are  twofold.  On  the 
first  part,  software  provided  by  CDC  Corp.  was  used  to  initialize  and  make 
entries  to  the  file  making  it  a  true  database.  The  second  part  of  it's  special 
nature  is  a  concept  of  an  index  and  dictionary  we  developed  and  built  into  the 
file.  The  index  is  akin  to  the  table  of  contents  one  finds  in  a  standard 
reference  book  in  that  it  contains  lists  of  categories  that  can  be  used  to 
selectively  accumulate  data  from  records.  The  dictionary  contains  all  the  data 
names,  their  location  and  length.  The  index  and  dictionary  items  are  prede¬ 
fined.  If  one  requests  a  data  category,  and/or  a  data  item  which  is  not 
defined,  the  software  will  stop  softly  describing  the  erroneous  entries. 

An  optional  product  of  our  data  base  interrogation  software  is  a  smaller 
file  containing  the  categorical  items  requested  by  the  user.  This  file  can  be 
sequential  access  with  records  of  varying  lengths,  or  it  can  be  random  access. 
The  random  access  option  produces  an  ANSI  standard  FORTRAN  '77  DIRECT  ACCESS* 
file;  it  is  portable  only  to  the  extent  that  the  CYBER  60-blt  words  can  be  pro¬ 
perly  converted  for  the  use  by  the  recipient  system.  The  random  access  file 
feature  is  Intended  for,  but  not  restricted  to,  use  on  our,  or  other,  CYBER 
installations. 

The  sequential  access  file  option  can  be  formatted,  or  unformatted.  An 
unformatted  sequential  file  contains  positionally  placed  binary  data  and  has  the 
advantage  that  precision  is  maximized,  but  the  disadvantage  that  conversion 
software  must  be  used  to  move  the  data  to  another  system.  Formatted  data 
however,  transfers  most  easily  using  a  wide  variety  of  commercially  available 
text  file  communications  software.  We  recommend  that  investigators  transferring 
the  data  to  a  smaller  personal  computer  use  the  formatted  sequential  access 
option. 

The  database  can  be  queried  for  named  data  belonging  to  categories  of  items 
that  are  defined  in  the  index.’  The  user  is  not  restricted  to  one  category,  or 

*  DIRECT  ACCESS  IN  FORTRANese,  but  we  avoid  this  terminology  because  of  the 
special  definition  of  DIRECT  ACCESS  files  when  one  refers  to  a  CYBER. 
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one  index.  There  are  seven  indices  each  with  several  categories.  There  is 
complete  flexibility  for  the  user  to  retrieve  data  from  a  single,  or  from  a 
range  of  categories  in  any  combination  of  the  seven  indices. 

Arrays  are  stored  in  contiguous  regions  of  the  record.  The  FORTRAN  default 
memory  organization,  columnwise  storage,  is  used.  Since  the  data  are  known  to 
the  system  by  name,  the  user  only  has  to  allocate  space  to  the  array  items  and 
let  the  interrogation  software  fill  the  arrays.  All  variable  sizes  are  spe¬ 
cified  in  Appendix  A. 

As  will  be  seen  below,  the  Fourier  transforms  are  multifaceted,  each  face 
with  it's  own  subset  of  items  within  a  namelist.  The  namelists  follow  patterns 
which  allows  the  user  to  reduce  his  inventory  of  variables.  This  feature  will 
often  improve  the  user's  computational  efficiency. 

Also  stored  in  the  records  are  normalized  Probability  Density  Functions 
(PDF's).  These  arrays  have  the  obvious  property  that  the  whole  PDF  must  be 
used.  There  are  thirteen  one-dimensional  PDF's  normalized  to  unit  area,  and  two 
joint  PDF's  normalized  to  unit  volume. 

Where  complex  quantities  are  stored,  we  have  elected  to  separate  the  real 
and  imaginary  parts  because  often  an  investigator  is  interested  in  one,  and  not 
the  other.  This  was  also  computationally  efficient  by  eliminating  the  overhead 
caused  by  computing  complex  quantities.  Thus  vectors  which  represent  complex 
quantities  are  stored  as  the  Re,  then  the  Im  parts  separately,  and  not  as  vec¬ 
tors  of  complex  pairs. 

We  have  included  for  convenience  an  Integer  function  which  returns  a  value 
that  defines  the  most  efficient  size  array  to  declare  for  reading  records  from 
the  database.  Using  this  feature  will  optimize  your  program's  I/O. 

At  the  simplest  level,  a  single  subroutine  call  will  retrieve  all  the 
requested  data  and  put  them  into  a  smaller  file.  The  user  needs  only  to  specify 
which  items  he  wants  transferred,  and  those  items  are  placed  in  records  in  the 
order  specified.  There  is  one  output  record  for  data  base  record  which  meets 
the  selection  criteria. 


More  sophisticated  use  of  the  data  base  interrogation  software  requires 
allocation  of  working  storage.  Since  our  design  gets  data  by  names  according  to 
some  rule,  there  are  seven  required  arrays  and  a  name  buffer  of  user  particular 
size. 


II. 2  The  GIT  mm  Wave  Propagation  Experiment 

The  GIT  millimeter-wave  data  consist  of  28  channels  recorded  at  100  Hz.  (In 
this  report  "millimeter-wave"  is  abbreviated  as  mm-wave.)  Figure  4  shows  the 
arrangement  of  receiver  antennas  as  viewed  from  the  transmitter;  the  antennas 
are  designated  by  a  number.  There  are  six  possible  horizontally  separated 
antenna  pairs,  each  with  a  different  spacing.  The  antenna  pair  numbers  and 


Figure  4.  The  instrument  positions  at  the  experiment  site.  The  dashed  and 
dotted  line  denotes  mm-wave  propagation  path  (1.4  km);  the  long-dashed  lines, 
the  optical  propagation  paths  (1  km  each);  and  the  short-dashed  lines,  the  opti¬ 
cal  rain  gauge  paths  (50  m  each).  Solid  lines  show  the  flow  of  micrometeorolo- 
gical  data  to  the  data  acquisition  system  in  the  receiver  trailer.  Antennas  are 
numbered  1  to  5  in  the  receiver  trailer.  The  relative  positions  shown  are  drawn 
approximately  to  scale. 
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their  spacings  are  given  in  Table  V.  For  a  few  experiment  runs,  the  fifth 
antenna,  below  antenna  2  was  used  in  place  of  antenna  4;  the  objective  was  to 
observe  anisotropy  and  multipath  effects.  As  described  in  Ref.  21,  the  data 
yield  the  intensity  at  all  four  antennas,  the  phase  difference  between  all  six 
antenna  pairs,  as  well  as  some  instrumentation  checks  known  as  figures  of  merit 


Table  V.  Antenna  pairs 


Pair 

Antenna 

Spacing 

number 

pair 

in  meters 

1 

1 

& 

2 

1.43 

2 

1 

& 

3 

4.29 

3 

1 

& 

4 

10.00 

4 

2 

6 

3 

2.86 

5 

2 

& 

4 

8.57 

6 

3 

& 

4 

5.71 

The  mm-wave  data  is  complemented  by  data  collected  from  two  micrometeoro- 
loglcal  (MET)  stations.  Each  station  contains  an  ensemble  of  instrumentation 
designed  to  determine  the  microstructure  of  turbulence  in  the  surface  boundary. 

The  NOAA  micrometeorologlcal  data  are  obtained  at  two  towers,  3.6  m  in 
height,  located  adjacent  to  the  propagation  path  at  approximately  one-third  and 
two-thirds  of  the  distance  from  the  mm-wave  receiver  to  the  transmitter.  These 
are  designated  as  stations  1  and  2  respectively,  as  shown  on  Fig.  4.  The 
instrument  arrangement  on  each  tower  in  shown  in  Fig.  5.  The  two  optical 
paths  were  located  50  m  and  65  m  east  of  the  mm-wave  path.  The  micrometeorolo¬ 
glcal  towers  were  80  m  east  of  the  mm-wave  path,  and  the  optical  rain 
gauge/dlsdrometer  paths  were  90  m  east.  Each  tower  has  a  three-axis  sonic  ane¬ 
mometer  with  a  platinum  resistance-wire  thermometer  in  the  center  of  the  sonic 
array,  a  nearby  Lyman-o  hygrometer  with  another  wire  thermometer  within  its  gap 
a  prop-vane  anemometer,  and  a  psychrometer  for  mean  temperature  and  humidity. 
The  wire  thermometers  and  Lyman-hydrometer  were  sampled  at  100  Hz.  The  sonic 


Figure  5.  MET  tower  mlcrometero Logical  instrumentation  at  each  station. 

The  optical  Instruments  associated  with  each  station  are  shown  schematically  in 
Fig.  4.  A  ■  Propeller-vane  anemometer;  B  -  Three-axis  sonic  anemometer  with 
wire  temperature  sensor  in  center;  C  »  Lyman-a  hygrometer  with  wire  temperature 
sensor  in  center;  D  »  Rotatable  boom  with  angle  encoder;  E  ■  Vaisala  humidity 
probe;  F  *  Aspirated  psychrometer ;  and  G  ■  Electronics  cabinet  containing  power 
conditioning  and  data  telemetry  systems. 

and  prop-vane  anemometers  were  sampled  at  25  Hz.  Beginning  with  the  0.39-Hz 
pulse,  the  psychrometer  gives  six  commutation  cycles  repeated  at  the  25-Hz  rate. 
A  commutation  cycle  consists  of  data  from  the  wet  and  dry  bulbs  and  two 
reference  values,  four  quantities  in  all.  The  psychrometer  then  returns  toward 
thermal  equilibrium  until  the  next  0.39-Hz  pulse  calls  for  data;  this  minimizes 
self-heatlng  in  the  psychrometer  due  to  power  dissipation.  Station  2  has  a 
Vaisala  mean  temperature  and  humidity  instrument  sampled  at  25  Hz,  but  the 
Vaisala  readings  of  humidity  are  judged  inferior  to  those  of  the  psychrometer. 
Station  2  included  an  up-looking  all  sky  pyranometer  for  solar  radiation  and 
cloud-shadow-occurrence  measurements,  and  station  1  included  a  barometer;  these 
were  sampled  at  25  Hz.  A  combined  optical  Cn^  and  optical  cross-wind  meter  was 
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operated  on  the  path  from  near  the  mm-wav,e  transmitter  to  station  1  and  another 
from  near  the  mm-wave  receiver  to  station  2;  these  optical  paths  overlap. 

Optical  rain  gauges  were  deployed  near  both  stations  and  six  weighing  bucket 
rain  gauges  were  deployed  along  the  path.  At  station  1  the  optical  rain  gauge 
included  a  disdrometer.  The  data  from  these  optical  instruments  were  sampled 
every  2.56  seconds.  The  reduction  of  the  data  from  the  rain  instruments  is 
discussed  in  Appendix  III. 

Additional  optical  instrumentation  was  deployed  during  the  fourth  experiment 
session.  A  10.6-pm  CO2  laser,  for  measurement  of  infrared  Cn^  in  clear  air  and 
for  phase-difference  spectra  in  rain,  was  operated  from  near  the  mm-wave 
transmitter  to  inside  the  mm-wave  receiver  trailer.  A  third  large-aperture 
optical  Cn2  and  cross-wind  meter  was  operated  over  the  northern  two-thirds  of 
the  range  directly  below  its  existing  counterpart.  The  purpose  was  to  measure 
vertical  shear  averaged  along  the  50  m  path.  A  He-Ne  laser  scintillometer  and 
LED  scintillometer  were  deployed  over  the  same  short  path  to  measure  the  inner 
scale  of  turbulence  by  its  effect  on  optical  scintillation.  At  station  1,  an 
up-looking  photocell  was  Installed  to  compare  cloud  cover  measurements  with 
those  of  the  pyranometer  at  station  2.  In  addition,  another  psychrometer  was 
installed  at  station  2  to  determine  the  instrument’s  reliability  by  comparison 
with  its  existing  counterpart.  Except  for  the  CO2  laser,  these  additional 
instruments  were  recording  data  during  the  second  half  of  the  fifth  session  as 
well. 

The  heights  of  the  instrumentation  as  well  as  subjective  estimates  of  the 
3o  errors  in  the  heights  are  as  follows:  The  mm-wave  beam  and  C02~laser  beam 
average  3.68  m  ±  0.01  m;  the  Lyman-a  gap,  horizontal  sonic  transducers,  prop- 
vane,  and  platinum  wires  are  at  3.61  m  ±  0.05  m;  the  psychrometer  inlets  are  at 
1.91  m  t  0.05  m;  the  upper-level  optical  wind  and  Cn^  meters  and  the  optical 
inner-scale  meter  are  at  3.78  m  t  0.1  m;  and  the  lower-level  optical  wind  and 
Cn^  meter  are  2m  t  0.02  m  below  the  upper-level  meters. 


-17- 


II. 3  Data  Quantities 


All  data  entered  are  accumulated  over  experiment  times  when  statistics 
are  stationary.  Specifically,  over  times  of  constant  variance  in  most,  or  all 
quantities.  Consequently,  a  given  data  tape  may  contain  several  stationary 
intervals.  In  such  a  case,  a  given  data  tape  may  be  used  to  enter  several 
records  in  the  data  base.  One  may  imagine  a  frontal  passage  during  a  run  where 
the  humidity  changes  drastically.  Figure  6  depicts  this  event,  and  shows  the 
two  resultant  stationary  intervals  which  would  be  entered.  Because  this  study 
is  based  upon  stationary  statistics,  transient  data  are  not  entered  in  the  data 
base.  This  data  is,  however,  available  to  interested  investigators  via  our 
archival  system.  All  data  described  in  this  section,  and  named  in  Appendix  A 
(this  section)  are  entered  as  a  single  record,  one  for  each  stationary  Interval. 


Figure  6.  Example  experiment  with  two  stationary  intervals. 

All  statistical  quantities  are  built  up  from  sets  of  local  time  averages. 
The  interim  files  MILSUMS  and  METSUMS  21  (Fig.  7)  contain  sums,  and  sums- 
of-squares  of  quantities  over  25.6  second  Intervals  over  ‘he  whole  data 
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tape.  They  also  contain  sums  of  parameters  for  the  MIL-  and  MET-  statistics  as 
described  here.  When  these  sums  files  are  reviewed  for  statistical  sta- 
tionarity,  the  intervals  are  selected,  and  the  sums  are  used  in  MILRDR  and 
METRDR  to  accumulate  and  normalize  statistics.  A11  probability  distributions 
are  included  in  this  process.  Both  the  MIL  and  MET  data  used  have  been  repaired 
[Ref.  (23)  and  (24)],  and  data  which  cannot  be  repaired  is  excluded  from  the 
sums  files  so  as  to  not  contaminate  the  statistics. 

The  resultant  statistics  are  incoherent  moments  about  the  mean  for  simple 
variables,  and  normalized  probability  distributions  tallied  over  the  interval  of 
stationarity.  The  interval  of  stationarity  will  also  be  referred  to  as  "the 
interval." 


Figure  7.  Overview  of  the  data  processing  software.. 


Fourier  analysis  of  both  data  sets  is  also  entered,  but  the  processing  is 
necessarily  quite  different.  High  frequency  (hf)  spectra  are  produced  by 
collecting  as  many  4096-point  blocks  of  contiguous  data  as  can  be  within  the 
predefined  stationary  interval.  When  data  are  of  good  quality,  and  the  interval 
of  stationarity  spans  more  than  81920  data,  the  FFT  program  divides  the  interval 
into  20  non-overlapping,  evenly  spaced  regions  from  which  twenty-4096  point 
spectra  are  taken.  This  rule  is  invoked  over  any  subrange  of  data  when  break¬ 
point  pairs  are  included  in  the  interval.  Any  number  of  4096  point  contiguous 
time  series,  up  to  20,  may  be  established  over  the  interval  of  stationarity, 
individually  transformed  and  averaged  together  to  produce  the  hf  power  and  cross 
spectra.  The  spectra  are  then  combined  into  36  bins  for  MET  data,  or  72  bins 
for  MIL  data  over  approximately  logarithmically  spaced  requency  bins.  In  each 
case,  the  first  bin  represents  the  DC  component  of  spectra. 

Low  frequency  (If)  spectra  are  produced  by  block  averaging  the  data  into 
longer  averaging  time  bins  and  transforming  these  bins.  The  length  of  the  bins 
depend  upon  the  length  of  the  interval,  but  have  an  upper  limit  of  0.6  sec  per 
bin  if  the  whole  tape  is  used.  Gaps  in  the  If  data  bins  due  to  break  points, 
f rame-or-block  loss  are  filled  using  a  4-point  LaGrangian  interpolation.  This 
ensures  a  contiguous  If  time  series  for  Spectral  Analysis. 

There  is  only  one  If  transform  over  the  entire  interval  of  stationarity. 

All  If  spectra  have  poor  statistical  significance  until  some  Intermediate  fre¬ 
quency  is  reached.  This  manifests  itself  as  noisy  spectra  in  the  low  frequency 
region  of  the  If  spectra.  Figure  8  shows  this  concept  for  a  MET  coherence 
spectra  where  it  is  most  evident.  Because  of  this  bin-only  averaging,  there  is 
a  frequency  below  which  all  coherence-spectral  components  are  identically  1. 

This  is  designated  fj  in  Fig.  8. 

The  list  of  cross  spectra  for  all  pairs  processed  are  as  follows: 


Co-spectra:  Given  two  quantities  xj  and  X2,  whose- spectral  transforms  are  ri 
and  r2  having  the  spectral  components  at  some  discrete  frequency  are  shown  in 
Fig.  9.  (The  r's  represent  vectors  in  the  complex  plane.) 


Met  Station  1  25H2  Data  FFT 
Lyman-a  Hygrometer  p  ‘Lyman-  a  Wire  T 
Coherence  Spectra  For  Tape  *S4  T16 
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Figure  8.  Temperature-humidity  coherence  spectra.  The  poor  statistics  at 
low  frequencies  is  evident  below  f  ^ . 


Figure  9.  Phase  relationship  between  two  transformed 
quantities  at  a  spectral  component. 
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We  define  the  co-component  as  the  projection  of  r2  on  ri  at  this  frequency. 
The  co-spectra  is  then  the  spectrum  of  all  the  r^^  cos(02  ~  0j)  quantities. 

It  represents  the  fraction  of  energy  in  r2  which  is  in  phase  with  rj. 

Quad-spectra:  The  quadrature  component  is  defined  as  the  projection  of  T2  on  a 
line  perpendicular  to  rj.  It  represents  the  energy  in  r2  which  is  90°  out  of 
phase  with  rj,  and  makes  no  contribution  to  the  energy  density  of  rj.  The 
quadrature  spectra  is  the  spectrum  of  all  the  i  sin(©2  -  0^)  quantities, 

where  i  «■  /-I. 

The  CYBER  tape  is  used  as  input  for  the  MILFFT  processing.  Although  the 
incoherent  averages  have  been  filtered,  the  data  used  in  the  FFT's  are  not. 
Trends  in  the  data  due  to  parametric  drift  in  the  amplifiers,  mixers,  etc.,  con¬ 
taminate  the  If  spectra,  and  the  If  region  of  the  hf  spectra. 

The  MET-tape  is  used  as  input  to  the  METFFT's.  Because  the  MET  data 
problems  are  dealt  with  at  METAPE  generation,  filters  are  not  necessary,  and  the 
METFFT's  are  generally  quite  good. 

The  quantities  calculated  from  the  millimeter-wave  data  and  the  equations 
used  for  the  computations  are  described  in  detail  in  Ref.  21.  Here  we  list  the 
calculated  quantity.  We  calculate  the  phase  difference  variance  and  intensity 
covariance  for  each  of  the  six  antenna  pairs  as  well  as  the  intensity  variance 
for  each  of  the  four  antennas.  We  calculate  the  largest  and  smallest  observed 
phase  difference  for  each  antenna  pair  and  intensity  for  each  antenna.  The  pro¬ 
bability  density  functions  are  calculated  for  each  of  the  six  phase  differences 
and  four  intensities.  The  cross  correlation  between  intensity  and  phase  dif¬ 
ference  is  calculated.  The  second  and  fourth  moments  of  the  electromagnetic 
field  are  determined.  We  calculate  the  power  spectra  and  cross  spectra  of  the 
intensities  and  phase  differences.  Quantities  known  as  figures  of  merit  are 
obtained.  - 

The  micrometeorologlcal  (hereafter  called  MET)  instrumentation  is  designed 
to  collect  data  which  will  provide  statistics  on  clear  air  turbulence.  From 
this  data  we  calculate  structure  functions  of  temperature,  humidity  and  reTr~c- 
tlve  index.  Also  calculated  are  Fourier  spectra  and  cross-spectra  from  data 
taken  at  25  and  100  Hz  sample  rates.  In  support  of  this,  we  calculate  averages, 
.variances,  and  cross-correlation  statistics.  All  MET  data  are  in  mgs  units. 
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During  clear  air  runs  for  which  the  MET  stations  were  fully  operating,  we 
calculate  average  and  rms  values  of  air  temperature  <T>,  humidity  <q>,  solar 
irradiance  <S>,  wind  speeds  <V>,  wind  velocity  components  <u>,  <v>,  <w>,  and 
wind  direction  <0>.  Variances  are  centered  about  the  mean.  Static  pressure 
averages  <P>  in  millibars  are  entered  without  variances.  We  have  several  ways 
of  measuring  temperature  and  humidity.  All  such  determined  values  are  entered 
in  the  data  base  except  the  Valsala  probe  because  of  it's  poor  performance. 

Cross  correlations  between  wind  components,  temperature  fluctuation  and 
humidity,  temperature  and  vertical  wind  (temperature  flux,  and  humidity  flux 
respectively)  are  calculated  about  the  mean.  An  instrument  offset  is  subtracted 
from  the  mean  to  correct  for  misalignment  of  the  sonic  anemometers.  When  the 
mean  wind  directions  in  the  MET-tower  coordinates  (boom  normal)  are  determined, 
wind  averages  are  rotated  into  the  streamwise  components.  Horizontal  and  ver¬ 
tical  momentum  fluxes  are  calculated  for  both  down-  and  cross-stream  components. 
In  addition,  we  calculate  these  fluxes  assuming  a  horizontally  homogeneous  sur¬ 
face  layer.  The  equations  used  for  calculating  the  quantities  are  given  in 
Ref.  21  (Appendix  IV),  as  is  the  definitions  of  geographic  and  instrumentation- 
boom  wind  angles  as  well  as  the  convention  for  the  optical  scintillometer  wind 
readings. 

We  choose  a  rotation  so  that  u@  becomes  the  downstream  component  of  wind, 
and  ve  becomes  the  cross-stream  components.  Selection  of  the  angle  through 
which  wind  related  data  is  rotated  depends  upon  the  performance  of  the  wind 
Instruments.  This  angle  (0)  may  be  the  average  vane  angle,  or  the  angle 
tan-1(<v>/<w>)  from  the  sonlcs.  The  latter  is  the  conventional  method  but  we 
prefer  the  former  for  most  runs  because  ’’he  sonics  offsets  are  not  well  known. 

In  addition,  we  may  select  0  assuming  a  horizontally  homogeneous  boundary  layer; 
0  »  tan”*  (<vw>/<uw>),  or  by  some  externally  derived  method,  0  -  input.  A  flag 
is  entered  in  the  data  base  which  specifies  how  0  was  determined  for  each  sta¬ 
tion.  For  each  station,  the  wind  selection  code,  THETSELECT  values  mean: 

THETSELECT  selection 

1  average  vane  angle 

2  average  sonics  angle 

3  horizontlly  homogeneous  wind  layer 

4  externally  derived  value 
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11.4.  Simple  Macroscopic  variables 


The  angle  encoders  at  each  MET-tower  give  the  geographic  direction  in  which 
boom  is  facing. 

Simple  variables  that  use  sonics  data  are  entered  in  the  data  base  in  both 
the  boom-normal  and  stream-normal  coordinates.  For  each  station,  data  include: 

average  wind  components  <u>,  <v>,  <w>;  <u@>,  <vq> 

variance  of  wind  components  u',  v',  w' ;  u'g,  v'g 

wind  stresses  <uw>,  <vw>;  <ugw>,  <vgw> 

cross  products  <uv>;  <ugvg> 

momentum  flux  <u'v'>;  <u'gv'g> 

The  simple  variables  which  are  directionally  independent  include  humidity 
corrected  dry  air  density,  buoyancy  flux,  radio  frequency  refractive  index,  dry 
bulb  air  temperature,  wet  bulb  temperature,  structure  parameters  of  optical 
refractive  index  Cn2,  optical,  Oj2  derived  from  optical  C^2,  the  boom's 
geographic  orientation,  atmospheric  barometric  pressure,  integrated  all¬ 
sky  solar  lrradlance,  absolute  humidity  calculated  from  both  Lyman-a  averages 
and  the  psychrometer  data,  the  horizontal  wind  speed  (<V>),  and  wind  direction 
(6).  Directionally  independent  cross  correlations  of  simple  variables  include: 

humidity  flux 
temperature  flux 
latent  heat  content 
fine-wire  cross  correlation 
horizontal  wind  stress 
dry  air  buoyant  stability 
moist  air  buoyant  stability 

II. 5.  Turbulence  data 

The  quantity  L  above  is  the  Monin-Obukov  stability  length  and  is  based  upon 
the  downstream  wind  stress  u*3,  known  as  friction  velocity 


L 


(II. 5.1) 
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where  g  -  9.8  m  s~2f  <  >  von  Karman  constant,  T  is  temperature  in  K,  and  only 
the  sonics  temperature  fluctuations,  Ts'  are  used,  we  use  k  =  .4.  We  also 
stored  in  the  file  for  each  MET  station  the  friction  velocity 

u*  -  <uew>1/2  •  (11.5.2) 


and  heat  flux  of  moist  air 
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(II. 5. 3) 


Qv  is  used  in  place  of  the  temperature  flux  in  (II. 5.1)  to  calculate  the  moist 
air  buoyant  stability  parameter,  z/L v.  In  each  case,  the  buoyancy  parameter  is 
also  calculated  and  stored  assuming  a  horizontally  homogeneous  surface  layer 
using  temperature  flux  calculations  as  though  THETSELECT  »  3. 

In  clear  air,  the  dry  bulb  and  wet  bulb  temperatures  at  each  station,  n  ■ 
1,2,  (TDn,  TWn)  are  determined  from  the  psychrometer  data.  Temperature  varian¬ 
ces  are  taken  from  the  fine  wires.  The  humidity  fluctuation  is  from  the 
calibrated  Lyman-a  data.  The  psychrometer  is  used  to  calibrate  the  Lyman-a, 
this  is  high  quality  humidity  data,  and  all  run  calibrations  are  scrutinized 
before  an  entry  is  made.  The  humidity  in  gm~3  is  stored  from  both  the  Lyman-a 
and  psychrometer. 


The  density  of  the  air  is  calculated  using  the  measured  pressure  and  tem¬ 
perature  in  the  ideal  gas  law.  We  apply  a  moisture  correction  term  using 
average  humidity  from  the  psychrometer  as  follows: 


pair 


P 

R  ,  T  . 
air  abs 


0.61 


(II. 5.4) 
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11.6.  Structure  parameters 


The  details  of  the  calculation  of  the  following  structure  parameters  is 
given  in  Appendix  IV.  The  space-lagged  structure  parameter  obtained  from  the 
temperature  difference  between  the  platinum  wire  thermometers  in  the  sonic  ane¬ 
mometer  and  Lyman-a  hygrometer  is  obtained.  Under  an  assumption  of  isotropic 
turbulence  the  space-lagged  temperature-humidity  cross-structure  parameter  is 
also  obtained.  Time-lagged  temperature  and  humidity  structure  parameters  as 
well  as  the  temperature-humidity  cross-structure  parameter  are  calculated  using 
three  different  time  lags  all  of  which  are  in  the  inertial  subrange.  The  radio 
ret - ive-index  structure  parameter  is  calculated  from  the  above  MET  structure 
parameters  using  sensitivity  coefficients  Ax  and  Aq  obtained  from  the  radio 
refractive-index  equation. 

These  data  are  stored  in  the  array  TIMELAGS  which  is  dimensioned.  The  first 
index,  n  -  1,3  refer  to  the  three  time  lags  described  above.  The  second  posi¬ 
tion  index  has  the  following  meaning: 

i^  time  lagged  parameter 

1  Cq2  from  Lyman-a  humidity 

2  Cj2  from  Lyman-a  wire 

3  Cx^  from  sonic  wire 

4  Cxq  from  Lyman-a  humidity 

and  its  wire 

The  actual  time  lags  are  not  interesting  and  are  not  stored. 

Of  the  path  averaged  optical  data  described  in  section  II. 1,  only  the 
optical  cross  wind  and  refractive  index  structure  parameter,  Cn2»  Is  averaged 
over  the  interval  and  entered  in  the  data  base,  we  also  calculate  Cx2  from  the 
optical  Cn2  using  - 


(79  x  io'6  P/T2)2 


(II. 6.1) 


Probability  distribution  functions  (PDF's)  are  calculated  for  the  pyrano- 
meter  solar  flux  (w  m~2)  and  wind  angles  at  each  station  (deg-arc)  from  the 


vanes.  Each  PDF  is  a  50  element  array.  The  central  48  elements  is  the  nor¬ 
malized  PDF,  and  the  1st  and  50th  are  catch  bins  for  values  outside  the  PDF 
limits.  One  must  also  retrieve  the  PDF  limits  from  the  file  so  the  data 
have  meaning. 

The  MET  station  PDF's  are  stored  as  arrays  of  real  numbers  which  represent 
the  probability  that  a  value  x  lies  in  the  range  (x,  x  +  Ax).  The  increment  Ax 
is  calculated  from  the  retrieved  limits  using 

Xmax  Xmin 
4x - 48 - 

The  useful  bins  are  for  i  *  2  to  48  inclusively;  thus  the  ith  bin  represents  the 
probability  that  a  value  x  lies  in  the  range 

[(xmin  +  (1"2)  *Ax)’  <Xmin  +  (1_1)  *Ax)]  * 


Joint  PDF  ( JPDF)  of  Lyman-a  humidity  and  its  fine  wire  temperature  fluc¬ 
tuation  is  also  calculated  for  each  station.  The  structure  of  the  JPDF  is  the 
two-dimensional  analog  of  the  PDF.  The  catch  bins  become  a  square  perimeter 
around  the  48  x  48  cell  JPDF  whose  limits  are  also  stored  on  file.  The  joint 
probability  that  the  temperature  lies  in  the  ith  bin  of  its  range  and  the  humi¬ 
dity  lies  within  the  jth  bln  of  its  range  is  stored  in  the  (l,j)  element  of  the 
JPDF  array. 

For  all  PDF's,  the  cell  widths  are  determined  by  reading  the  limits  and 
dividing  the  spread  by  48. 


i.e.,  AT' 


(I»«  -  W/48 


II. 7.  Fourier  spectra  of  MET  data 

We  calculate  the  power  spectra  of  the  fluctuating  temperature  and  humidity 
from  the  resistance-wires  and  the  Lyman-a  hygrometer,  as  well  as  the 
temperature-humidity  cospectra  and  quadrature  spectra  from  these  instruments. 
We  also  calculate  the  power  spectra  of  all  three  orthogonalized  components  of 
the  fluctuating  wind  from  the  sonic  anemometers,  as  well  as  the  cospectra  and 
quadrature  spectra  of  these  components  taken  in  pairs.  The  cospectra  and 
quadrature  spectra  are  calculated  for  temperature  and  vertical  velocity,  and 
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for  humidity  and  vertical  velocity. 

FFT's  are  combined  to  produce  sets  of  cross  spectra  between  pairs  of  quan¬ 
tities.  The  cross  spectra  produced  include  cospectra,  quadrature  spectra, 
coherence,  phase  spectra,  normalized  co-  and  quadrature-spectra.  Power,  and 
cross  spectra  on  the  interval  of  stationarity  come  in  two  parts.  The  high  fre¬ 
quency  (hf)  spectra  are  transformed  from  a  4096-point  time  series.  This 
corresponds  to  40.96  sec  of  data  if  only  100-Hz  quantities  are  involved  or 
163.84  sec  of  data  of  it  at  least  one  quantity  is  of  the  25-Hz  sampling  rate. 
When  data  are  of  mixed  sample  rate,  the  higher  data  rate  is  block,  averaged  into 
a  time  series  at  the  lower  data  rate  is  averaged  into  a  time  series  at  the  lower 
rate.  The  FFT's  are  then  combined  into  36  bins  and  entered  in  the  file. 

In  addition  to  the  cross  spectra  described  in  section  11.3  the  MET  spectra 
includes : 

Normalized  co-,  and  normalized  quad-spectra:  The  co-  and  quadrature-spectra 
described  in  11.3  are  normalized  and  stored. 

Phase  spectra:  We  also  calculate  the  spectrum  of  phase  angles  between  the 
Fourier  components  (see  Fig.  9): 

(02  -  9|) 


Coherence  spectra:  We  define  the  coherence  spectra  as  the  spectra  of 
<rlr2>/<ri><r2>.  described  in  Fig.  9,  where  averages  are  taken  over  contribu¬ 
tions  in  each  of  the  log-averaged  bins  and  the  r's  represent  vectors  in  the 
complex  plane.  It  offers  a  measure  of  how  well  frequency  bins  were  chosen,  and 
how  well  the  bln  average  spectra  mimics  the  full  transform.  The  low  frequency 
components  .of  coherence  are  reduced  to  the  trivial  identity,  1,  because  only  one 
component  of  the  original  t.ransiorm  of  x^  contributed  to  those  If  bins* 

For  each  MET  station,  n  ■  1,2,  both  25  MHz  and  100  MHz  spectra  are  calcu¬ 
lated.  The  pairs  of  MET  quantities  use  in  cross  spectra  are: 
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Because  the  fourier  transform  produces  a  complete  set  of  basis  vectors,  we  leave 
rotations  and  other  orthogonal  transforms  on  the  FFT's  to  the  user. 

During  inclement  weather  experiments,  the  fine  wires  and  Lyman-a  hygrometers 
are  covered.  Rain  drops  in  the  volume  sampled  by  Sonic  anemometers  cause 
spurious  signals  not  related  to  wind  speeds  to  contribute  to  measurement  errors 
there.  Furthermore,  in  session  V,  ice  buildup  was  observed  at  the  inlet  of  the 
psychrometer  at  station  2  which  may  have  changed  its  psychrometrlc  constant  used 
in  data  reduction,  and  on  several  runs,  the  Prop-vane  had  froze  producing  no 
data  at  all.  These  tapes  were  not  used  to  make  entries  in  the  data  base,  but  it 
alerts  us  to  problems  with  the  data  that  rain  and  snow  can  cause.  We  advise 
users  to  be  cautious  when  using  MET  quantities  for  other  than  clear-air  data. 

We  do  not  make  the  METAPE  for  inclement  weather.  All  MET  data  is  derived 
from  the  2.56s  data  base  for  these  tapes.  The  primary  purpose  of  Inclement 
weather  MET  data  entries  is  to  support  the  interesting  mm-wave  data. 

For  MET  station  n,  n  -  1,2,  the  horizontal  components  of  wind  are  derived 
from  the  prop-vane  using  <un>  -  <V  cos  0>  and  <vn>  ■  <V  sin  0>  where  V  is  the 
prop  speed,  and  0  is  the  vane  angle  returned  from  the  2.56s  MET  data  base.  We 
set  <wn>  -  0  . 

For  completeness,  we  rotate  the  wind  components  as  described  in  section 
11.3.  Clearly  only  an  externally  derived  value  for  <0>,  or  the  average  vane 
angles  can  be  used  to  select  the  rotation  angles. 

We  include  for  comparative  purposes  the  average  cross-path  wind  speed 
measured  by  the  optical  cross  wind  instruments.  The  two  Instruments  look 


towards  each  other  so  they  must  have  opposite  signs  if  there  is  agreement. 

We  suggest  that  the  user  resolve  the  prop-vane  data  into  cross-path  wind 
components  and  compare  the  answers  to  the  optical  data.  This  will  yield  a 
measure  of  data  reliability.  We  have  observed  the  following: 

a.  In  heavy  rain,  dense  fog,  or  snow,  the  optical  cross  wind  data 

can  get  contaminated  from  scintillation  due  to  scattering  by  the  par¬ 
ticles;  or  attenuated  due  to  absorption.  In  this  case  the  prop-vanes 
are  probably  more  reliable. 

b.  In  case  of  extreme  cold  and  high  humidity,  the  prop-vane  may  have  accu¬ 
mulated  ice  causing  errors  in  both  wind  speed  and  direction.  If  the 
precipitation  is  not  too  heavy,  the  optical  data  is  probably  better, 
but  only  one  component  of  wind  is  known;  and  not  well  at  that. 

c.  In  all  cases  of  wind  below  about  2  m  s~l,  the  vane  angle  is  usually 
wrong  due  to  the  long  lag  in  response  to  directional  changes.  When  the 
wind  calms  to  below  1  m  s“l  the  prop  becomes  suspect.  (In  clear-air 
data  we  rely  on  the  sonics). 

The  psychrometer  is  used  to  determine  the  humidity,  dry  bulb  temperature  and 
wet  bulb  temperature.  This  data  is  written  on  the  file  where  Lyman-a  hygrometer 
data  would  normally  be.  Temperature  variances  are  based  upon  the  dry  bulb  tem¬ 
perature  from  the  psychrometer.  Atmospheric  pressure  is  derived  from  barometer 
as  usual.  The  moisture-corrected  air  density  is  also  calculated  and  entered 
using  water  vapor  density  from  the  psychrometer  and  the  ideal  gas  law.  A  point 
of  caution  regarding  the  pyranometer  data:  Rain  and  snow  tend  to  cool  this 
detector  so  that  the  reported  solar  lrradiance  is  generally  low. 

Perhaps-  the  most  interesting  data  to  correlate  with  the  mm-wave  intensity 
fluctuations  is  rain  rate  And  'drop  size  distribution.  Because  of  the  special 
characteristics  of  the  optical  rain  gauges,  and  need  for  per-event  calibration. 

The  optical  rain  gauge  output  is  in  volts-dc,  and  disdrometer  average  counts  are 

25 

stored  unnormalized.  A  separate  publication  is  being  prepared  which  will 

unravel  the  mystery  of  the  optical  rain  gage  as  we  produce  calibration  data  for 
the  several  runs  of  interest. 
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The  disdrometer  data  is  somewhat  simpler  to  reduce,  but  determination  of 
drop  size  distribution  is  left  to  the  user  because  of  several  models  in  the 
literature.  The  data  in  the  array  DROPSIZES  are  129  integers  representing  the 
correlation  for  slit  crossing  times.  The  first  65  cells  represents  a  correla¬ 
tion  time  from  0.0  to  +0.64  s  by  0.0*  seconds,  and  the  next  64  cells  represent 
correlation  times  from  -0.64  s  to  -0.01  s  by  0.01  seconds.  The  number  in  each 
cell,  Nc,  is  related  to  the  correlation  of  the  average  rain  drop  passing  between 
slits  with  a  correlation  time  tc  defined  above.  The  normalized  correlation 
coef f icient(s)  directly  from  the  correlator  is  given  by 

N 

rt  ■  8lnti  sifi*  (n.7.1) 

The  slit  crossing  times  are  related  to  rain  drop  vertical  velocity,  and 
therefore  size;  but  just  how  depends  upon  the  user's  choice  of  a  model  of  a 
droplet  falling  at  terminal  velocity,  i.e.,  is  it  spherical,  or  does  it  have  a 
flattened  bottom?  We  present  the  correlation  function  only.  Use  of  this  data 
*3  user's  discretion.  Users  of  this  data  are  encouraged  to  review  articles 
referenced  in  Appendix  III.  Figure  10  depicts  the  data  stored  in  the  array 
DROPSIZES  after  folding  the  correlation  time  scale  and  normalizing  using 
(II. 7.1). 


CORRELATION  TinE  IN  Sc-OnCS 

Figure  10.  Disdrometer  output  folded  back  and  normalized 
using  Eq.  II-7.1. 


II. 8.  Accessing  the  Data 


There  Is  complete  freedom  for  the  user  to  build  the  subroutines  into  a 
program  which  uses  the  retrieved  data,  or  to  retrieve  data  and  transfer  it  to  a 
smaller  local  file  for  later  investigation.  The  latter  is  the  simplest  use  of 
the  file,  and  the  steps  are: 

1.  Establish  the  rule  by  which  data  is  to  be  retrieved. 

2.  Define  which  data  to  retrieve. 

3.  Get  the  data. 

We  refer  to  the  data  base  interrogation  software  from  here  on  as  "the 
system.” 

II8.1  Establishing  the  keys  of  the  records:  making  the  rules 

A  rule  by  which  records  are  read,  (or  ignored),  is  established  by  setting 
keys  in  seven  categories.  Table  VI  shows  the  categories  and  the  maximum 
number  of  possible  values  that  can  be  set  in  each  category.  Each  category  is  an 
integer  variable,  or  array  of  size  up  to  the  limit  shown.  Records  may  be 
collected  in  any  combination  of  categories,  and  any  combination  of  categorical 
values.  An  example  follows  where  two  weather  conditions  are  selected. 

One  must  also  tell  the  system  how  many  keys  are  set  in  each  category.  If 
the  user  specifies  zero  keys  selected,  the  system  does  not  use  that  category  to 
select  records  (in  other  words,  it  will  return  records  having  all  values  of  keys 
in  that  category)  a  bug,  therefore,  is  inherent  in  the  system.  If  a  user 
selects  1  range  of  values  in  a  category,  and  sets  the  key  value  of  that  category 
to  zero,  the  system  will  look  for  records  which  have  a  zero  value  in  that  cate¬ 
gory.  Except  for  the  stationary  interval  category  described  above,  no  records 
will  be  found.  Section  II. 9  describes  a  quick  and  easy  way  of  getting  data 
where  we  have  dealt  with  this  "bug  in  the  system. 


Table  VI.  Record  Key  Categories 


Key 

Max  length 

Explanation 

Weather 

8 

Weather  condition  during  experiment 

Cn2-radio 

8 

A  range  of  radio  frequency 
refractive  index  structure  param¬ 
eters 

Frequency 

8 

A  range  of  mm-wave  transmitter 
frequencies 

Z/L 

8 

A  range  of  bouyant  stability  values 
in  the  surface  layer 

Session 

8 

One  of  five  experimental  sessions 

Tape 

100 

One  of  hundred  of  experiment  data 
tapes  taken  in  the  5  sessions 

I  ST  A 

8 

One  of  several  stationary  intervals 
possible  on  each  tape 

II. 8.  2  Weather  condition 

The  weather  categories  for  the  weather  key  values  are  in  Table  VII.  If  a 
person  wished  to  collect  data  for  all  kinds  of  snow  data  only,  he  would  select 
weather  codes  to  be  4  and  5.  An  example  of  code  follows: 


INTEGER  WX(2) 

DATA  WX/4,5/,  NWX/2/ 

the  variable  NWX  tells  the  system  that  two  weather  keys  are  selected  and  to 
retrieve  data  from  both  keys  specified  in  the  array  WX.  To  collect  data  in  all 
weather  conditions  regardless  of  the  values  in  WX,  the  user  will  set  NWX  *  0. 


Table  VII.  Weather  Condition  Codes 


Key 

Value 

Weather 

Condition 

1 

Clear  Air 

2 

Fog  &  Rain 

3 

Rain 

4 

Snow 

5 

Blowing  Snow 

n 

h 

n 

;  II. 8. 3  Radio  frequency  refractive  index  structure  keys 

Table  VIII  shows  the  range  of  rf  refractive  index  structure  parameters 
u  selected  by  it's  key.  If  a  user  wants  to  collect  data  for  all  Cn^-tad  > 

h  3.16  lO-*1  they  would  code  the  following: 

E 


INTEGER  CNZRAD  (3) 

DATA  CNZRAD/ 1,2, 3/,  NCNZR/3/ 


The  system  will  return  all  records  for  which  the  radio  frequency  refractive 
index  structure  parameter  is  greater  than  3.16  *  ICT11. 


Table  VIII.  Cn2-rad  Nominal  Range  Key 


Key 

Value  Range  of  Cn2-rad 

C^-rad  >  3.16  x  LCT10 

3.16  x  10-H  <  Cn2-rad  <  3.16  x  10*1° 

3.16  x  10-12  <  Cn2-^  <  3.16  x  10-11 

3.16  x  10-I3  <  Cn2-rad  <  3.16  x  10“ 1 2 

3.16  x  io-14  <  Cn2-rad  <  3.16  x  10"13 

3.16  x  10"15  <  Cn2-rad  <  3.16  x  10“ 

Cn2-rad  <  3.16  x  10“ I4 

To  select  all  values  of  Cn2-rad,  he  would  set  NCNZR  =  0.  The  ranges  and  values 
returned  are  clarified  in  Appendix  B. 

II. 8. 4  mm-wave  transmitter  frequency  range  keys 

Table  IX  gives  the  mm-wave  frequency  range,  f,  in  GHz  for  the  key  values 
in  the  frequency  category.  All  frequency  ranges  will  be  selected  when  the  user 
requests  no  categories  similar  to  setting  NCNZR  *  0  above. 


1 

2 

3 

4 

5 

6 
7 
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Table  IX.  mm-Wave  Frequency  Range  Keys 


Key 

Value 

mm-Wave 

frequency  range,  GHz 

Propagated  mm-wave  fre¬ 
quencies  at  Flatville 

1 

f  <  122 

116.30,  118.75 

120.75 

2 

122  <  f  <  157 

142.00 

3 

157  <  f  >  200 

172.91,  173.00,  172.62 

4 

f  >  200 

229.25,  230.44 

II. 8. 5  Buoyant  stability  condition,  Z/L,  keys 

The  parameter,  Z/L,  is  the  ratio  of  the  height  above  ground  of  the  measure- 
meat  (3.68m)  to  the  Monin-Oboukhov  scale  length,  L.  When  the  surface  layer  is 
unstable  (Z/L  <  0),  the  turbulence  is  driven  mainly  by  bouyant  plumes  rising, 
these  plumes  are  driven  by  upward  heat  flux  from  the  surface.  This  condition 
prevails  in  sunny  (afternoon)  weather.  The  ranges  of  this  stability  parameter 
are  given  in  Table  X.  As  in  the  previous  key,  if  no  ranges  are  selected, 
all  Z/L  key  values  will  be  returned. 


Table  X.  Z/L  Range  Key 


Key  Stability 

Value  Parameter  Range  Condition 


1 

2 

3 

4 

5 


Z/L  <  -0.5 
-0.5  <  Z/L  <  -0.1 
-0.1  <  Z/L  <  0.1 

0.1  <  Z/L  <  0.5 

Z/L  >  0.5 


unstable 


near  neutral 


stable 


II. 8. 6  Session,  tape,  and  interval  keys 

Because  occasionally  an  investigator  may  wish  to  access  data  from  a  par¬ 
ticular  experiment,  we  include  these  categories.  Not  all  tapes  are  processed 
through  to  the  data  base  because  of  a  variety  of  problems  not  unusual  to  an 
experiment  of  this  scope.  As  described  earlier,  one  will  generally  specify  that 
the  number  of  session  id's,  the  number  of  tape  id's,  and  the  number  of  stationary 
interval  id's  are  each  zero  so  the  system  will  not  use  this  unscientific  data  as 
selection  criteria  i.e.,  set  NSES  -  NTPE  »  NSTA  *  0  below. 

For  investigators  who  are  interested  in  one  or  more  particular  runs,  the 
range  of  session  data  and  tape  values  are  given  in  Table  XI.  For  all  sessions 
and  all  tapes,  the  valid  values  of  the  interval  keys  are  0-9  representing  sta¬ 
tionary  Intervals  A-J  respectively.  For  example,  if  one  wished  to  review  the 
data  shown  in  Fig.  6,  they  would  code 

INTEGER  SES,  WSES,  TPE,  WTPE,  ISTA(2),  NSTA 

DATA  SES/4/..NSES/1/,  TPE/ 16/,  NTPE/-1/,  ISTA/0,1/,  NSTA/2/  • 

This  would  tell  the  system  that  one  session,  one  tape  and  2  intervals  are  to  be 
returned:  namely  intervals  A  and  B  from  session  4's  tape  16. 


Table  XI.  Session-Tape  Key  Values:  An  Experiment  Overview 


Session 

Key  value 

Valid  Tape 

Key  values 

Session 

Data 

Season 

1 

8,9,14,15,22-25,27 

15  July  '83  -  20  July  '83 

Summer 

2 

2-45 

12  Nov  '83  -  6  Dec 

co 

00 

Fall 

3 

1-93 

11  Feb  '84  -  6  Mar 

00 

Winter 

4 

3-75 

8  Jun  '84  -  3  July 

00 

Summer 

5 

1-77 

8  Feb  '85  -  18  Mar 

in 

00 

Winter 

11.8.7  Data  names 

There  are  538  named  data  items  in  the  Dictionary.  Each  name  has  up  to  ten 
characters  (1  CYBER  word).  The  list  of  names,  and  their  meaning  is  detailed  in 
Appendix  A. 

11.8.8  Accessing  the  data;  method  I:  The  easy  way 

We  provide  a  single  specific  subroutine  which  will  act  Independently  to  open 
the  data  base,  get  requested  data  by  name  according  to  your  rule  as  described  in 
section  II. 8.1,  transfer  the  data  to  a  smaller  local  file,  then  close  the  data 
base.  We  require  that  each  of  the  specification  categories  be  imensioned  to 
their  maximum  size  shown  in  Table  VI,  and  an  array  of  variable  names  be 
declared  large  enough  to  hold  all  the  data  names  requested.  When  all  elements 
of  a  category  are  zero,  the  system  will  not  use  this  category  in  selection. 

The  Format : 

CALL  READATA  (LUNAF.,.  LUNOUT,  I  TYPE,  WX,-  SES,  ZOVRL,  F,  TAPE, 

CNZRF,  ISRA,  VARRAY,  NITEMS,  NUMID) 

The  inputs  to  READATA  are: 

LUNAF  -  The  DATA  BASE  logical  unit  specifier  as  declared  in  a 
program  card. 

LUNOUT  *  The  logical  unit  specifier  for  the  output  file  created 
by  this  subroutine. 


ITYPE 


The  type  of  data  file  to  create  as  follows: 

1  »  Sequential  unformatted  (Binary) 

2  ■  Sequential  formatted-list  directed 

3  -  Random  access  unformatted 

WX  »  8  element  integer  array  containing  the  weather 

categories  of  interest,  or  all  zeros,  ref.  section 
II. 8. 2. 

SES  -  8  element  integer  array  containing  the  session 

numbers  of  interest,  or  all  zeros,  ref.  section 

11. 8. 6. 

ZOVRL  -  8  element  integer  array  containing  the  stability  keys 

of  interest,  or  all  zeros,  ref.  section  II. 8. 5. 

F  ”8  element  integer  array  containing  all  mm-wave 

frequency  categories,  or  all  zeroes,  ref.  section 
II. 8. 4. 

TAPE  =*  100  element  integer  array  containing  all  tape  numbers 

of  interest,  or  all  zeroes,  ref.  section  II. 8. 6. 

CNZRF  »  8  element  integer  array  containing  all  rf  Cn2  of 

interest,  or  all  zeroes,  ref.  section  II. 8. 3. 

ISTA  -  8  element  integer  array  containing  all  stationary 

intervals  of  interest,  or  all  zeroes,  ref.  section 

11 .8.6. 

VARRAY  -  An  array  of  CHARACTER* 1 0 ,  NITEMS  in  length  containing 
the  names  of  data  as  described  in  the  dictionary  in 
Appendix  A. 

NITEMS  -  The  number  of  data  items  requested,  and  the  length  of 
VARRAY. 


The  output  of  READATA  is: 


NUMID  “  The  number  of  records  found  which  meet  the  criteria 
specified  in  the  category  arrays;  WX,  SES,  ZOVRL,  F, 
TAPE,  CNZRF,  and  ISTA. 


If  the  user  specified  all  zeroes  in  all  category  arrays,  READATA  will  put  the 
data  lists  into  the  output  file  for  every  record  entered  in  the  data  base.  This 
is  the  opposite  of  the  bug  described  in  section  II. 8. 

The  local  file  created  by  READATA  on  unit  LUNOUT  is  called  NEWFIL.  Each 

record  contains  the  requested  data  in  VARRAY  in  the  order  of  items  in  VARRAY. 

There  are  NUMID  records  in  the  file  NEWFIL. 
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We  provide  an  example  of  this  usage  in  Appendix  B 


II. 9  Accessing  the  Data:  Method  II,  the  Programmer's  Way 

This  section  describes  the  intermediate  level  subroutines  available  to 
the  user.  READATA  described  in  section  II. 8. 8  uses  these  routines  to  create 
NEWFIL.  The  advantage  to  using  this  method  is  that  an  interim  file  is  not 
needed,  and  the  user  may  work  with  the  data  in  memory.  The  disadvantage  is 
that  the  data  base  will  normally  remain  connected  to  the  job  and  cannot  be  used 
by  others  until  it  is  released.  We  encourage  users  of  this  method  to  make  a 
local  copy  of  the  data  base,  then  return  the  file  before  using  its  clone. 
Programs  in  Appendix  C  is  an  example  of  this  method  of  data  acquisition. 

We  now  describe  the  requirements  for  using  this  method. 

II. 9.1  Method  II:  required  work  space 

All  required  work  buffers  used  by  the  system  and  not  by  the  user  are  allo¬ 
cated  "behind-the  scenes"  for  this  particular  data  base.  The  user  need  only 
allocate  space  for  his  purposes  as  follows: 


11. 9. 2  Category  key  buffers 

Each  category,  described  in  section  II. 8,  requires  an  integer  space  of 
length  equivalent  to  the  number  of  keys  selected  in  that  category.  If  no  keys 
ate  selected,  then  a  simple  Integer  must  be  specified  for  that  category  anyway. 
The  system  must  also  be  told  how  many  keys  are  selected  in  each  category. 

There  are  seven  variables  or  arrays,  plus  seven  key  counters. 

11. 9. 3  Record  identifier  array 

An  integer  array  of  record  Identifiers  which  meet  the  selection  criteria 
must  be  provided.  The  length  of  this  buffer  depends  upon  how  many  records 
satisfy  the  criteria.  256  is  a  good  number,  but  if  your  program  requires  more, 
then  increase  it.  The  maximum  size  that  will  ever  be  needed  is  equal  to  the 
total  number  of  records  in  the  file.  See  below. 
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II. 9. 4  Record  space 


A  real  array  must  be  allocated  of  sufficient  length  to  read  the  requested 
data.  We  supply  an  integer  function,  IWSRTR,  which  returns  the  most  efficient 
length  of  a  record  buffer  to  specify.  Its  arguments  are  the  variable  names 
list,  and  the  number  of  items  requested,  i.e., 

ISIZE  -  IWSRTR( VARRAY, NITEMS) . 

11. 9. 5  Variable  space 

There  are  three  arrays  required  to  access  the  data.  These  array  do  not 
store  the  data.  For  the  sake  of  discussion,  suppose  the  user  is  accessing 
NITEMS  from  the  data  base  by  name. 

11. 9.6  Variable  names  list 

A  character  *  10  array  of  length  NITEMS  must  be  allocated.  This  array  is 
filed  with  valid  data  names.  The  names  of  data  are  given  in  the  Dictionary  in 
Appendix  A. 

11. 9. 7  Variable  locations  list 

An  integer  array  of  length  NITEMS  must  be  allocated  to  store  the  address  of 
the  named  variables.  Subroutine  GETVAR,  described  below  will  return  the  loca¬ 
tions  so  the  user  need  not  do  this. 

11. 9. 8  Variable  sizes  list 

An  Integer  array  of  length  NITEMS  must  be  allocated  to  store  the  lengths  of 
the  variable  named.  Subroutine  GETVAR  will  return  these  lengths. 


II. 9. 9  Data  space 


A  list  of  variables  or  array  which  store  the  named  data  must  be  allocated 
with  sufficient  space  to  store  the  data  returned.  The  variable  sizes  are  listed 
in  the  Dictionary  in  Appendix  A.  No  error  will  result  if  the  lengths  are  wrong. 
The  user  will  merely  get  unpredictable  results. 


11. 9. 10  Subroutine  calls 

We  now  describe  the  subroutines  used  for  acquiring  data  from  the  file  by 
Method  II.  They  are  presented  in  the  usual  calling  sequence. 

11. 9. 11  OPENDB:  Open  the  data  base 

The  Format  is: 

CALL  OPENDB  (LUNAF) 

where  LUNAF  is  the  logical  unit  number  assigned  by  the  user  to  the  data  base. 

11. 9. 12  GETVAR:  Check  the  names  list 

This  program  checks  the  named  items  in  character  *  10  array  items  (NITEM)  to 
see  if  they  are  valid  dictionary  names.  If  they  are,  it  returns  their  locations 
and  sizes.  If  one  item  is  not  in  the  dictionary,  the  program  stops. 

The  Format  is : 

CALL  GETVAR  (ITEMS,  NITEMS,  LOCATN ,  NSIZES) 

Inputs : 

ITEMS  «  Character  *  10  array.  Section  II. 9. 6. 

NITMES  =  Number  of  items  requested,  Section  II. 9.5. 

Outputs : 

LOCATN  -  Integer  array  of  locations,  Section  II. 9. 7. 

NSIZES  -  Integer  s.rray.  of  sizes.  Section  II. 9. 8. 

There  is  a  1:1:1  correspondence  in  the  element  locations  of  the  data  names 


in  ITEMS  to  their  locations  and  sizes 


II. 9. 13  IDS:  Return  record  identifiers  which  meet  selection  criteria 


This  subroutine  assembles  a  list  of  data  base  records  which  meet  the 
selection  criteria  rule(s)  described  in  section  II. 8. 

Format : 

CALL  IDS  (IWX,  NWX,  ISESN,  NSESN,  ISTAB,  NSTAB,  IFREQ,  NFREQ, 

1  ITAPES,  NTAPES,  ICN2RF,  NCN2RF,  ISTAT,  NSTAT,  ID ARRAY, 

1  NUMIDS) 

■Inputs : 

IWX  -  Integer  array  of  weather  keys,  section  II. 8. 2. 

NWX  -  The  number  of  weather  keys  selected.  If  zero,  all 

weather  cortditions  are  selected. 

ISESN  -  Integer  list  of  sessions,  section  11.8.6. 

NSESN  -  The  number  of  session  keys  selected.  If  zero  all 
sessions  will  be  selected. 

ISTAB  -  Integer  array  of  Bouyant  Stability  keys,  section 

11. 8. 5. 

NSTAB  *  The  number  of  Z/L  keys  selected.  If  zero,  all 
stability  conditions  will  be  selected. 

IFREQ  -  Integer  array  of  mm-wave  transmitter  frequencies  range 
keys,  section  II. 8. 4. 

NFREQ  »  The  number  of  frequency  range  keys  selected.  If  zero, 
all  frequency  ranges  will  be  selected. 

ITAPES  ■  Integer  list  of  session  tape  numbers  to  select,  section 

11. 8. 6. 

The  number  of  session  tapes  selected.  If  zero,  all 
tapes  will  be  selected. 

Integer  array  of  radio-frequency  refractive-index 
structure  parameter  range  keys  selected,  section 
II. 8. 3. 

2 

The  number  of  C  keys  selected.  If  zero,  all 
ranges  will  be  selected. 

Integer  list  of  stationary  intervals  to  select,  section 

II. 8. 6. 

NSTAT  ■  The  number  of  stationary  intervals  selected.  If  zero, 
all  intervals  will  be  selected. 


NTAPES  - 

ICN2RF  - 

NCN2RF  - 

ISTAT 


Outputs : 

IDARRAY  ■  An  integer  list  of  records  which  meet  the  combined 
selection  criteria  on  input. 

NUMID  -  The  number  of  records  entered  in  IDARRAY. 

A  potential  pitfall  is  if  the  user  sizes  IDARRAY  to  be  less  than  NUMID 
returned  by  IDS.  If  this  occurs,  adjust  the  size  of  IDARRAY  to  be  greater  than 
or  equal  to  NUMID,  and  rerun  your  job. 

11. 9. 14  READDB:  Read  a  data  base  record 

This  routine  reads  a  single  record  from  the  data  base. 

Format : 

CALL  READDB  (LUNAF ,  RECORD,  LENREC,  ID) 

Inputs : 

LUNAF  -  Logical  unit  number  assigned  to  the  data  base. 

LENREC  -  Length  of  the  record  to  read.  This  is  usually  the 
value  returned  from  IWSRTR,  section  11.9.4. 

ID  -  Record  identifier  to  read.  This  will  normally  be  an 

element  of  IDARRAY  in  section  II. 9. 12. 

Outputs; 

RECORD  -  An  array  of  length  LENREC,  and  mixed  type. 

11. 9. 15  GETVAL:  Retrieve  data  from  the  data  record 

This  program  gets  the  requested  data  from  a  data  base  record. 

Format ; 

CALL  GETVAL  (RECORD,  LENREC,  LOCAL,  ISIZE,  VALUE) 

Inputs : 

RECORD  »  Real  array  of  length  LENREC  returned  from  READDB. 

LENREC  -  Length  of  record  above.  This  is  usually  the  value 
returned  from  IWSRTR,  section  II. 9. 4. 
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LOCAL  =  Location  of  the  variable  being  returned.  This  is 

normally  an  element  of  the  array  LOCATN  returned  from 
GETVAR,  section  II. 9. 12. 

ISIZE  *  Length  of  the  named  variable  being  returned.  This  is 
normally  an  element  of  the  array  NSIZES  returned  in 
GETVAR,  section  II. 9. 12. 

Output: 

VALUE  =  Real  array  of  size  ISIZE  containing  the  returned  data. 


II. 9. 16  CLOSEDB:  Close  the  file  and  release  the  resource 

This  program  closes  the  data  base  and  disconnects  it  from  the  user's  program. 

Format : 

CALL  CLOSEDB  (LUNAF) 

Input : 

LUNAF  3  Logical  unit  number  assigned  to  the  file. 
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APPENDIX  A 


THE  DICTIONARY  OF  DATA 

Table  A-l  is  a  listing  of  the  dictionary  giving  the  variable  names,  posi¬ 
tion,  type,  and  description.  For  brevity  "A"  means  the  average  of  a  value  over 
the  interval.  "S.D."  is  the  standard  deviation  of  a  value  over  the  interval,  or 
square  root  of  its  variance. 

NOTES  ON  THE  LISTING 

Note  1:  References  to  all  intensity  related  arrays  in  the  mm-wave  data  contain 
values  of  intensity  from  antennas  1  through  4  inclusively  in  that 
order.  References  to  6  element  antenna  pair  data  refer  to  pairs  1&2, 
1&3,  1&4,  2&3,  2&4,  and  3&4  for  indexes  =  1  to  6  respectively. 

Note  2:  The  real  and  imaginary  parts  of  M4  in  II. 3  element  indexes  1-33  refer 

to  ijnm  values  for  cases  1-33  given  in  Table  2  of  Appendix  IV  (Ref.  21) 
respectively. 

Note  3:  The  first  6  elements  of  array  FIGAR  contain  the  type  figures  of  merits 
for  antenna  parls  1  through  6.  The  second  6  elements  contain  the 
corresponding  variances. 

Note  4:  The  first  6  elements  of  the  array  FIGVAR  contain  the  E-type  figures  of 
merit  for  antenna  pairs  1-6,  the  second  6  elements  contain  their 
variances.  The  next  4  elements  contain  the  phase-type  figures  of  merit 
for  antenna  triples,  and  the  last  element  contains  the  phase-type 
figure  of  merit  from  4  antenna  pairs.  These  figures  of  merit  are 
explained  in  Appendix  IV  (Ref.  21). 

Note  3:  Generic  FFT  entries.'  Each  FFT  is  proceeded  by  a  processing  header. 

The  spectra  are  entered  as  low  frequency  spectra  ”LF....”,  then  high 
frequency  spectra  "HF...." .  The  first  item  is  the  frequency  bin 
assignments,  next  are  the  power  spectra,  then  the  cross  spectra. 

Note  6:  MIL  FFT  cross  spectra  consist  only  of  co-  and  quadrature  spectra 

designation  "...COSPC”  and  ”....QDSPC”.  The  "Cij"  mean  "cross  spectra 
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between  channel  i  and  j"  where  Table  A-l  defines  the  data  in  "channel  i 
or  j". 


TABLE  A-l. 


i  or  j 

channel  data 

1 

antenna  1  intensity 

2 

antenna  2  intensity 

3 

antenna  3  intensity 

4 

antenna  4  intensity 

5 

phase  difference  pair 

1&2 

6 

phase  difference  pair 

1&3 

7 

phase  difference  pair 

1&4 

8 

phase  difference  pair 

2&3 

9 

phase  difference  pair 

2&4 

0 

phase  difference  pair 

3&4 

Thus,  LFC17QDSPC  at  position  3624  decodes  fo  low  frequency  quadrature 
spectra  between  antenna  1  and  phase  difference  pair  1-4.  This  sequen¬ 
ces  repeats  for  MIL  data  high  frequency  spectra,  ”HF....”. 
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APPENDIX  B 


EXAMPLE  OF  ACQUIRING  DATA  THE  SIMPLE  WAY 

Figure  B-i  is  the  FORTRAN-V  (’77+)  source  listing  of  a  demonstration 
program  which  reads  all  clear-air  data  and  retrieves  the  solar  pyranometer  out¬ 
put,  and  from  each  station;  the  optical  and  rf  refractive-index  structure  para¬ 
meters,  and  the  temperature  structure  parameter  calculated  from  the  optical 
measurements . 

Figure  B-2  is  a  print  of  the  contents  of  the  local  file,  NEWFIL,  created  by 
READATA.  The  first  record  contains  the  number  of  items  requested,  and  the 
number  of  records  found  that  satisfy  selection  criteria. 

As  annotated  in  the  source  listing,  we  collect  records  for  all  clear  air 
tapes  for  which  the  rf  refractive  index  structure  parameter  lies  between  3.12 
1(T14  and  3.16  10~H. 


Figure  B-l.  FORTRAN- ' 7 7  compiler  listing  of  a  program  which  accesses  the  data  by  method  I,  II. 8 
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APPENDIX  C 

EXAMPLE  OF  ACQUIRING  DATA  BY  METHOD  II 

Figure  C-l  and  C-2  is  the  source  listing  of  program  QRATIO  and  its  output 
which  we  used  to  determine  how  well  we  have  calibrated  the  Lyman-oFigure  C-3  is 
a  program  which  we  used  to  see  how  closely  the  intensity  PDF's  match  the  theore¬ 
tical  distribution  of  intensity.  Figure  C-4  is  a  plot  from  the  program  of  Fig. 
C-3. 
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Figure  C-l.  FORTRAN-' 77  source  listing  of  a  program  to  validate  the  MET  sta¬ 
tions  Lyman-a  hygrometer  data. 


THE  NO.  Of  CLEAR  AIR  TAPE  ENTRIES  IS  21 
ALL  VARIABLES  REQUESTED  WERE  FOUND 
FROM  TAPE  *S4  T9 

FROM  STATION  1  Q • PSYCH, Q- LYMAN  A,  PSYC/LYA  =  11.91510  12.74966  .93454 

FROM  STATION  2  Q-PSYCH,Q-LYMAN  A,  PSYC/LYA  =  11.53722  12.94938  .89095 

FROM  TAPE  *S4  T25 

FROM  STATION  1  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  =  9.91492  9.79565  1.01218 

FROM  STATION  2  Q • PSYCH, Q- LYMAN  A,  PSYC/LYA  =  9.61510  9.76911  .98424 

FROM  TAPE  *S4  T3 

FROM  STATION  1  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  =  10.08415  10.10029  .  99840 

FROM  STATION  2  Q-PSYCH.Q-LYMAN  A,  PSYC/LYA  =  9.82749  10.19999  .96348 

FROM  TAPE  *S4  T4 

FROM  STATION  1  Q - PSYCH, 0- LYMAN  A,  PSYC/LYA  =  10.43429  10.40422  1.00289 
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FROM  STATION  1  Q • PSYCH, Q- LYMAN  A,  PSYC/LYA  =  11.91510  12.39897  .96097 

FROM  STATION  2  Q • PSYCH, Q- LYMAN  A,  PSYC/LYA  =  11.53722  12.74030  .90557 

FROM  TAPE  *S4  T16 

FROM  STATION  1  Q • PSYCH, Q- LYMAN  A,  PSYC/LYA  =  13.22113  13.21760  1.00027 

FROM  STATION  2  Q • PSYCH, Q* LYMAN  A,  PSYC/LYA  =  12.71640  13.17422  .96525 

FROM  TAPE  *S4  T16 

FROM  STATION  1  Q • PSYCH, Q- LYMAN  A,  PSYC/LYA  =  12.55008  12.54014  1.00079 

FROM  STATION  2  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  =  12.19556  12.67743  .96199 

FROM  TAPE  *S4  T16 

FROM  STATION  1  0 -  PSYCH, Q- LYMAN  A,  PSYC/LYA  =  12.51995  12.43589  1.00676 

FROM  STATION  2  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  =  12.41343  12.82425  .96797 

FROM  TAPE  *S4  T18 

FROM  STATION  1  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  =  12.89936  12.91992  .  99841 

FROM  STATION  2  O-PSYCH, O-LYMAN  A,  PSYC/LYA  =  12.56615  13.05664  .96243 

FROM  TAPE  *S4  T22 

FROM  STATION  1  Q -  PSYCH, Q- LYMAN  A,  PSYC/LYA  =  10.27393  10.26783  1.00059 

FROM  STATION  2  Q- PSYCH, Q- LYMAN  A,  PSYC/LYA  =  9.98448  10.29055  .97026 

FROM  TAPE  *S4  T23 

FROM  STATION  1  Q-PSYCH,Q-LYMAN  A,  PSYC/LYA  =  12.07097  12.07119  .99998 

FROM  STATION  2  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  «  11.63918  11.88928  .97896 

FROM  TAPE  *S4  T23 

FROM  STATION  1  Q - PSYCH, Q- LYMAN  A,  PSYC/LYA  =  12.14681  12.10343  1.00358 

FROM  STATION  2  O-PSYCH, Q-LYMAN  A,  PSYC/LYA  *  11.86800  12.07404  .98294 

FROM  TAPE  *S4  T26 

FROM  STATION  1  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  10.28450  10.28912  .99955 

FROM  STATION  2  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  9.75101  10.29356  .94729 

FROM  TAPE  *S4  T26 

FROM  STATION  1  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  10.21738  10.31783  .99026 

FROM  STATION  2  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  9.90139  10.54962  .93855 

FROM  TAPE  *S4  T46 

FROM  STATION  1  O-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  11.83544  11.85860  .99805 

FROM  STATION  2  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  11.60790  11.87200  .97775 

FROM  TAPE  *S4  T46 

FROM  STATION  1  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  12.20823  12.28132  .99405 

FROM  STATION  2  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  11.92418  12.26805  .97197 

FROM  TAPE  *S4  T47 

FROM  STATION  1  Q-PSYCH, Q-LYMAN  A,  PSYC/LYA  =  9.58335  9.59166  .99913 

FROM  STATION  2  Q-PSYCH, O-LYMAN  A,  PSYC/LYA  *  9.30255  9.57059  .97199 

FROM  TAPE  *S4  T31 

-  FROM  STATION  1  Q-PSYCH.Q-LYMAN  A,  PSYC/LYA  =  17.14909  17.12146  1.00161 

FROM  STATION  2  Q-PSYCH.Q-LYMAN  A,  PSYC/LYA  *  17.14756  17.22267  .99564 

FROM  TAPE  *S4  T59 

FROM  STATION  1  Q-PSYCH.Q-LYMAN  A,  PSYC/LYA  *  9.88226  9.89075  .99914 

FROM  STATION  2  Q-PSYCH.Q-LYMAN  A,  PSYC/LYA  *  9.61503  9.89534  .97167 


Figure  02.  Copy  of  printer  output  from  the  program  listing  of  Figure  C— 1 . 
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Figure  C-3.  (continued) 
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Figure  C-3.  (continued) 
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Abstract 

This  system,  used  to  measure  the 
effects  of  atmospheric  turbulence  at 
frequencies  between  116  and  173  GHz , 
comprises  a  front-fed,  off-axis  parabola 
transmitting  antenna  located  1.4  km  from  a 
receiver  array.  The  five  antennas  in  the 
latter  receive  local  oscillator  power  from 
a  phase-locked  klystron  source  via  a  low- 
loss  lens  beam  waveguide  distribution 
system. 


A  key  part  in  the  Georgia  Tech  and 
National  Oceanic  and  Atmospheric 
Administration  study  of  millimeter  wave 
atmospheric  turbulence  effects  (1-31  has 
been  the  development  of  an  appropriate 
antenna  system.  The  equipment  may  be 
transported  to  sites  of  interest  and  has 
been  used  at  a  specially  selected  site  (11 
near  Flatville,  Illinois.  The  use  of  an 
array  of  receivers  allows  phase  differences 
in  the  wavefront,  and  hence  angle  of 
arrival  scintillation,  to  be  measured.  The 
method  used  is  described  elsewhere  in  these 
proceedings  (1-21. 

Tittnaaittci 

Three  objectives  met  in  the  design  of 
the  transmitter  antenna  are:  (1)  it  fills 
the  receiver  array  for  propagation 
distances  of  1  km  or  more;  (2)  it  avoids 
overfilling  the  receiver  array  and  has 
minimal  sidelobes;  and  (3)  it  gives  far 
field  conditions  at  the  receiver.  A  picture 
of  the  design  solution  is  shown  in  Fig.  1. 

The  primary  antenna  is  an  off-axis, 
rectangular  (610  X  152  mm)  section  of  a 
paraboloid  reflector.  As  a  result  of  the 
unequal  height  and  width,  the  antenna 
pattern  is  spread  out  horizontally  and  is 
narrow  vertically,  efficiently  filling  the 
receiver  array  which  is  1.5  ■  in  vertical 
extent  and  10  m  in  horizontal.  The  mirror 
was  copy-turned  from  solid  aluminum. 

The  primary  is  front  fed  by  an  E-plane 
corrugated  rectangular  horn,  without 
obscuration.  The  corrugations  have  been 
included  to  minimize  sidelobes  in  the  feed 
pattern,  and  further  reduction  is  effected 
by  application  of  an  absorbing  foam  to  the 
periphery  of  the  mirror. 
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As  shown  in  Fig.  1,  the  antenna  is 
mounted  on  a  pedestal  which  goes  through 
the  floor  of  the  surrounding  van,  down  to  a 
concrete  pad  which  provides  stability  in 
pointing  the  antenna.  The  van  and  antenna 
have  been  sited  on  a  berm  such  that  the 
transmitter  is  at  the  same  height  above  the 
plain  as  the  receiver  array  (namely  3.5m). 

RBCtiXAI  Array 

Each  receiving  antenna  in  the  array 
comprises  a  conical  corrugated  horn  feed,  a 
305  mm  diameter  polyethylene  lens  and  a 
large,  optical  quality  diagonal  mirror  for 
beam  steering,  as  may  be  seen  in  Fig.  2. 
All  are  mounted  on  a  large  I-beam  that 
rests  on  two  pillars  anchored  to  concrete 
ground  pads  for  stability  and  isolation 
from  the  surrounding  semi-trailer  van. 

The  array  may  be  termed  non-redundant, 
since  any  combination  of  two  antennas, 
gives  a  unique  spacing  for  which  the 
relative  phases  of  the  signals  can  be 
measured.  The  six  intervals  in  the 
horizontal  direction  are  1.43,  2.86,  4.29, 
5.71,  8.57  and  10  meters,  intervals  that 
are  not  only  unique,  but  nearly  uniformly 
distributed.  Preliminary  results  til 
indicate  transverse  coherence  lengths  in 
the  propagated  radiation  that  are  within 
the  dimensions  of  the  array.  Thus  the  array 
size  is  well  matched  to  the  measurement 
requirement.  One  antenna  shown  in  Fig.  2  is 
1.43  m  below  the  horizontal  row  of 
antennas.  This  allows  differences  in 
vertical  and  horizontal  scintillation  to  be 
measured. 
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Provided  that  all  receivers  use  the 
same  LO,  it  is  possible  to  measure 
meaningful  phase  differences.  The  mixers  in 
the  receivers  are  of  the  harmonic  type  so 
the  LO  frequencies  are  in  the  range  58  to 
87  GHz.  The  LO  sources  are  klystrons  with 
available  powers  of  a  few  hundred 
milliwatts.  The  waveguide  size  used  for  LO 
components  throughout  this  range  is  WR-12. 

The  combination  of  the  limited  power 
available  and  the  large  array  size  over 
which  the  L0  signal  must  be  distributee 
(takes  it  important  to  have  high  efficiency 
in  the  LO  transmission  medium.  However,  the 
losses  in  long  runs  of  KR-12  waveguide  may 
l  <•  r'Vkr  t  :  re  t :  •  k  . 


Instead,  a  secies  of  relay  lenses, 
termed  a  lens  beam  waveguide  f 4 J ,  was 
chosen  to  distribute  the  LO  signal,  as 
illustrated  in  Fig.  3.  Twenty-one  lenses, 
152  mm  in  diameter  and  made  of  high-density 
polyethylene  were  used,  but  the  longest  run 
between  the  klystron  and  a  receiver 
involved  fifteen.  Corrugated  conical  horns 
couple  the  LO  signal  into  and  out  of  the 
beam  waveguide  and  produce  a  beam  with  a 
Gaussian-shaped  field  profile  (41. 

The  lenses  have  a  curvature  and 
spacing  (737  mm)  such  that  they  are 
confocal  (41.  The  horn  feeds  were  designed 
such  that  a  beam  waist  occurs,  for  a 
frequency  of  55  GHz,  half  way  between  the 
lenses  as  shown  in  Fig.  3.  At  a  higher 
frequency,  the  feed  pattern  changes  and  the 
beam  waist  shifts,  but  at  every  odd 
numbered  lens  in  a  series,  the  beam  returns 
to  the  size  that  couples  efficiently  to  the 
output  horns.  Since  all  legs  of  the  beam 
waveguide  have  an  odd  number  of  lenses,  it 
works  well  over  a  wide  frequency  range. 

The  splitting  of  LO  power  between  the 
paths  to  the  various  receivers  is 
accomplished  by  stretched  electrof ormed 
nickel  meshes.  Selections  were  made  of 
meshes  of  different  reflectance  to  balance 
the  power  to  all  receivers  at  crucial 
frequencies.  It  is  calculated,  from 
measured  lens  and  mesh  parameters  that 
about  ten  percent  of  the  klystron  power 
launched  into  the  beam  wavequide  is 
delivered  to  each  receiver.  Satisfactory 
mixer  performance  has  been  obtained  with 
this  system. 

Conclusion 

The  first  field  tests,  described 
elsewhere  in  these  proceedings  (11,  have 
demonstrated  the  success  of  the  designs 
presented  in  this  paper.  A  notable 
achievement  is  the  lens  beam  wavequide  LO 
distribution  system  with  its  wide  frequency 
range . 
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Fig.  1  Cutaway  of  transmitter  van 
showing  transmitting  antenna 
mounted  on  a  pedestal. 
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Receiver  array  showing  four 
antennas  above  the  I-beam  support 
and  one  below. 
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Fig.  3  A  representative  portion  of  the 
lens  bean  waveq-ide  j 
distribution  system. 
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Atmospheric  fluctuations  in  temperature,  humidity,  and  total  pressure  cause  fluctuations  in  the 
absorption  and  refraction  of  electromagnetic  waves.  We  consider  only  the  absorption  and  refraction 
fluctuations  caused  by  variations  in  the  width  and  strength  of  water  monomer  resonances.  We 
differentiate  the  absorption  coefficient  and  refractive  index  with  respect  to  temperature,  humidity, 
and  total  pressure  in  order  to  obtain  the  coefficients  that  relate  the  fluctuations  in  absorption 
and  refraction  to  these  atmospheric  parameters.  The  resulting  coefficients  as  well  as  the  absorption 
coefficient  and  the  refractivity  are  calculated  by  summing  over  all  water  vapor  infrared  resonances 
and  are  displayed  graphically.  In  the  limit  of  low  frequencies,  these  calculated  coefficients  are 
in  good  agreement  with  the  formula  for  radio  wave  refractive  index.  The  effects  of  turbulent 


fluctuations  in  total  pressure  are  usually  negligible 
INTRODUCTION 

Fluctuations  of  the  amplitude  and  phase  of 
electromagnetic  waves  propagating  through  the 
atmosphere  are  caused,  in  part,  by  local  fluctuations 
in  absorption  and  refraction.  We  consider  the  ab¬ 
sorption  and  refraction  caused  by  the  water  vapor 
monomer  and  calculate  the  relationship  between 
local  fluctuations  in  absorption  and  refraction  and 
the  local  fluctuations  in  temperature  T,  humidity 
Q,  and  total  pressure  P.  Similar  results  have  been 
determined  for  the  microwave  resonances  of  02 
by  Liebe  and  Hopponen  (1977],  The  combination 
of  their  results  and  ours  gives  a  complete  view 
of  fluctuations  in  microwave  resonance  absorption 
and  refraction  caused  by  fluctuations  in  T,  Q,  and 
P.  Our  results  are  reported  in  detail  here  for  fre¬ 
quencies  less  than  16  cm-1  (480  GHz)  because  of 
their  interest  to  radio  scientists,  but  our  calculations 
are  complete  from  radio  wavelengths  to  nearly  5 
pim  and  are  reported  by  Hill  et  al.  (1980a,  b\. 
In  a  related  study,  Zhevakin  and  Naumov  [1967a] 
calculate  the  derivatives  of  the  real  part  of  the 
refractive  index  with  respect  to  temperature  and 
water  vapor  partial  pressure  for  frequencies  in  the 
range  1  cm*1  <  v  <  20  cm'1,  where  v  =  l/\  is 
inverse  wavelength. 

There  are  mechanisms  other  than  the  H20  and 
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02  resonances  that  absorb  electromagnetic  waves. 
For  instance,  there  is  extinction  by  aerosols  as  well 
as  the  observed  and  as  yet  unexplained  excess 
attenuation.  We  do  not  consider  the  absorption 
fluctuations  caused  by  these  mechanisms.  For  in¬ 
stance,  consideration  of  aerosol  or  particulate  con¬ 
tributions  would  necessitate  considering  fluctua¬ 
tions  in  particle  concentration  and  particle  size 
distribution.  Fluctuations  in  T and  Q  cause  fluctua¬ 
tions  in  aerosol  extinction  by  varying  the  particle 
size  distribution  [Nilsson  1979] .  The  observed  ab¬ 
sorption  (for,  say,  v  s  50  cm"1)  is  greater  than 
that  calculated  from  the  sum  of  all  molecular  reso¬ 
nances.  This  excess  absorption  is  greater  than  that 
from  resonances  only  in  the  windows  between  the 
absorption  resonances.  Including  the  excess  ab¬ 
sorption  in  the  present  analysis  would  require  ac¬ 
curate  knowledge  of  its  variation  with  T,  Q,  P, 
and  v.  Waters  [1976]  describes  the  model  of  excess 
absorption  proposed  by  Gaut  and  Reifenstein  which 
is  proportional  to  humidity,  frequency  squared,  and 
T~* Recent  measurements  have  shown  a  quadrat¬ 
ic  rather  than  linear  dependence  on  humidity  at 
both  microwave  ( Hogg  and  Guiraud,  1979]  and 
higher  frequencies  (for  13.8  cm'1  s  v  s  50  cm"1 
by  Bohlander  [1979]  and  at  v  =  7. 1  cm"1  by 
Llewellyn-Jones  et  al.  [1978]).  Recent  measure¬ 
ments  [ Bohlander ,  1979;  Rice  and  Ade,  1979]  have 
also  shown  a  more  gradual  increase  of  excess 
absorption  with  frequency  for  v  >  7  cm"1  than 
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that  given  by  the  Gaut  and  Reifenstein  model.  The 
observations  by  Bohlander  [1979]  in  the  tempera¬ 
ture  range  295  K ,<T<  355  K  indicate  a  temperature 
dependence  of  excess  absorption  more  nearly  like 
T  4  for  13  cm-1  <  v  <  20  cm-1,  increasing  to 
T~6  for  20  cm'1  <  v  <  35  cm-1.  For  v  =  7.1 
cm-1,  Llewellyn- Jones  et  al.  [1978]  found  excess 
absorption  to  have  a  temperature  dependence  like 
T~ 4  for  T2  300  K  but  a  much  more  rapid  decrease 
with  increasing  temperature  for  T  £  300  K.  We 
conclude  that  the  dependence  of  excess  absorption 
on  v,  T,  Q,  and  P  at  microwave  to  submillimeter 
wavelengths  is  not  sufficiently  well  known  to  in¬ 
clude  in  the  calculations  that  follow.  As  excess 
attenuation  has  a  negligible  impact  on  refraction, 
this  problem  does  not  affect  our  values  of  refraction 
or  the  contributions  of  fluctuations  in  T,  Q,  or  P 
to  refraction  fluctuations.  In  fact,  measurements 
by  Kemp  et  al.  [1978]  of  refraction  by  water  vapor 
in  the  range  10  cm'1  <  v  <  50  cm-1  show  good 
agreement  with  refraction  calculated  from  the  sum 
of  water  monomer  lines  using  the  Gross  line  shape. 


THEORY 

We  write  the  refractivity  of  air  as  its  value  at 
optical  frequencies  Nopt  plus  the  refraction  Na  due 
to  the  infrared  resonances  of  water  vapor: 

N=Na„  +  N„ 

where 


*op,  =  Nd+Nw 
Nd  =  103.29  P/T 
1 V„  =  - 1 .765  X  10-"e 

N.  -  2 

I 

and  where  P  is  total  atmospheric  pressure  in  torr, 
Tis  temperature  in  degrees  Kelvin,  Q  is  water  vapor 
concentration  in  molecules  per  cubic  centimeter, 
and  Nal  is  the  refraction  of  the  1th  water  vapor 
resonance.  The  expressions  for  Nd  and  Nw  are 
derived  from  formulae  given  by  Goody  [1965]. 
Since  P  is  the  total  pressure  (not  the  partial  pressure 
of  dry  air),  Nd  contains  contributions  from  H20 
as  well  as  02,  N2,  and  A Nw  is  just  a  correction 
to  account  for  the  fact  that  at  visible  frequencies 
the  refraction  per  molecule  is  less  for  H20  than 
for  02  or  N2 .  The  continuum  refraction  Nd  +  Nm 


is  caused  by  the  ultraviolet  resonances  of  0:,  N, , 
and  H20  and  is  essentially  independent  of  frequency 
for  \  a  10  ixm.  We  have  neglected  the  near  infrared 
resonances  of  02  and  the  refraction  due  to  C02; 
their  relatively  small  effect  is  discussed  by  Zhevakin 
and  Naumov  [  1967a] .  As  noted  previously,  the 
microwave  resonances  of  02  are  discussed  by  Liebe 
and  Hopponen  [  1977] .  The  real  part  of  the  refractive 
index  n„  is  related  to  refractivity  by  N  =  (nR  - 
1)  x  10‘. 

The  absorption  coefficient  p  of  the  water  vapor 
monomer  resonances  is  found  by  summing  the 
contributions  P,  of  all  lines: 


P-£P. 

i 


The  imaginary  part  of  the  refractive  index  is  related 
to  P  by  n,  =  p/4irv. 

We  use  the  Gross  line  shape  so  that  the  refraction 
and  absorption  of  a  single  resonance  with  line 
frequency  v ,  is  given  by 


N., 


S,Q  v,  - 

2irJ  (v?  -  yJ)2  + 


X  10‘ 


S,Q  4v2a, 
it  (wj  -  vJ)J  +  4w2a,1 


where  S,  is  the  line  strength,  and  a,  is  the  line 
width,  lie  line  strength  depends  on  T,  and  line 
width  depends  on  T,  P,  and  Q.  The  functional 
dependence  that  we  use  is  given  by  Hill  et  al. 
[1980a]. 

We  write  T.  Q,  and  P  as  the  sum  of  their  mean 
values  ( T } ,  (Q  ),(?)>  and  their  fluctuations  caused 
by  turbulence  ST,  8(7, 81®.  For  instance,  T  -  { T )  + 
8  T.  We  can  now  relate  the  fluctuations  in  n„  and 
n,  to  those  of  T,  Q,  and  P: 


bn„ 

bn, 


bT  bP 

—  At - +  A  - +  A„ 

<r>  (P) 


bT  bP 

B  - +  Br - +  Bq 

< T )  (P)  0 


bQ_ 

<Q> 

sg 

<g> 


(1) 

(2) 


We  call  the  coefficients  Ar,  Ar,  AQ,  Br,  Bg,  and 
Bf  the  sensitivity  coefficients,  because  they  are 
measures  of  the  change  in  the  refraction  and  ab¬ 
sorption  caused  by  fractional  fluctuations  in  T,  Q, 
and  P.  This  nomenclature  is  similar  to  that  of  Liebe 
and  Hopponen  [1977],  The  dimensionless  sensitiv¬ 
ity  coefficients  are  given  by 
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f)  xl° 

1  /  P.Q 

(dl 

A,=  </»>(  - 
\di 

/ 

x  10 

P/r.e 

A°={Q\1 

N  \ 

-  x  10 

Q/er 

Bt  =  (4irv)-1  (T) 


Br=  (4irv)',</>> 


(-) 

Va7'/ pg 

\dP  )  TQ 

(-) 

/  T.P 


Bq  =  (At<v)-'(Q) 


where  the  partial  derivatives  are  to  be  evaluated 
at  the  mean  values  (T),  (Q).  and  (P),  and  the 
subscripts  denote  the  variables  held  fixed  during 
the  partial  differentiation.  The  analytic  form  of  the 
derivatives  is  given  by  Hill  et  al.  [1980a,  h]  and 
need  not  be  repeated  here.  For  instance,  for  AQ 
we  differentiate  the  explicit  Q  dependence  in  Nw 
and  Nal  as  well  as  an  implicit  Q  dependence  in 
«,(T,  P,  Q). 


CALCULATED  SENSITIVITY  COEFFICIENTS 

There  are  over  42,000  infrared  water  vapor  reso¬ 
nances  in  the  most  recent  Air  Force  Geophysics 
Laboratory  (AFGL)  line  parameters  compilation 
(we  use  the  compilation  discussed  by  Rothman 
[1978]).  For  the  frequencies  discussed  here,  name¬ 
ly,  v<  16  cm'1,  only  the  resonances  within  the 
rotational  band  at  v,  s  330  cm  ~ 1  need  be  considered. 
The  negligibility  of  the  resonances  in  the  vibration¬ 
al-rotation  bands,  v,  a  1000  cm"1,  is  discussed  by 
Zhevakin  and  Naumov  [1967a].  Nevertheless,  for 
each  frequency  v  at  which  we  compute  0 ,  Na,  A r, 
Ar,  Ag,  Br,  Bf,  and  BQ,  we  have  performed  the 
sum  over  all  the  resonances. 

We  calculate  the  sensitivity  coefficients  for  the 
mean  conditions  IT)  =  300  K,  (P)  =  760  torr,  and 
80%  relative  humidity,  so  (Q)  =  6.82  x  1017  mol¬ 
ecules  per  cubic  centimeter  (or  28.2  mbar  partial 
pressure  of  water  vapor).  The  results  are  shown 
graphically  for  v  <  16  cm'1  (480  GHz)  in  Figure 
1.  Separate  graphs  of  Na,  Ar,  and  for  v  <  3 
cm'1  are  shown  in  Figure  2  so  that  the  details 
at  low  frequencies  are  apparent.  Similar  graphs  for 


frequencies  as  high  as  2060  cm  1  are  given  by  Hill 
et  al.  [1980a,  6], 

The  graphs  of  absorption  coefficient  show  the 
water  vapor  resonances  at  v  <  16  cm  “ 1  and  a  broad 
underlying  absorption  from  the  wings  of  the  very 
strong  lines  at  v  >  16  cm'1.  The  absorption  coeffi¬ 
cient  is  given  in  nepers  per  centimeter.  To  convert 
to  decibels  per  kilometer,  multiply  our  values  by 
4.34  x  10s.  Of  course,  the  microwave  resonances 
of  02  and  the  excess  absorption  are  not  shown. 
There  are  43  water  vapor  resonances  at  v  <  16  cm  ~ 1 
and 5  at  v  <  3  cm'1.  We  have  numerically  resolved 
all  these  resonances,  although  only  the  seven  stron¬ 
gest  resonances  are  evident  in  the  graphs. 

The  graphs  of  refractivity  Na  show  the  usual 
effect  of  a  resonance  superposed  on  a  decreasing 
background.  The  decreasing  behavior  of  this  back¬ 
ground  is  dominated  by  the  strong  resonance  at 
SS6  GHz.  A  constant,  i.e.,  nondispersive,  behavior 
of  this  background  is  obvious  at  the  lowest  frequen¬ 
cies.  This  nondispersive  value  of  Na  for  v  — »  0  is 
caused  by  the  sum  of  all  the  strong  resonances 
at  v  >  16  cm"1.  As  v  — *  0,  each  line  contributes 
a  refraction: 

Nal—*  (S,C/2h2i\)  X  106 

The  sum  of  all  these  contributions  adds  to  the 
relatively  Q  insensitive  refraction  at  optical  fre¬ 
quencies  to  provide  the  relatively  strong  Q  depen¬ 
dence  at  radio  frequencies.  We  see  that  the  strongest 
contributions  come  from  the  rotational  band  where 
S,  is  largest  and  v,  is  smallest. 

We  graph  the  contribution  of  the  resonances  to 
Ar,  namely,  Ar  +  Nd  x  10~6.  To  obtain  Ar,  one 
must  add  the  continuum  contribution  -Ndx  10  6 
=  -2.62  x  10'“  to  the  values  graphed. 

The  pressure-broadened  02  band  near  60  GHz 
has  a  peak  absorption  about  1/20  as  large  as  the 
prominent  water  vapor  line  at  183  GHz  (6.1  cm'1) 
for  our  ( T ),  (/*),  and  {Q).  The  results  of  Liebe 
and  Hopponen  [1977]  show  that  the  contribution 
of  this  02  band  to  refractivity  and  /tris  also  smaller 
than  that  shown  in  Figure  1  for  the  183-GHz  water 
vapor  line.  The  02  band  near  119  GHz  is  much 
weaker  than  that  at  60  GHz  and  has  an  even  smaller 
effect  on  refractivity  and  Ar. 

The  graphs  of  A  0  show  that  the  infrared  reso¬ 
nances  dominate  the  A  Q  value,  since  the  continuum 
contribution  is  only -1.2  x  10*‘ for  our  (T),  (/*), 
(Q).  Except  for  the  change  in  scale,  the  AQ  graphs 
look  like  those  of  Na .  This  is  not  surprising,  because 


Fig.  1.  Absorption  coefficient,  the  refractivity  and  the  sensitivity  coefficients  for  v  <  16  cm  '  and  mean 
conditions:  T  «  300  K,  P  —  1  atm,  80%  relative  humidity,  i.e.,  partial  pressure  of  water  vapor  is  28.2  mbar. 


neglecting  the  continuum  and  line  width  contribu¬ 
tions  to  Ag  gives  Ag  proportional  to  Nm. 

Th  graphs  for  Ar  are  not  shown  here.  They 
are  given  by  Hill  el  al.  ( 1980b  ] .  The  line  contribution 
to  A  r  is  localized  to  each  line  so  that,  unlike  Na , 
it  has  no  ‘background’  value  from  strong  resonances 
at  v  >  16  cm'1.  This  line  contribution  to  Ar  oscil¬ 
lates  about  zero,  taking  very  small  values  between 
resonances,  with  peak  values  of  only  ±6  x  10~* 
for  v  <  16  cm"1.  This  is  to  be  compared  with  the 


Ar  value  from  the  continuum  of  Nd  x  10~6  =  2.62 
x  10~4  for  our  ( T ),  (P),  (Q).  Thus  ^ris  essen¬ 
tially  the  same  for  v  <  16  cm'1  as  it  is  for  visible 
frequencies. 

A  test  of  our  computations  is  to  compare  our 
values  of  refraction,  Ar,  A,,  and  Ag  with  known 
values  at  radio  wavelengths.  We  compare  our 
computed  values  at  v  =  10 ~J  cm'1,  i.e.,  K  =  1  m, 
with  the  radio  refractive  index  given  by  Thayer 


[1974] .  Our  Hd  and  Nw  both  differ  by  0.16%  from 
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Fig.  2.  Details  of  the  refractivity  N.  and  A  T  and  A„  for  v  <  3 
cm'1. 

those  derived  from  Thayer’s  equation.  Our  Na  value 
is  2.5%  smaller  than  that  derived  from  Thayer’s 
equation.  Zhevakinand  Naumov  [  1967*]  calculated 
Na  by  adding  contributions  of  strong  lines;  their 
value  differed  from  that  at  centimeter  wavelengths 
by  2.7%.  The  contributions  of  the  resonances  to 
Ap  is  -7  x  10'12,  compared  with  the  dry  air  con¬ 
tribution  2.61  x  10~4;  so  the  0.16%  difference  of 
our  A,  from  that  obtained  from  Thayer’s  equation 
is  caused  by  the  different  Nd  values.  The  contribu¬ 
tion  of  water  vapor  resonances  to  A  T  is  5%  smaller 
in  absolute  value  than  \A  T  +  Nd  x  10  6 1  calculated 
from  Thayer’s  equation.  However,  we  have  not 
added  contributions  of  infrared  resonances  of  02 


or  C02  or  other  species  to  our  evaluation  of  Ar 
or  the  value  of  Nd  at  radio  frequencies.  The  con¬ 
tribution  of  all  resonances  to  AQ  is  2.4%  smaller 
than  that  obtained  from  Thayer’s  equation.  We  do 
not  know  whether  or  not  the  agreement  between 
our  calculated  Na  and  the  known  value  at  radio 
wavelengths  implies  that  the  strengths  and  reso¬ 
nance  frequencies  of  strong  rotational  lines  in  the 
AFGL  compilation  are  correct  or  if  the  agreement 
is  fortuitous  due  to  cancelling  errors. 

The  graphs  of  BQ  appear  to  be  similar  to  those 
of  the  absorption  coefficient.  This  is  because  of 
the  approximation  BQ  «  3/4-irv,  obtained  by  ne¬ 
glecting  the  Q  dependence  of  the  line  width.  We 
see  that  humidity  fluctuations  cause  the  strongest 
absorption  fluctuations  where  absorption  is  stron¬ 
gest.  Kanevskii  [1972]  considered  the  effect  of 
absorption  fluctuations  on  the  amplitude  of  a  plane 
wave  for  frequencies  10  cm'1  <  v  <  14  cm"1.  He 
shows  large  effects  on  the  amplitude  variance  at 
frequencies  where  resonances  produce  large  ab¬ 
sorption,  and  practical  propagation  paths  would  be 
very  short.  Otherwise  the  effect  is  small. 

We  show  a  graph  of  BT/BQ  to  illustrate  the 
relative  importance  of  temperature  and  humidity 
fluctuations  to  absorption  fluctuations.  Between 
lines,  Bt/  Bq  has  a  value  of  -2.1  with  excursions 
to  smaller  absolute  values  at  the  resonances.  The 
results  of  Liebe  and  Hopponen  [1977]  show  that 
the  pressure-broadened  02  band  dominates  BT  near 
60  GHz.  Of  course,  Br  includes  the  effect  of  changes 
in  02  concentration,  because  the  ideal  gas  law  shows 
that  for  fixed  total  pressure,  water  vapor  concentra¬ 
tion,  and  fixed  ratio  of  02  concentration  to  those 
of  N2  and  Ar,  any  temperature  change  causes  an 
opposite  change  in  02  concentration. 

We  do  not  show  graphs  of  Br  here.  Graphs  of 
Bp/  Bq  are  given  by  Hill  ei  al.  [  19806  ] .  These  graphs 
show  Brl Bq  close  to  0.8  between  resonances  with 
excursions  to  nearly  -1.0  at  the  stronger  reso¬ 
nances. 


DISCUSSION 

In  atmospheric  turbulence,  a  typical  order  of 
fractional  fluctuations  is 


5 P  bT  bQ 

- c - s - 

< p )  <r>  (Q) 

In  view  of  the  relative  values  of  Af  and  Ar  and 
those  of  Bp  and  BQ,  we  neglect  the  contributions 
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of  total  pressure  fluctuations  to  both  the  real  and 
imaginary  parts  of  refractive  index.  In  the  case 
of  very  strong  humidity  fluctuations,  (8 T/(T))  <c 
(&Q/(Q )),  one  may  neglect  the  temperature  fluc¬ 
tuations  relative  to  the  humidity  fluctuations  for 
fluctuations  in  both  real  and  imaginary  parts  of 
the  refractive  index,  because  Ar/  A Q  -  -3  for  v  < 
16  cm-1,  and  the  absolute  value  of  BTj BQ  does 
not  exceed  2.2.  One  must  be  careful  for  the  case 
of  absorption  fluctuations,  because  inclusion  of 
excess  absorption  may  significantly  increase  the 
magnitude  of  Br/  BQ  in  the  windows  between  ab¬ 
sorption  lines.  At  submillimeter  and  lower  frequen¬ 
cies,  say,  v  <  50  cm"1,  no  contribution  by  excess 
absorption  is  added  to  our  calculated  values  of  (3, 
Bg,  Bt/Bq,  or  BPI  Bq  for  reasons  discussed  pre¬ 
viously.  Readers  may  use  their  favorite  formula 
for  excess  absorption  as  a  function  of  v,  T,  Q, 
and  P  to  add  these  contributions  to  the  values  in 
our  graphs. 

Neglecting  pressure  fluctuations.  Hill  et  al. 
[1980a]  give  the  relationship  between  structure 
parameters  of  both  real  and  imaginary  parts  of  the 
refractive  index  and  structure  parameters  of  tem¬ 
perature  and  humidity.  They  also  consider  the 
variance  of  log-amplitude  to  show  that  refraction 
fluctuations  give  much  larger  amplitude  fluctuations 
than  do  absorption  fluctuations  for  most  cases  of 
practical  interest. 

Acknowledgments.  We  (hank  Dr.  Hans  Liebc  for  his  helpful 
discussions. 
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A  laser  rain  gauge 
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Detection  of  the  scintillations  caused  by  raindrops 
falling  through  a  laserbeam  can  serve  as  the  basis 
of  a  rain  gauge.  The  laser  technique  measures  rain¬ 
fall  averaged  over  paths  several  hundred  meters 
long.  Because  the  measurements  can  be  made  on  line 
without  extensive  data  reduction,  the  technique  ap¬ 
pears  promising  for  practical  use. 

Precise  measurements  of  rainfall  are  nowhere 
near  as  simple  as  you  might  imagine.  The  problem 
is  that  rainfall  is  highly  variable  from  point  to  point, 
and  conventional  rain  gauges  simply  measure  the 
amount  of  rain  that  falls  at  a  single  point.  However, 
for  agricultural  and  flood-control  purposes,  it  is 
total  rainfall  over  an  area  that  is  important.  To 
measure  that  quantity,  a  dense  network  of  conven¬ 
tional  gauges  is  needed.  We  were  able  to  perform 
such  a  measurement  by  observing  the  fluctuations 
in  intensity  of  a  laser  beam  passing  through  several 
hundred  meters  of  rain-filled  atmosphere  near  the 
ground. 

We  have  shown  that  the  scintillations  produced 
by  raindrops  falling  through  a  collimated  laser  beam 
can  be  used  to  measure  the  drop-size  distribution 
and  the  rainfall  rate.' 2  both  averaged  over  the  path. 
We  have  also  shown  that  the  variance  of  the  scin¬ 
tillations  measured  by  a  linear  detector  at  fre¬ 
quencies  near  one  kilohertz  is  closely  related  to  rain 
rate2  and  is  nearly  independent  of  drop-size  distri¬ 
bution.  If  only  rain  rate  is  desired,  the  variance- 
type  optical  rain  gauge  is  simple  to  construct.  It 
can  use  a  diverging  beam,  thus  eliminating  the 
practical  difficulties  in  maintaining  adjustment 
and  pointing  of  a  collimated  beam.  Furthermore, 
it  is  less  sensitive  to  the  presence  of  updrafts  and 
downdrafts  along  the  beam  and  can  thus  be  used 
over  rough  terrain.  Here  we  describe  the  basic  prin¬ 
ciple  and  the  design  of  an  optical  rain  gauge  using 
the  variance  of  the  irradiance  scintillations  induced 
by  raindrops. 

For  a  plane  wave  incident  on  a  raindrop,  the 
scattered  light  can  be  approximated  by  two  com¬ 
ponents;  the  light'diffracted  by  the  sharp  boundary 
of  the  drop  and  the  light  refracted  by  passage 
through  the  drop.4  The  refracted  light  is  much 
weaker  than  the  diffracted  light  because  the  wave¬ 
length  of  the  incident  wave  is  much  shorter  than 
the  radius  of  the  drop,  and  the  refractive  index  is 
large  (i.e.,  fcaAn»l,  where  k  is  the  wavenumber,  o 
is  the  radius  of  the  drop,  and  1  +  An  is  the  re'*active 

Robert  S  Lawrence entf  Ting-i  Wang  are  both  special's ts  m  cct'cat 
remote  sensing  at  the  National  Oceanic  ana  Atmospheric  Adminisva 
tion's  Wave  Propagator  Laboratory  in  Boulder  Lawrence  received  a 
master's  degree  m  physics  in  I960  from  Yale  University  Wang  re 
ceived  a  doctorate  m  electrical  engineering  m  1973  from  Dartmouth 


Fig  1  The  interference  pattern  produced  by  a  spherical  rmndrop 
The  fine  fringes  are  caused  by  phase  interference  of  the  incident 
wave  and  scattered  wave 


index  of  water).  Thus  the  total  scattered  light  is 
approximately  equal  to  the  diffracted  light. 

The  scattered  field  can  be  considered  as  a  spheri¬ 
cal  wave  emitted  from  the  center  of  the  raindrop 
because  we  are  only  interested  in  the  far-field  solu¬ 
tion  (i.e.,  where  the  drop  radius  is  much  less  than  the 
Fresnel-zone-size.a  «Ux)'",  where  z  is  the  distance 
between  the  drop  and  the  receiver,  and  the  trans¬ 
verse  coordinate  is  much  less  than  the  longitudinal 
coordinate).  The  interference  at  the  receiving  plane, 
caused  by  the  incident  plane  wave  and  the  scattered 
spherical  wave,  gives  the  fine  fringes  shown  in  Fig. 
1.  The  envelope  of  the  interference  pattern  is  the 
well-known  Airy  diffraction  pattern  of  a  sphere. 
Path-averaged  rain  rates  can  be  obtained  by  observ¬ 
ing  the  movement  of  these  raindrop-induced  inter 
ference  patterns. 

We  have  shown2  that,  for  a  line  detector  extend 
ing  in  the  horizontal  direction,  the  variance  of  the 
irradiance  scintillations  is 

o2~A  (1  -  0.045IA0  21  (/, 2  +  /,/,  +  /,’IJ  (1) 

where  h  Is  the  path-averaged  rain  rate  (milli 
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meters  per  hour),  and  /,  and  f,  are  the  low  and  high 
cutoff  frequencies  (kHz)  respectively.  Obviously 
when  /,  and  f,  are  small,  the  variance  is  proportional 
to  the  path-averaged  rain  rate.  However,  if  the 
observed  frequency  is  too  low,  contamination  by 
turbulence-induced  scintillations  cannot  always 
be  ignored,  particularly  when  the  wind  is  strong.  If 
the  observed  frequency  is  too  high,  the  measure¬ 
ment  is  sensitive  to  drop-size  distribution  which 
usually  varies  with  location  and  time.  We  found 
that  a  narrowband  filter  around  1  kHz  is  an  optimum 
choice.  In  this  case  if,  »/j  =  1  kHz),  Eq.  1  becomes 

[1  -  0.135A0  J,J  (2) 

For  rain  rates  of  0.1  mm/hour  to  100  mm/hour,  we 


can  approximate  Eq.  2  by 


Eq.  3  is  derived  based  on  the  assumption  of  a  plane 
wave  such  as  from  a  collimated  beam.  It  is  not  dif¬ 
ficult  to  show  that  for  a  spherical  wave  (point 
source),  the  relationship  between  rain  rate  and  the 
variance  of  the  irradiance  scintillations  is 


For  a  diverging  beam,  the  dependence  will  be  be¬ 
tween  ol  01  and  oJ '. 

A  schematic  diagram  of  the  system  is  shown  in 
Fig.  2.  The  transmitter  is  a  helium-neon  laser.  In 
front  of  the  detector  is  a  thin  slot,  20  centimeters 
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Fig  3  A  compamon  shows  excellent  agreement  between  the  opti¬ 
cally  measured  ram  rate  I upper  diagram!  and  the  average  ram  rate 
measured  by  two  tipping- bucket  rain  gauges  near  the  ends  of  the 
optical  path ! lower  diagram  l 

long  and  0.1  cm  high,  extending  horizontally.  A 
narrowband  interference  filter  excludes  back¬ 
ground  light-  Two  lenses  focus  the  light  on  two  p,-». 
silicon  detectors.  Logarithmic  amplifiers  normalize 
the  fluctuating  signals  with  respect  to  unperturbed 
light  intensity  so  that  a  change  in  transmitter  in¬ 
tensity  or  a  movement  of  the  beam  will  not  affect 
the  measurements.  The  subtraction  of  the  signals 
from  the  log  amplifiers  cancels  the  effect  of  noise 
common  to  both  detectors  (such  as  120-hertz  ripple 
due  to  imperfect  filtering  of  the  power  supply),  and 
makes  the  instrument  less  sensitive  to  vibrations 
and  changes  in  the  transmitted  intensity.  The  out¬ 
put  is  fed  into  a  narrow-bandpass  electronic  filter 
with  center  frequency  1  kHz  and  rejection  ratio 
greater  than  60  decibels  between  1  kHz  and  200  Hz. 
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Fig  4  A  linear-regression  scatter  diagram  of  the  rain  rates  measured 
by  the  laser  system  and  ram  rates  measured  by  the  average  of  two 
conventional  tipping- bucket  ram  gauges 


The  standard  deviation  of  the  filtered  signals  is 
calculated  by  an  Analog  Device  root-mean-square 
integrated  circuit  (AD536)  with  a  10-second  time 
constant.  A  log  and  antilog  (exponential)  integrated 
circuit  (AD433)  calculates  the  rain  rates.  To  prevent 
false  measurements  due  to  low  signals,  the  outputs 
of  the  photodetectors  are  amplified  and  clipped. 
When  either  signal  is  below  the  threshold,  the  com¬ 
parator  flips  the  switch  so  that  a  voltage  indicating 
"signal  off’  is  shown  on  the  display. 

A  comparison  of  the  rainfall  rates  measured  by 
the  optical  technique  and  conventional  tipping- 
bucket  rain  gauges  is  shown  in  Fig.  3.  Because  of 
the  poor  time  resolution  of  the  tipping-bucket  rain 
gauges,  we  use  a  one-minute  time  constant  on  both 
records.  Two  tipping-bucket  rain  gauges  were 
placed  near  the  ends  of  the  200-m  optical  path.  The 
average  of  the  two  rain  gauges  is  shown  in  the  lower 
diagram.  The  agreement  between  the  optical  mea 
surements  and  the  conventional  rain  gauge  mea¬ 
surements  is  excellent,  with  a  correlation  coefficient 
of  0.97.  The  linear  regression  of  these  two  quantities 
is  shown  in  Fig.  4.  For  rain  rates  above  about  10 
millimeters  per  hour,  the  results  are  plotted  for  a 
fixed  one  minute  time  constant  (shown  by  dots  in 
Fig.  4).  For  low  rain  rates,  the  time  resolution  of 
the  tipping  bucket  rain  gauge  is  insufficient  to 
compare  with  the  optical  measurements.  We  de 
cided  to  take  the  average  of  the  optical  data  be¬ 
tween  tips  of  the  tipping-bucket  gauge  and  to  plot 
this  against  the  tipping-bucket  gauge  reading 
(shown  by  big  dots  in  Fig.  4).  The  agreement  is  very 
good  and  might  have  been  further  improved  by- 
using  more  rain  gauges  along  the  path. 
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Data  Analysis  of  the 

NOAA-GIT  Millimeter-Wave  Propagation  Experiment 

Near  Flatville,  Illinois 


R.  J.  Hill,  W.  P.  Schoenfeld1,  J.  P.  Riley,  J.  T.  Priestley, 

S.  F.  Clifford,  S.  P.  Eckes1,  R.  A.  Bohlander2,  R.  W.  McMillan2 


ABSTRACT.  The  authors  discuss  the  quantities  calculated  from  the 
millimeter-wave  and  microraetorological  data  base.  The  data  base  was 
obtained  from  the  millimeter-wave  propagation  experiment  conducted 
near  Flatville,  Illinois,  by  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  and  Georgia  Institute  of  Technology  (GIT). 

The  repair  of  data-base  errors  by  digital  computation  is  also  dis¬ 
cussed,  and  the  overall  structure  of  the  data-processing  programming 
is  described. 


1.  INTRODUCTION 

The  millimeter-wave  propagation  experiment  consisted  of  five  data-taklng 
sessions.  The  dates  for  these  sessions,  including  on-site  setup  and  takedown, 
were  June  3  to  July  15,  1983;  October  30  to  December  11,  1983;  February  5  to 
March  9,  1984;  May  29  to  July  5,  1984;  and  January  29  to  March  26,  1985. 

Figure  1  is  a  self-explanatory  diagram  of  the  experiment.  Georgia  Institute 
of  Technology  (GIT)  was  responsible  for  the  mil lime ter- wave  data,  and  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA)  for  the  mlcrometeorolo- 
glcal  data.  The  methods  of  analysis  for  the  clear-air  portion  of  the  microme- 
teorological  data  are  described  here;  analyses  of  the  meteorological  data  for 
fog,  rain,  and  snow  are  not  discussed.  On  the  other  hand,  the  methods  of 
analysis  for  the  millimeter-wave  data  given  here  are  identical  for  clear  air, 
fog,  rain,  or  snow. 

^Cooperative  Institute  for  Research  in  Environmental  Sciences,  University 
of  Colorado  and  National  Oceanic  and  Atmospheric  Administration,  Boulder,  CO 
80309. 

2 

Engineering  Experiment  Station,  Georgia  Institute  of  Technology,  Atlanta, 

GA  30332. 


The  philosophy  underlying  the  data  processing  is  as  follows.  We  calculate 
all  quantities  of  interest  using  as  few  computer  programs  as  possible.  For 
instance,  we  avoid  programming  that  has  the  special  purpose  of  calculating 
intensity  covariance  or  the  mutual  coherence  function.  All  such  diverse  sta¬ 
tistics  are  calculated  through  one  program,  thereby  reducing  the  complexity  of 
the  program  library  and  reducing  input /output  costs.  We  use  visual  inspection 
of  the  data  from  computer-generated  microfilm  stripcharts  that  reveal  full 
data  time  series.  Quantities  are  also  averaged  into  time  intervals  of  2.56 
and  25.6  s  and  graphed  to  reveal  trends  and  other  problems.  We  examine  graphs 
of  power  spectra  and  cross-spectra  to  detect  noise  in  the  data.  Only  by  such 
means  of  data  visualization  can  we  be  sure  that  the  myriad  errors  within  the 
data  can  be  identified  and  avoided.  This  visualization  of  the  data  is  thus 
essential  to  our  data  validation.  Finally,  the  end  product  of  the  data  analy¬ 
sis  programs  is  what  we  call  answer  files  that  contain  the  statistics  from 
statistically  stationary  periods  on  each  data  tape.  Thus  investigators  are 
not  limited  in  their  interpretation  of  the  data  by  any  a_  priori  decisions  as 
to  how  the  results  are  to  be  presented.  For  instance,  the  answer  files  may  be 
Interrogated  to  produce  plots  of  the  mutual  coherence  function  for  a  variety 
of  surface-layer  stabilities,  each  stability  case  arising  from  different  data 
runs.  There  remains  great  flexibility  in  graphing  and  tabulating  the  answer 
files  for  the  purpose  of  making  scientific  discoveries.  Investigators  are 
then  limited  only  by  their  insight  and  ingenuity. 

2.  QUANTITIES  COMPUTED  FROM  THE  GIT 
MILLIMETER-WAVE  DATA 

2.1  Introduction 

The  GIT  millimeter-wave  data  consist  of  28  channels  taken  at  100  Hz.  (In 
this  report  "millimeter-wave"  is  abbreviated  as  mm-wave.)  Figure  1  shows  the 
arrangement  of  receiver  antennas,  housed  in  a  trailer,  as  viewed  from  the 
transmitter;  the  antennas  are  designated  by  a  number.  There  are  six  possible 
antenna  pairs,  each  with  a  different  spacing.  The  antenna  pair  numbers  and 
their  spacings  are  given  in  Table  1.  For  a  few  experiment  runs,  the  fifth 
antenna,  below  antenna  2  (see  Fig.  1),  was  used  in  place  of  antenna  4;  the 
objective  was  observation  of  anisotropy  or  multipath  effects. 


Figure  1. — The  instrument  positions  at  the  experiment  site.  The  dashed 
and  dotted  line  denotes  ram-wave  propagation  path  (1.4  km);  the  long- 
dashed  lines,  the  optical  propagation  paths  (1  km  each);  and  the  short- 
dashed  lines,  the  optical  rain  gauge  paths  (50  m  each).  Solid  lines 
show  the  flow  of  micrometeorologlcal  data  to  the  data  aqulsltion  system 
in  the  receiver  trailer.  Antennas  are  numbered  1  to  5  in  the  receiver 
trailer. 


Table  1. — Antenna  pairs 
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1  &  2 

1.43 
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1  &  3 

4.29 
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1  &  4 

10.00 
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2  &  3 

2.86 
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2  4  4 

8.57  | 
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3  4  4 

5.71  ! 
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For  each  antenna  pair,  the  GIT  ram-wave  data  contain  the  square  of  the 
voltage  from  four  phase  shifts  of  0°,  90°,  180°,  and  270°.  For  the  antenna 
pair  formed  by  antennas  i  and  j  these  are  denoted  by 

<Vn*ij  *  ^V90^ij  *  (V180^ij  *  and  ^V270*ij  *  ^ 

The  data  also  contain  the  intensity  at  each  antenna,  denoted  by  1^  for  the  ith 
antenna. 

It  is  important  to  keep  in  mind  that  the  intensities  are  measured  in 
arbitrary  units  and  that  differences  exist  between  the  time-averaged 
(hereafter  referred  to  as  mean)  intensities  of  different  antennas.  Thus  any 
statistics  calculated  from  the  intensities  must  be  normalized  by  the 
appropriate  mean  intensities.  Because  of  hardware  characteristics,  the  phase 
difference  between  antenna  pairs  has  an  arbitrary  offset  that  is  different 
for  each  pair.  Thus  any  phase-difference  statistics  must  be  calculated  from 
the  phase  difference  with  its  mean  value  subtracted. 

Drift  in  the  mm-wave  receivers  causes  the  measured  intensities  and  squared 
voltages  in  (1)  to  change  slowly  with  time.  Thus  it  is  necessary  to  low-pass 
filter  the  mm-wave  data;  the  filter  frequency  can  be  varied  as  needed.  We 
redefine  the  phase  difference  as  its  value  with  its  low-pass-filtered  value 
subtracted.  To  excellent  approximation,  the  mean  phase  difference  is  then 
zero.  We  redefine  the  intensity  as  its  value  divided  by  its  low-pass-filtered 
value.  Then,  to  excellent  approximation,  the  mean  intensity  is  unity.  In  the 
following  we  discuss  various  statistical  quantities  for  which  mean  phase  dif¬ 
ferences  are  subtracted  and  mean  intensities  are  divided  out.  One  must  keep 
in  mind  that  if  the  low-pass  filtering  is  used,  then  a  subsequent  subtraction 
of  the  mean  phase  difference  and  scaling  by  the  mean  intensity  has  essentially 
no  effect.  However,  if  the  low-pass  filtering  is  not  used,  then  the  computer 
code  continues  to  operate  correctly  by  subtracting  the  nontrivial  mean  phase 
difference,  and  scaling  by  the  nonunity  mean  intensities. 


2.2  Phase-Difference  Statistics 

The  phase  difference  for  antenna  pair  i  and  j  is  computed  using  the 
following: 
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The  Fortran  ATAN2  function  is  used  in  the  computations.  This  gives  0^  in  the 
range  -it  to  it.  That  range  is  extended  to  — 2 tt  to  2tt  by  following  the  tran¬ 
sition  in  the  value  of  ATAN2  as  the  phase  difference  crosses  -ir  or  tt.  It  is 
possible  to  extend  the  range  to  -•  to  «  by  slightly  more  complicated  program¬ 
ming;  this  has  not  been  done  since  it  is  expected  that  the  range  from  -2n  to 
2k  will  never  be  exceeded.  In  any  case,  we  will  be  aware  from  microfilm 
stripcharts  if  this  range  is  exceeded. 


For  each  antenna  pair  we  calculate  mean  phase  difference  <0^>  and  the 
mean  of  the  squared  phase  difference  <(6jj)2>  from  which  we  obtain  the  phase- 
difference  variance 


<(6ij  “  <Pij»2>  *  <(Pij)2>  “  <0ij>2  *  (3) 

These  six  variances  produce  the  phase  structure  function,  when  graphed  as  a 
function  of  antenna  spacing. 


2.3  Intensity  Statistics 

For  each  antenna,  we  calculate  the  mean  intensity  <1^  and  the  mean  of  the 
squared  intensity  <I12>.  From  these  we  obtain  the  normalized  intensity 
variance 
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For  each  of  the  antenna  pairs  i  and  j ,  we  calculate  the  cross-product  average 
<1^  Ij>;  this  gives  the  normalized  covariance 
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When  graphed  versus  antenna  spacing,  the  six  are  the  normalized  spatial 
covariance.  On  the  basis  of  inertial-range  theory  and  no  absorption  fluc¬ 
tuations,  the  minimum  antenna  separation  of  1.43  m  implies  that  the  normalized 
covariance  would  be  only  about  0.1,  and  be  very  much  smaller  at  larger  separa¬ 
tions.  Thus  the  effects  of  absorption  fluctuations  might  be  seen  in  the 
large-spacing  covariance  values  because  the  absorption  fluctuations  affect 
much  larger  spatial  scales  in  the  scintillation  pattern  than  do  refractive 
fluctuations. 


2.4  Mixed  Intensity  and  Phase-Difference  Statistics 

For  each  antenna  pair,  we  calculate  the  intensity-phase-difference  cross¬ 
correlations  <I1  and  <1^  8^>.  These  are  used  to  obtain  the  quantity 


1  <pL _ iL.)  b  >  „  I  /<Ii.Al>  _  i1!  V\ 

2  ^<it>  <ij> )  2  \  <it>  <ij>  ) 


define  the  log-amplitude  by  x^  =  An(A1/<Ai>).  Since 
<I>  =  <A2>  -  <(A  -  <A>)2>  +  <A>2  , 


and  for  our  case  of  weak  scintillation  the  amplitude  variance  is  much  smaller 

2 

than  the  square  of  the  mean  amplitude,  then  <I>  *  <A>  •  So  the  log-amplitude 
can  be  approximated  by 
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Therefore,  to  first  order  in  the  smallness  of  intensity  fluctuations  relative 
to  their  mean,  which  is  an  excellent  approximation,  we  have  the  quantity 


<(xi  -  Xj)  ®ij>  =  2  <(<i1>  '  oj)  eij: 


(12) 


2.5  Second  Moment  of  the  Field:  The  Mutual 
Coherence  Function 

Let  TI^  denote  the  electromagnetic  field  at  the  ith  antenna,  and  let  U^* 
denote  the  field's  complex  conjugate  at  antenna  j.  The  second  moment  of  the 
field,  known  as  the  mutual  coherence  function,  MCF,  is  <U^  U  *>.  We  do  not 

JL  *' 

calculate  the  case  i  *  j  because  <U^  >  ■  <I^>  is  the  average  intensity, 

which  is  already  calculated  as  in  Sec.  2.3.  If  we  write  the  field  as 


Aie 


(13) 


where  A^  -  /T^  is  its  amplitude  and  ^  is  its  phase,  then  the  MCF  is 


*  i($i  "  ♦<) 

<Ut  Uj  >  -  <At  e  J  > 

i(Bi  -  <B» 

-  <Kt  AjC  1J  >  . 

Here  the  phase  difference  is 
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since  from  Eq.  (2)  is  the  phase  difference  to  within  an  arbitrary  offset, 
this  offset  being  removed  by  subtracting  The  unnormalized  MCF  is  not 

meaningful  because  the  intensities  are  known  only  in  arbitrary  units.  The 
normalized  MCF  is  denoted  by  M2  and  is  given  by 
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We  use  the  relationship 
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e  J  «  cos(8.  .)  +  i  sin( 8. .)  • 
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We 


then  average  the  real  and  imaginary  parts  of  the  numerator  of  M^;  namely. 


we  calculate  the  averages  of  and  J  where 


Rij  =  Xj  cos(61j) 


(20) 


Jij  E  /xi  sin<ei;J)  • 


(21) 


Note  that  we  must  take  the  sine  and  cosine  of  8^ j »  which  itself  is  an  arc¬ 
tangent,  rather  than  using  identities  for  trignometric  functions  of  the  arc¬ 
tangent  because  we  must  use  a  temporal  filter  on  8jj  (see  Sec,  4.3)  but  we 
must  not  filter  the  quantities  in  (1).  The  inverse  of  the  denominator  of 
M2  can  be  written  as  RQ  -  iJD,  where 


Rd  -  cos(<Sij>)//<I1><Ij> 

(22) 

JD  -  8in(<8lj>)//<Ii><Ij>  , 

(23) 

so 

^  -  ( <R1 j >  +  i<Jij>)( Rq  -  iJD)  . 

(24) 

The  normalization  of  the  MCF  is  performed  in  real  arithmetic;  the  real  and 
imaginary  parts  of  are  determined  from 


Re(M2)  -  <Rij>RD 

+  <Jtj>JD 

(25) 

Im(M2)  =  <Jij>R0 

<Rij>JD  * 

(26) 

The  modulus  and  phase  of  M2  are  also  calculated  using 


r  -  {[Re(M-)]2  +  (Im(M  )]2}l/2 


♦  =  arctanf Im(M2)/Re(M2) J 


such  that 


*  r  e 


i* 


(27) 

(28) 

(29) 


2.6  Fourth  Moment  of  the  Field 

Having  four  receiving  antennas  allows  us  to  calculate  the  fourth  moment  of 
the  field,  which  is  given  by 


<U.  U.*  0  U  *>  . 
1  j  n  m 
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Ifi“j»n*m  then  ■  <1^  >,  which  is  the  mean  of  the  squared  intens.ity 
and  has  been  calculated  in  Sec.  2.3.  If,  for  instance,  i  -  j  and  n  ■  m,  then 
the  fourth  moment  becomes  <1.  I  which  has  also  been  calculated  in  Sec.  2.3. 

Thus,  not  all  possible  combinations  of  i,  j,  n,  and  m  are  of  interest.  In 
fact,  certain  choices  of  i,  j,  n,  and  m  merely  produce  the  complex  conjugate 
fourth  moment  arising  from  a  different  choice  of  indexes  or  produces  the  same 
fourth  moment  because  of  multiplicative  commutation.  Thus  we  calculate  only 
those  fourth  moments  that  are  distinct  from  intensity  products  and  are 
distinct  under  complex  conjugation  and  multiplicative  commutation.  The  sets 
of  i,  j,  n,  m  that  produce  these  distinct  values  of  are  given  in  Table  2. 

t 

Similar  to  the  case  of  the  second  moment  of  the  field,  only  the  normalized 
fourth  moment  is  meaningful  because  the  intensities  are  known  only  in 


Table  2. — Indexes  for  producing  distinct  values  of  M. 


arbitrary  units  and  the  phase  differences  are  known  only  to  within  additive 
constants.  The  normalized  fourth  moment  Is  defined  by 


<U.  U.*  U  U  *> 
1  J  n  m 


/<U.  U.*XU.  U.*XU  U  *><U  U  *> 
11  jj  nn  mm 


By  analogy  with  the  second  moment,  this  gives 


<A  rj  "PUCSu  *  6„„»> 


/<I  XI  XI  XI  >  exp[l(<8  >  +  <6  »] 

1  j  n  m  lj  nm 


In  the  computation  we  average  the  real  and  Imaginary  parts  of  the  numerator  of 


;  namely,  we  obtain  the  averages  <R^>  and  <J^>  where 


<R,>  =  </I.  I,  I  I  cos(8.  +  0  )> 

4  1  j  n  m  lj  nm 


<V  =  </ri  li  rn  rm  sln(8lj  +  0nm» 


For  the  first  case  In  Table  2,  R^  and  are  calculated  by  noting  that 


R4  "  R12  R34  "  J12  J34 


J4  *  J12  R34  +  R12  J34 


where  R^  and  were  previously  calculated  values  of  the  real  and  Imaginary 


parts  of  the  Instantaneous  second  moment  and  are  given  In  Eqs.  (20)  and  (21). 
The  second  and  third  cases  In  Table  2  are  calculated  similarly.  Cases  4  to  27 
in  Table  2  have  one  pair  of  equal  indexes ,  1  *  j .  These  are  calculated  by 
noting  that 


r4  *  h  \m 


J4  ’  !1  J™ 


where  1^  Is  the  Intensity  at  antenna  1  (and  also  j  since  i  *  j);  Rnm  and 


J  are  the  real  and  imaginary  parts  of  the  instantaneous  second  moment  for 
nm 


the  remaining  antenna  pair  as  given  in  Eqs.  (20)  and  (21).  Finally,  cases 
28  to  33  in  Table  2  are  the  squares  of  instantaneous  second  moments  (since 


a  -  i  and  m  ■  j)  and  are  thus  calculated  using 


\  -  <V2 '  (V2 


J4  "  2  Rij  Jij  * 


The  inverse  of  the  denominator  of  can  be  written  as  -  iJ4D,  w*'ere 


R4D  5  +  <S„„»/d 

J4D  5  +  <*„„»'« 


<*  *  Ml>  <IJ>  <V  <V  • 


«4  -  «V  ♦  1<V><R4D  -  lJ4D>  * 


For  cases  1  to  3  in  Table  2  these  definitions  of  R^D  and  J^D  are  used 
literally.  For  cases  4  to  27,  in  which  i  *»  j ,  we  have  Bjj  ■  0  and  1^  ■  I  j ,  so 
the  following  simplified  expressions  are  used: 


R4D  -  cos«Bnm»/d 


J4D  *  8ln(<6nm>)/d 


with  d  -  <Ij>  /<In>  <Im>.  These  trigonometric  evaluations  and  the  square  root 
in  d  are  the  same  as  were  used  in  normalizing  the  second  moment,  so  they  need 
not  be  computed  a  second  time.  For  cases  28  to  33  in  Table  2,  in  which  i  *  J 
and  n  *  m,  we  use 


R4D  *  “d*  -  V 


J4D  “  2RD  JD 
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where  and  JD  are  obtained  using  Eqs.  (22)  and  (23).  Finally,  the  real  and 
imaginary  parts  of  the  normalized  fourth  moment  are  obtained  using 


Re(M4) 

-  <**>  «4D  *  <V  J«n 

(48) 

Im(M4) 

■  <V  R4D  -  <V  J4D  • 

(49) 

2.7  Extrema 

For  all  four  antennas,  the  maximum  and  minimum  observed  intensities  are 
obtained.  These  extremal  intensities  are  normalized  by  dividing  by  the  mean 
Intensity,  because  the  intensity  is  measured  in  arbitrary  units.  Likewise  the 
maximum  and  minimum  observed  phase  difference  is  obtained  for  each  of  the  six 
antenna  pairs.  These  extremal  phase  differences  are  normalized  by  subtract¬ 
ing  the  mean  phase  difference  of  the  antenna  pair,  because  this  mean  is  an 
arbitrary  offset.  Besides  being  intrinsically  interesting,  these  extrema 
will  take  on  any  greatly  erroneous  values  that  appear  in  the  intensities  and 
phase  differences.  Hence,  the  extrema  reveal  erroneous  data  that  are  Included 
in  all  of  the  other  mn-wave  statistics  described  in  this  report. 

2.8  Probability  Distribution  Functions 

We  calculate  the  probability  distribution  function  (PDF)  of  each  of  the 
four  intensities.  We  have  control  over  the  range  of  intensity  values  included 
in  these  PDFs,  and  foreknowledge  of  the  intensity  extrema  allows  us  to  effi¬ 
ciently  divide  the  legitimate  range  of  observed  values  into  bins,  thereby 
discarding  obviously  erroneous  data.  The  PDFs  are  calculated  as  histograms  by 
adding  unity  to  the  histogram  value  in  a  given  bin  each  time  the  intensity 
lies  within  the  given  bln.  These  histograms  are  normalized  after  the  desired 
amount  of  data  has  been  included  within  them.  In  the  same  manner  we  obtain 
the  six  phase  difference  PDFs,  one  for  each  antenna  pair. 

2.9  Fourier  Transforms 

We  calculate  the  power  spectra  of  each  of  the  four  Intensities  and  six 
phase  differences.  The  computer  programs  allow  us  to  also  calculate  cross- 


.k* i.BWWi I  LfiXVi 


spectra  between  only  two  quanties,  e.g.,  between  intensities  from  different 
antennas,  or  phase  differences  from  different  antenna  pairs,  or  between  an 
intensity  and  phase  difference.  By  cross-spectra  we  mean  cospectra,  quadra¬ 
ture  spectra,  coherence  spectra,  or  phase  spectra  of  two  different  quantities. 

The  fast  Fourier  transform  algorithm  is  used  on  4096-point  time  series, 
corresponding  to  40.96  s  of  data.  A  number  of  these  transforms,  nominally  10 
of  them,  are  averaged  together  to  increase  the  statistical  reliability.  The 
resulting  transforms  are  combined  to  produce  the  aforementioned  power  and 
cross-spectra,  which  are,  in  turn,  each  averaged  to  36  spectral  values  in  36 
approximately  logarithmically  spaced  frequency  bins. 

Fourier  components  of  yet  lower  frequency  are  obtained  by  block,  averaging 
the  time  series  to  4096  points  per  data  tape  (or  per  time-stationary 
Interval).  The  data  tapes  are  roughly  2400  s  in  duration,  so  averages  of 
roughly  0.6  s  suffice.  The  fast  Fourier  transforms  of  these  0.6-s  averaged 
time  series  are  also  combined  to  produce  the  aforementioned  power  and  cross¬ 
spectra  and  are  averaged  to  36  logarithmically  spaced  frequency  bins. 

Finally,  these  low-frequency  spectra  are  combined  with  their  high-frequency 
counterparts  to  form  composite  spectra  having  72  frequency  bins.  There  is 
adequate  spectral  overlap  between  the  lowest  frequency  in  the  high-frequency 
part  and  the  highest  frequency  in  the  low-frequency  part. 

2.10  Figures  of  Merit 

Certain  redundancies  in  the  non-wave  data  make  it  possible  to  calculate 
quantities  that  would  have  known  values  if  the  measurements  were  perfect.  We 
call  these  quantities  figures  of  merit,  since  their  deviations  from  their  ideal 
values  Indicate  the  level  of  precision  in  the  mm-wave  instrumentation. 

The  first  such  figure  of  merit  is  the  p  type  given  by 

p  =  l^V0*ij  “  ^V180^ij^  +  ^'Wij  "  *V270*ijl  }/Xi  Xj  * 

There  is  one  p- type  figure  of  merit  for  each  of  the  six  antenna  pairs.  The 
p  type  is,  in  effect,  the  squared  modulus  of  the  MCF  of  the  non-wave  instru¬ 
mentation.  Ideally  It  would  have  the  value  unity  at  each  instant,  and  its 
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deviations  from  unity  place  a  limit  on  the  accuracy  of  measuring  the  propaga¬ 
tion  MCF  at  each  antenna  separation.  We  calculate  the  mean  and  central 
variance  for  these  six  p-type  figures  of  merit. 

The  second  figure  of  merit  is  the  Z  type  given  by 


l  5 1(Vij 


+  (y  +  (y 
'  90;ij  1  lfto'ij 


(51) 


This  type,  one  for  each  of  six  antenna  pairs,  should  also  be  unity  at  each 
Instant.  We  calculate  the  mean  and  central  variance  of  each  of  the  six  E's. 


The  last  figure  of  merit  is  an  obvious  one  called  the  "phase"  type.  If  i, 
j ,  and  k  are  three  distinct  antennas  such  that  the  k  antenna  lies  between  the 
i  and  j  antennas,  then  the  phase-difference  sum  Bjk  +  Bj^  -  B^j  should  ideally 
be  a  constant  determined  by  the  arbitrary  phase-difference  offsets  of  the 
antenna  pairs.  Except  for  these  offsets  this  phase-difference  sum  would 
Ideally  be  zero.  Because  of  these  offsets,  the  mean  value  of  this  figure  of 
merit  is  meaningless,  but  the  variance  about  the  mean  is  a  measure  of  the  pre¬ 
cision  of  the  phase-difference  measurements.  This  variance  is  the  phase-type 
figure  of  merit.  The  phase-difference  variance  from  each  of  Blk,  ,  and 
Bjj  should  be  much  larger  than  this  figure  of  merit,  otherwise  at  least  one 
of  the  phase  differences  is  suspect.  There  are  four  such  phase-type  figures 
of  merit  when  the  antennas  are  taken  three  at  a  time  as  above,  and  one  when 
all  four  antennas  are  Included.  The  latter  is  the  variance  of  B^  +  ®23  + 

834  -  0j4*  We  calculate  these  five  figures  of  merit. 


3.  QUANTITIES  COMPUTED  FROM  THE  NOAA 
MICROMETEOROLOGICAL  DATA 

3.1  Introduction 

The  NOAA  micrometeorologlcal  data  are  obtained  at  two  towers,  3.6  m  in 
height,  located  adjacent  to  the  propagation  path  at  approximately  one-third 
and  two-thirds  of  the  distance  from  the  ram-wave  receiver  to  the  transmitter. 
These  are  designated  as  stations  1  and  2  respectively,  and  are  shown  on  Fig. 

1.  The  two  optical  paths  were  located  50  m  and  65  m  east  of  the  mm-wave  path. 
The  micrometeorologlcal  towers  were  80  m  east  of  the  ram-wave  path,  and  the 
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optical  rain  gauge/dlsdrometer  paths  were  90  m  east.  Each  tower  has  a  three- 
axis  sonic  anemometer  with  a  platinum  resistance-wire  thermometer  In  the 
center  of  the  sonic  array,  a  nearby  Lyman-a  hygrometer  with  another  wire  ther¬ 
mometer  within  its  gap,  a  prop-vane  anemometer,  and  a  psychrometer  for  mean 
temperature  and  humidity.  The  wire  thermometers  and  Lyman-a  hygrometer  were 
sampled  at  100  Hz.  The  sonic  and  prop-vane  anemometers  were  sampled  at  25  Hz. 
Beginning  with  the  0.39-Hz  pulse,  the  psychrometer  gives  six  commutation 
cycles  repeated  at  the  25-Hz  rate.  A  commutation  cycle  consists  of  data  from 
the  wet  and  dry  bulbs  and  two  reference  values,  four  quantities  in  all.  The 
psychrometer  then  returns  toward  thermal  equilibrium  until  the  next  0.39-Hz 
pulse  calls  for  data;  this  minimizes  self-heating  in  the  psychrometer  due  to 
power  dissipation.  Station  2  has  a  VAistflff  mean  temperature  and  humidity 
instrument  sampled  at  25  Hz,  but  the  Vtflsfflg  readings  of  humidity  are  judged 
inferior  to  those  of  the  psychrometer.  Station  2  included  an  up-looking  pyra- 
nometer  for  solar  radiation  and  cloud-shadow-occurrence  measurements,  and  sta¬ 
tion  1  included  a  barometer;  these  were  sampled  at  25  Hz.  A  combined  optical 
Cn  and  optical  cross-wind  meter  was  operated  on  the  path  from  near  the  mm- 
wave  transmitter  to  station  1  and  another  from  near  the  mm-wave  receiver  to 
station  2;  these  optical  paths  overlap.  Optical  rain  gauges  were  deployed 
near  both  stations  and  six  weighing  bucket  rain  gauges  were  deployed  along  the 
path.  At  station  1  the  optical  rain  gauge  included  a  disdrometer.  The  data 
from  these  optical  instruments  were  sampled  at  0.39  Hz.  The  reduction  of  the 
data  from  the  rain  instruments  is  not  discussed  in  this  report. 

Additional  optical  instrumentation  was  deployed  during  the  fourth  experl- 

2 

ment  session.  A  10.6-pm  CO.  laser,  for  measurement  of  infrared  C  in  clear 

2  n 

air  and  for  phase-difference  spectra  in  rain,  was  operated  from  near  the 

mm-wave  transmitter  to  inside  the  mm-wave  receiver  trailer.  A  third  large- 

2 

aperture  optical  CR  and  cross-wind  meter  was  operated  over  the  northern  two- 
thirds  of  the  range  directly  below  its  existing  counterpart.  The  purpose  was 
to  measure  vertical  shear  of  the  path-averaged  optical  cross  wind.  A  He-Ne 
laser  scintillometer  and  LED  scintillometer  were  deployed  over  the  same  short 
path  to  measure  the  inner  scale  of  turbulence  by  its  effect  on  optical  scin¬ 
tillation.  At  station  1  an  up-looking  photocell  was  installed  to  compare 
cloud  cover  measurements  with  those  of  the  up-looking  pyranometer  at  station 
2.  In  addition,  another  psychrometer  was  installed  at  station  2  to  determine 
the  Instrument's  reliability  by  comparison  with  its  existing  counterpart. 


GxcepC  for  the  CC>2  laser,  these  additional  Instruments  were  recording  data 
during  the  second  half  of  the  fifth  session  as  well. 

The  heights  of  the  instrumentation  as  well  as  subjective  estimates  of  the 
3 a  errors  in  the  heights  are  as  follows:  The  mm-wave  beam  and  CC^-laser  beam 
average  3.68  m  ±  0.1  m;  the  Lyman-a  gap,  horizontal  sonic  transducers,  prop- 
vane,  and  platinum  wires  are  at  3.61  m  ±  0.05  m;  the  psychrometer  inlets  are 
at  1.91  m  ±  0.05  m;  the  upper-level  optical  wind  and  Cn  meters  and  the  opti¬ 
cal  inner-scale  meter  are  at  3.78  m  ±  0.1  m;  and  the  lower-level  optical  wind 
2 

and  Cn  meter  are  2  m  ±  0.02  m  below  the  upper-level  meters. 

3.2  Conversion  to  Units 

Unlike  the  mm-wave  data,  the  micrometeorological  data  must  be  converted  to 
physical  units  like  degrees  Celsius  or  meters  per  second.  This  conversion 
depends  on  mean  temperature  for  the  sonic  anemometers  and  on  mean  humidity  for 
the  Lyman-a  hygrometers.  The  0.39-Hz  data  is  converted  to  units  giving  a 
datum  from  each  instrument  every  2.56  s.  The  middle  four  of  the  six  com¬ 
mutation  cycles  from  the  psychrometer  are  averaged  to  provide  wet-  and  dry- 
bulb  temperatures  at  the  0.39-Hz  rate.  The  temperatures  from  the  platinum 
resistance-wire  thermometers  and  the  voltage  counts  from  the  Lyman-a  hygrome¬ 
ters  are  averaged  over  2.56-s  bins  to  determine  their  trends.  Comparison  of 
these  Lyman-a  voltage  counts  averaged  over  the  entire  tape  with  the  similarly 
averaged  psychrometer  humidity  establishes  the  Lyman-a  calibration. 

The  resistance-wire  signals  are  converted  to  degrees  Celsius  and  Lyman-a 
signals  are  converted  to  absolute  humidity  in  grams  of  water  vapor  per  cubic 
meter,  all  at  their  100-Hz  data  rate.  The  prop-vane  signals  are  converted  at 
their  25-Hz  data  rate  to  wind  speed  in  meters  per  second  and  angle  in  degrees. 
This  wind  angle  is  zero  if  the  wind  is  perpendicular  to  the  front  of  the 
instrument  boom  and  increases  if  the  vane  turns  clockwise  as  viewed  from 
below.  The  wind  components  along  the  sonic  axes  are  converted  to  meters  per 
second,  at  their  25-Hz  rate.  The  solar  pyranometer  output  is  expressed  as  a 
fraction  of  full  sun  (1.395  W/m2)  at  its  25-Hz  rate.  The  barometer  signal  is 
converted  to  millibars  at  its  25-Hz  rate. 

The  angle  between  the  horizontal  sonic  axes  is  120°  to  allow  a  greater 
acceptance  of  wind  angle  variation.  Consequently  the  measured  horizontal 


sonic  wind  components  are  not  orthogonal.  Our  sonic  configuration  is  shown  by 
Kairaal  (1980a;  Figs.  2  and  3).  Orthogonal  components  are  obtained  from  the 
nonorthogonal  ones  such  that  the  u-coraponent  is  perpendicular  to  the  instru¬ 
ment  boom  and  is  positive  when  the  wind  is  toward  the  boom  and  the  v-component 
is  parallel  to  the  boom,  and  is  positive  when  the  wind  is  blowing  right  to 
left  of  a  person  facing  the  boom.  If  we  let  and  Vg  denote  the  wind  com¬ 
ponents  along  the  "A"  and  "B"  sonic  axes  as  depicted  by  Kaimal  (1980a;  Figs.  2 
and  3) ,  then  u  and  v  are  given  by 


“  ■  VA  +  VB 

(52) 

V  -  (VB  -  Va)//3  . 

(53) 

The  vertical  axes  of  the  sonic  anemometer  are  maintained  to  within  0.1°  of 
vertical  so  that  the  measured  vertical  wind  component,  w,  is  negligibly  con¬ 
taminated  by  horizontal  wind  fluctuations.  The  u,  v,  and  w  components  then 
form  a  right-handed  orthogonal  triad.  Although  the  boom  is  rotated  into  the 
wind  at  the  beginning  of  each  data  run,  there  remains  a  slight  difference  bet¬ 
ween  our  notation  for  u,  v,  w  and  the  standard  micrometeorologlcal  notation 
that  the  u-component  is  along  the  time-averaged  direction  of  the  wind  (i.e., 
it  is  the  streamwise  component)  and  that  v  is  the  cross-stream  component. 

Angle  encoders  at  both  meteorological  towers  give  the  angle  of  the  instru¬ 
ment  boom  relative  to  geographic  directions.  These  encoders  give  angles 
changing  from  0°  to  90°,  to  180®,  to  270°  as  the  boom  rotates  from  facing 
north  to  facing  east,  to  facing  south,  to  facing  west  (N,  E,  S,  W) ,  respec¬ 
tively.  (As  a  result  of  an  error,  these  encoder  angles  are  for  directions  W, 
N,  E,  S  for  the  first  two  data  tapes  of  the  third  experiment  session.)  We 
define  the  geographic  wind  angle  as  the  encoder  angle  minus  the  prop-vane 
angle.  As  mentioned  before,  the  vane  angle  is  relative  to  the  boom;  it 
Increases  clockwise  when  one  looks  up  at  the  boom  from  below,  and  is  zero  for 
wind  blowing  into  the  face  of  the  boom.  The  geographic  wind  angle  is  then  0®, 
90®,  180®,  and  270®  for  wind  blowing  from  the  N,  E,  S,  and  W,  respectively. 
(Except  for  the  first  two  data  tapes  of  the  third  experiment  session  when 
these  angles  are  for  wind  from  the  W,  N,  E,  and  S,  respectively.)  The  optical 

cross  wind  telemetered  from  station  1  is  obtained  over  the  southern  two-thirds 

of  the  mm-wave  range  and  is  positive  for  wind  from  the  east.  The  optical 

cross  wind  telemetered  from  station  2  is  obtained  over  the  northern  two  thirds 

of  the  range  and  is  positive  for  wind  from  the  west. 


3.3  Averages,  Variances,  and  Cross-Correlations 


Table  3  lists  the  averages  performed  la  each  computation.  An  average  Is 
denoted  by  the  angle  brackets.  The  symbols  for  the  quantities  are  as  follows: 
q  Is  the  Lyman-a  humidity  in  grams  per  cubic  meter;  T^  is  the  temperature  in 
Celsius  from  the  resistance  wire  within  the  Lyman-a  gap;  Tg  is  the  temperature 
in  Celsius  from  the  resistance  wire  within  the  sonic  anemometer;  u,  v,  w  are 
the  orthogonalized  sonic  anemometer  wind  components  in  meters  per  second  as 
discussed  in  Sec.  3.2;  V  and  9  are  the  prop-vane  wind  speed  in  meters  per 
second  and  angle  in  degrees,  respectively;  S  is  the  pyranometer  solar  flux. 

Table  3. — Kinds  of  averages  performed  on  the 
NOAA  micrometeorological  data 


Averages 

<q>,  <Tl>,  <Ts>,  <u>,  <v>, 

<w>,  <V>,  <0>,  </u2  +  v2>,  <S> 

Mean  squares 

<q2>,  <Tl2>,  <Ts2>,  <u2>,  <v2>, 

<w2>,  <V2>,  <02>,  <u 2  +  v2>,  <S 

Cross-correlations 

<u  w>,  <v  w>,  <u  v>,  <q  w>,  <Tl 

<Ts  w>,  <q  Tl>,  <q  Ts>,  <Tl  Tg> 

</u2  +  v2  w> 

Note:  See  text  for  an  explanation  of  the  symbols. 


The  averages  listed  in  Table  3  are  obtained  for  both  micrometeorological  sta¬ 
tions.  In  the  case  of  cross  products  like  <Tg  w>,  which  have  a  100-Hz  quan¬ 
tity,  Tg,  multiplied  by  a  25-Hz  quantity,  w,  we  average  four  consecutive 
values  of  the  100-Hz  quantity  before  performing  the  product  and  including  the 
term  in  the  summation. 

The  mean  squares  and  cross-correlations  given  in  Table  3  are  not  of 
interest  in  themselves.  After  the  averaging  is  completed  the  mean  squares  are 
converted  to  variances  taken  about  the  mean  (centralized  variances)  and  the 
cross-correlations  are  likewise  converted.  The  following  equations  give  the 
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desired  central  moment  on  the  left  side  and  the  required  calculation  on  the 
right  side.  For  the  variances,  where  Q  may  be  w,  u,  v,  V,  Tg,  T^,  or  q,  we  use 

<(Q  -  <Q»Z>  -  <Q2>  -  <Q>2  .  (54) 

For  the  centralized  cross-correlations,  where  the  pair  Qj »Q£  may  be  u,w  or 
v,w  or  u,v  or  q,w  or  T^,w  or  Tg,w  or  TL,q,  we  use 

<(Qt  -  <Q1»(Q2  -  <Q2»>  “  <Qj  Q2>  -  <Q!><Q2>  •  (55) 

We  subtract  the  mean  of  quantities  like  w,  which  should  have  a  zero  mean 
value,  because  instrumental  offsets  cause  the  measured  values  to  have  nonzero 
means,  and  these  offsets  must  be  subtracted. 


3.4  Turbulence  Fluxes 

All  of  the  quantities  discussed  in  this  section  are  obtained  from  averages 
previously  performed.  The  correlations  discussed  in  the  previous  section  form 
the  basis  for  the  calculation  of  the  fluxes. 

The  quantities  <(Tg  -  <Tg>)(w  -  <w>)>  and  <(q  -  <q>)(w  -  <w>)>  are  the 
vertical  fluxes  of  temperature  and  humidity,  respectively.  Temperature  flux 
can  be  converted  to  heat  flux  by  multiplying  by  the  constant-pressure  heat 
capacity  of  air,  and  humidity  flux  can  be  converted  to  latent  heat  flux  by 
multiplying  by  the  latent  heat  of  vaporization.  The  above  temperature  flux  is 
obtained  by  correlating  the  fluctuating  temperature  from  the  resistance-wire 
thermometer  within  the  sonic  array  with  the  fluctuating  vertical  velocity  from 
the  sonic  anemometer.  The  Lyman-a  hygrometer  is  spatially  separated  from  the 
sonic  array  so  that  it  does  not  disturb  the  turbulent  wind  field  at  the  array. 
The  spatial  separation  slightly  reduces  the  humidity  flux  from  its  true  value. 
However,  we  also  calculate  the  temperature  flux  using  the  resistance  wire 
within  the  Lyman-a  gap  from  the  cross-correlation  <TL  w>.  The  ratio  of  the 
two  temperature  fluxes, 

<(TL  -  <TL»(w  -  <w»>/<(Tg  -  <Tg> )(w  -  <w>)>  ,  (56) 

is  a  correction  factor  for  this  spatial  separation.  Because  of  the  high 
correlation  and/or  anticorrelation  between  temperature  and  humidity  fluc- 
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tuations,  the  humidity  flux  can  be  multiplied  by  this  ratio  to  correct  it  for 
the  spatial  separation. 


To  obtain  the  vertical  flux  of  horizontal  momentum  we  must  cross-correlate 


fluctuating  vertical  wind  speed  with  the  fluctuating  horizontal  wind  component 
along  the  time-averaged  wind  direction  (the  streamwise  component).  Given  the 
mean  wind  angle  <9>  from  either  the  vane  or  sonic  anemometer,  the  streamwise 
component  Ug  and  cross-stream  component  Vg  are  given  in  terms  of  the  u  and  v 
components  by 


Ug  =«  u  cos<9>  +  v  sin<9> 


v.  »  -u  sin<9>  +  v  cos<9>  . 

O 


Letting  primes  denote  zero-mean  fluctuations,  e.g.,  u  '  -  u.  -  <u_>  and  v' 

0  tj  o 


v  -  <v> ,  then  we  have 


Ug'  *  u'  cos<9>  +  v'  sin<0> 


v0’  ■  -u'  sin<9>  +  v'  cos<8>  . 


The  vertical  momentum  flux  is  given  by  multiplying  the  quantity 


<uQ'  w’>  «  <u'  w’>  cos<9>  +  <v'  w*>  sin<0> 


by  the  mass  density  of  the  air.  Calculation  of  the  quantities 


<u'  w'>  »  <u  w>  -  <u><w> 


<v'  w’>  •  <v  w>  -  <v><w> 


has  already  been  discussed  in  the  previous  section.  The  vertical  flux  of 


cross-stream  momentum 


<Vg'  w*>  »  -  <u*  w’>  sin<9>  +  <v'  w'>  cos<9>  , 


and  the  horizontal  flux  of  momentum  obtained  from 
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<V  V> 


/2(<v'2>  -  <u,2»  sin(  2<8>)  +  <u'  v’>  cos(2<9>) 


(65) 


should  both  be  zero  for  a  horizontally  homogeneous  surface  layer,  but  they  are 
computed  anyway.  In  addition,  for  station  2  we  calculate  the  speed  flux  given 
by 


-  </u‘ 


v2>^(w  - 


<w>  )> 


We  can  calculate  the  momentum  flux  using  (61)  in  a  variety  of  ways 

depending  on  how  we  choose  <0>.  We  could,  for  instance,  take  <9>  to  be  the 

average  of  the  vane  angles.  The  conventional  method  is  to  set  the  average 

cross-stream  component  to  zero,  i.e.,  <v  >  -  0,  so  from  (58)  <9>  =* 

0 

arctan(<v>/<u>)  as  given  by  Kaimal  (1980a).  However,  this  <9>  is  sensitive 
to  the  poorly-known  offset  voltages  of  the  sonic  anemometer.  On  the  other 
hand,  we  can  assume  a  horizontally  homogeneous  surface  layer  and  set  <v  '  w’> 

“  0  In  (64)  to  define  <0>;  this  definition  of  <8>  is  independent  of  sonic  off¬ 
set  voltages.  Using  this  latter  <9>,  (61)  gives  for  the  momentum  flux 


<u  '  w’>  *  /<u'  w’>2  +  <v '  w’>2  . 

This  result  is  independent  of  sonic  offset  voltages,  but  depends  on  horizontal 
homogeneity.  We  calculate  this  latter  momentum  flux,  but  we  also  use  the 
mean  vane  angle  to  calculate  the  expressions  in  (61),  (64),  and  (65). 


3.5  Spatial  Structure  Parameters 

The  spatial  separation  of  the  two  resistance-wire  thermometers  and  of  the 
Lyman-»a  hygrometer  and  the  resistance  wire  within  the  sonic  anemometer  enable 
us  to  calculate  temperature  structure  parameters  and  temperature-humidity 
cross-structure  parameters.  The  temperature  structure  parameter  is  given  by 

°T2  *  <(V  -  V)2>/rSL2/’  (67> 

-  «Tg2>  -  <Tg>2  +  <Tl2>  -  <Tl>2  -  2<T§  Tl>  +  2<TgXTL>)/rSL2/3  (68) 


where  primes  denote  zero-mean  fluctuations  and  rg^  is  the  distance  between  the 
resistance-wires  in  the  sonic  anemometer  and  the  Lyman-a  gap.  The  second  for¬ 
mula  (68)  shows  that  this  structure  parameter  can  be  calculated  from  quan¬ 
tities  already  discussed  in  Sec.  3.3. 

If  we  had  additional  Lyman-a  hygrometers  positioned  within  the  sonic  ane¬ 
mometer  yielding  fluctuating  humidities  qg'»  then  we  could  calculate  the 
temperature-humidity  cross-structure  parameter  defined  by 

CT,  ’  <(V  -  V><’s'  -  ’')>/rSl.2/3  •  «»> 

However,  if  we  assume  that  this  hypothetical  hygrometer  would  give  <T  '  q  '>  • 

<Tl'  q’>  and  <Tg'  q’>  *  <TL'  qg'>  (the  latter  requires  horizontal  Isotropy), 
then  we  have 

CTq  -  2(<Tl*  q’>  -  <Tg*  q*>)/rSL2/3  (70) 

-  2«Tl  q>  -  <TLXq>  -  <Tg  q>  +  <Tg><q»/rSL2/3  ,  (71) 

where  (71)  shows  that  the  calculation  obtains  from  quantities  already 
discussed  in  Sec.  3.3. 


3.6  Time<Lagged  Structure  Paraucters 

2  2 

The  time-lagged  structure  parameters  of  temperature,  C  ,  humidity,  C_  , 

2  ^ 

C  ,  and  the  temperature-humidity  cross-structure  parameter,  C_  ,  are  defined 


CT2  -  <(Tg(t)  -  Tg( t  +  t)J2>/(<V>t)2/3 

s 

CT2  -  <[TL(t)  -  TL(t  +  t) ] 2>/(<V>t)2/3 
L 

C  2  -  <[q(t)  -  q(t  +  t)]2>/(<V>t)2/3 
<1 


CT  q  "  <[TL(t)  ‘  Vc  +  t)Hq(t)  -  1<t  «■  t)]>/«V>t)4 
L 


where  t  is  time,  r  is  time  lag,  and  <V>  is  the  average  wind  speed.  The  data 
are  digitized  at  100  Hz,  so  t  must  be  a  multiple  of  0.01  s.  Let  the  height 
above  ground  be  denoted  by  Z,  then  Z/5  should  be  in  the  Inertial  range  of 


v.  v.  r.  tr*  «  ■  - 


sizes.  We  define  a  dimensionless  factor  by 

F  =  (Z/5X100  samples/s)/<V>  .  (76) 

We  then  define  three  integers  *  F,  J2  ■  F/2,  and  ■  F/4,  where  the 
integers  are  rounded  up  from  the  real  numbers.  These  integers  are  the  number 
of  0.01-s  data  samples  to  skip  to  calculate  the  above  four  structure  parame¬ 
ters  for  both  stations  at  three  values  of  time  lag.  The  time  lags  used  are 
thus  *  Jj^O.Ol  s),  Tj  *  ^(0.01  s),  and  *  J^CO.Ol  s);  these  correspond 
approximately  to  spatial  lags  of  Z/5,  Z/10,  and  Z/20.  The  integers  and  lags 
can  be  different  at  the  two  stations  because  <V>  can  differ  between  the  sta¬ 
tions. 

Of  course,  the  time-lagged  structure  parameters  and  cross-structure  para¬ 
meters  can  also  be  obtained  from  the  spectral  level  of  the  appropriate  power 
spectra  and  cospectra,  which  are  discussed  in  Sec.  5.9. 


3.7  Extrema 

We  obtain  the  minimum  and  maximum  observed  values  of  several  quantities 
for  each  averaging  period.  These  quantities  are  the  wind  speed  and  wind  angle 
from  the  prop-vane  anemometer  and  the  solar  flux  from  the  pyranometer. 

3.8  Probability  Distribution  Functions 

The  probability  distribution  functions  (PDFs)  are  calculated  for  the  pyra¬ 
nometer  solar  flux  and  the  wind  angles  from  the  two  prop-vane  anemometers. 

The  joint  probability  distribution  functions  (JPDFs)  of  Lyraan-o  humidity  and 
its  associated  resistance-wire  temperature  is  calculated  for  both  stations. 

The  PDFs  of  the  fluctuating  humidity  or  temperature  can  be  obtained  by  summing 
the  JPDF  over  all  temperature  or  humidity  bins,  respectively.  These  PDFs 
obtained  from  JPDFs  are  sometimes  called  marginal  PDFs.  Since  for  clear  air 
the  fluctuations  of  real  and  imaginary  parts  of  the  refractive  index  depend 
almost  entirely  on  only  temperature  and  humidity  fluctuations,  we  can  obtain 
the  JPDF  of  these  real  and  Imaginary  parts  from  the  temperature-humidity  JPDF. 
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Of  course,  the  marginal  PDFs  of  refraction  and  absorption  can  be  obtained  from 
this  latter  JPDF  as  well. 


* 


The  method  of  calculating  these  PDFs  and  JPDFs  as  histograms  using 
foreknowledge  of  the  extrema  to  divide  efficiently  the  data  into  bins  is  the 
same  as  was  discussed  for  the  mm-wave  PDFs  in  Sec.  2.8. 

3.9  Fourier  Transforms 

We  calculate  the  power  spectra  of  the  fluctuating  temperature  and  humidity 
from  the  resistance-wires  and  the  Lyman-ot  hygrometer,  as  well  as  the 
temperature-humidity  cospectra  and  quadrature  spectra  from  these  instruments. 
We  also  calculate  the  power  spectra  of  all  three  orthogonalized  components  of 
the  fluctuating  wind  from  the  sonic  anemometers,  as  well  as  the  cospectra  and 
quadrature  spectra  of  these  components  taken  in  pairs.  The  cospectra  and 
quadrature  spectra  are  calculated  for  temperature  and  vertical  velocity,  and 
for  humidity  and  vertical  velocity. 

The  fast  Fourier  transform  (FFT)  algorithm  is  used  on  4096-point  time 
series.  This  corresponds  to  40.96  s  of  data  if  only  100-Hz  quantities  are 
involved  or  163.84  s  of  data  if  at  least  one  quantity  is  of  the  25-Hz  sampling 
rate.  A  number  of  these  FFTs  are  averaged  together,  nominally  10  of  them,  to 
Increase  the  statistical  reliability.  As  with  the  mm-wave  Fourier  transforms, 
the  FFTs  are  combined  to  produce  power  spectra  and  cross-spectra  that  are 
similarly  averaged  into  29  logarithmically  spaced  frequency  bins. 


3.10  Refraction  and  Absorption  Structure  Parameters 

The  zero-mean  fluctuations  in  the  real  and  Imaginary  parts  of  the  refrac¬ 
tive  index,  denoted  n^'  and  n^.',  respectively,  are  related  to  those  of  tem¬ 
perature  and  humidity  by 
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where  <T>  is  the  mean  Kelvin  temperature,  the  effects  of  pressure  fluctuations 
are  neglected,  and  these  equations  are  valid  for  fluctuations  much  smaller 
than  the  means.  The  dimensionless  sensitivity  coefficients  A^.,  A^,  B^.,  and 
Bq  depend  on  mean  temperature,  humidity,  pressure,  and  the  electromagnetic 
wave  length;  they  are  calculated  from  the  radio-frequency  refractive-index 
equation,  as  well  as  the  sum  of  dispersion  and  absorption  contributions  from 
water  vapor  and  oxygen  resonances,  oxygen  and  nitrogen  nonresonant  absorption, 
and  the  water  vapor  anomalous  absorption.  The  computer  programs  for  this 
calculation  are  not  discussed  here. 

Equations  (77)  and  (78)  provide  statistics  of  refraction  and  absorption 
from  those  of  temperature  and  humidity.  For  instance,  the  refractive-index 
structure  parameters  are  given  by 
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If  these  structure  parameters  are  replaced  by  their  corresponding  power 
spectra  and  cospectra,  then  these  equations  give  the  power  spectra  and 
cospectt'a  of  refraction  and  absorption  in  terms  of  those  of  temperature  and 
humidity.  Likewise,  by  analogous  replacement,  these  equations  express  the 
variances  and  cross-correlations  of  refraction  and  absorption  in  terms  of 
those  of  temperature  and  humidity. 


Clearly,  all  the  relevant  statistics  of  refraction  and  absorption  fluctua¬ 
tions  can  be  obtained  from  the  already-calculated  mlcrometeorological  sta¬ 
tistics  once  the  sensitivity  coefficients  are  known.  In  particular,  the 
linear  relationships  (77)  and  (78)  make  it  easy  to  calculate  the  joint  proba¬ 
bility  distribution  function  of  n^'  and  n^'  from  that  of  fluctuating  tem¬ 
perature  and  humidity. 


4.  PROBLEMS  WITHIN  THE  DATA 

4.1  Introduction 

The  greatest  fraction  of  personnel  time  expended  on  the  data  processing 
software  has  been  used  to  fix  problems  In  the  data.  The  discovery  and  diagno¬ 
sis  of  these  problems  has  been  challenging  and  time  consuming.  During  the 
July  1983  session,  the  data  aquisltlon  system  dropped  occasional  bytes, 
dropped  data  frames  of  0.01-s  duration  as  well  as  entire  blocks  of  data  of 
0.16-s  duration,  and  wrote  parity  errors  on  the  data  tapes.  The  aquisltlon 
system  was  modified  for  subsequent  sessions.  The  mm-wave  data  contain  the 
expected  data  dropouts  during  periods  of  phase-lock  loss,  but  they  also  con¬ 
tain  occasional  spikes,  step  changes  in  signal  gain,  harmonic  noise  due  to 
vibration  of  the  steel  I-beam  holding  the  receivers,  and  other  noise  of 
unknown  origin.  During  the  first  session,  the  intensity  at  antenna  4  is 
Incorrect.  The  micrometeorological  data  had  several  Instrument  failures 

during  that  session.  The  Lyman-a  hygrometer  at  station  1  was  not  working;  one 

2 

of  the  optical  wind  and  meters  was  not  working;  the  resistance-wire  ther¬ 
mometer  signals  exhibited  cold  spikes  and  hot  spikes  that  are  characteristic 
of  contamination  of  the  wire  by  hydroscopic  particles;  and  the  sonic  anemo¬ 
meter  signals  contained  numerous  spikes  (one  per  several  seconds  per  channel). 
Most  of  these  problems  have  been  solved  by  clever  but  time-consuming  software 
development. 

4.2  Cures  for  Data  Aquisltlon  Errors 

A  data  frame  is  the  digitized  data  from  all  Instruments  In  each  one- 
hundredth  of  a  second.  The  first  byte  In  each  frame  is  synchronization  (sync) 
information.  It  was  found  that  some  of  these  sync  bytes  are  missing  from  the 
data  tape.  Since  the  sync  bytes  should  arrive  with  every  frame,  it  was 
possible  to  insert  the  missing  byte  when  the  data  tapes  were  rewritten  for  our 
data  processing  computer. 

Entire  data  frames  were  found  to  be  missing  with  no  gap  in  frame  sync 
bytes.  The  first  experiment  session  had  a  10-Hz  clock  written  to  tape,  so 
precise  frame  loss  detection  is  Impossible.  The  25-Hz  data  are  written  to 
tape  by  repeating  the  same  value  four  times  for  each  25-Hz  channel;  missing 


data  frames  manifest  themselves  by  fewer  than  four  repetitions  of  25-Hz  data. 

A  special  flag  is  set  to  unity  in  the  first  data  frame  of  this  25-Hz  repeti¬ 
tion  and  sec  to  zero  for  the  last  three  frames  in  the  repetition.  Using  this 

special  flag  it  is  possible  to  identify  how  many  frames  are  missing,  assuming 

not  more  than  three  frames  were  lost.  We  then  replace  the  missing  frame  or 
frames  by  duplicating  the  frame  preceding  the  occurrence  of  lost  frames  and 
setting  the  special  flags  when  necessary.  This  restores  the  integrity  of  the 
25-Hz  and  0.39-Hz  data  and  the  temporal  continuity  of  the  data.  Beginning 
with  tape  12  of  session  2  through  session  5,  the  lost  frame  detection  is  much 

easier  since  the  100-Hz  clock  channel  shows  each  frame  that  is  out  of 

sequence.  Changes  in  the  data  aquisition  system  for  runs  after  the  tape  12  of 
session  2  have  apparently  stopped  the  frame-loss  problem. 

The  data  acquisition  system  writes  data  to  magnetic  tape  one  block  at  a 
time;  a  block  consists  of  16  frames,  i.e.,  0.16  s  of  data.  During  session  1, 
if  the  read-af ter-write  check  failed,  then  the  block  was  rewritten  on  tape. 

The  time  consumed  in  the  rewrite  caused  the  next  block  of  data  to  be  lost. 

Thus  entire  data  blocks  are  missing  from  session  1  data  tapes.  The  rewrite 
was  disabled  for  subsequent  experiment  sessions  to  avoid  this  block-loss 
mechanism.  We  cannot  recreate  lo&t  blocks.  However,  their  approximate  posi¬ 
tions  on  the  data  tapes  can  be  Identified  using  the  system  time  (accurate  to 
0.1  s  In  session  1  and  0.01  s  In  subsequent  sessions)  that  Is  written  Into 
each  frame.  The  positions  of  these  lost  blocks  are  written  to  the  "pedigree" 
disk  file  so  that  these  data  gaps  can  be  avoided  during  Fourier  transfor¬ 
mation.  So  far,  no  block  losses  have  been  found  following  changes  in  the  data 
acquisition  system  that  were  implemented  just  before  tape  12  of  session  2. 

The  data  aquisition  system  wrote  parity  errors  onto  the  data  tapes. 

Roughly  one-half  of  all  data  tapes  contain  a  parity  error.  By  turning  off  the 
error  processing,  we  can  read  the  data  tape  In  an  unsafe  mode.  To  date  this 
has  been  tried  on  only  two  data  tapes.  The  problem  of  the  parity  errors 
exists  in  all  five  experiment  sessions. 

4.3  Cures  for  Millimeter-Wave  Data  Problems 

All  28  mm-wave  channels  require  a  correction  as  a  result  of  nonlinearity 
of  the  detectors.  (The  crystal  detectors  in  the  signal  combiner  are  not 


28 


perfect-square-law  detectors  in  the  radio-frequency  power  regime  used.)  The 
extent  of  nonlinearity  was  measured  for  each  of  the  28  channels  and  fitted  by 
a  polynomial.  In  the  data  analysis  each  datum  for  each  mm-wave  channel  and 
for  every  0.01  s  was  corrected  by  evaluating  a  fourth-order  polynomial 
(roughly  7  million  polynomial  evaluations  per  data  tape).  One  set  of  polyno¬ 
mial  coefficients  is  used  for  session  1,  another  set  for  session  2  through 
tape  10  of  session  3,  and  yet  another  set  thereafter.  This  is  necessary 
because  different  RF  input  power  was  used. 

The  intensity  at  antenna  4,  1^,  was  not  measured  during  session  l  because 

of  a  hard-wiring  error.  Fortunately,  the  redundancy  in  the  mm-wave  data 

2 

allows  us  to  calculate  I,  as  follows.  If  p  in  Eq.  (50)  is  obtained  from  an 

4  2 

antenna  pair  Including  antenna  4,  then  setting  p  -  1  allows  us  to  calculate 

I,  from  the  other  channels.  Likewise  setting  E  *  1  in  Eq.  (51)  gives  the  same 

result.  We  wrote  a  separate  computer  program  to  evaluate  these  methods  of 

synthesizing  an  intensity.  The  program  compared  the  normalized  variance  of 

the  three  intensities  1^,  I2>  and  1^  with  their  synthesized  values  from  every 

relevant  antenna  pair  as  well  as  averages  of  the  synthetic  Intensities  over 

every  combination  of  such  antenna  pairs.  We  determined  that  setting  E  ■  1 

2 

produced  better  results  than  setting  p  *  1  and  that  averaging  the  synthesized 
Intensities  improved  the  reliability.  For  the  test  case  the  difference  bet¬ 
ween  normalized  variances  from  measured  and  synthesized  intensities  ranged 
from  13.5%  to  a  few  tenths  of  a  percent  for  the  E  *  1  method,  and  the  nor¬ 
malized  variance  from  the  synthesized  1^  was  17%  smaller  than  that  from  Ij, 

Ij,  and  Ij.  There  is  no  guarantee  that  such  good  results  will  be  obtained 
from  all  of  session  1  data  tapes.  At  present  we  use  the  E  *  1  method  to 
synthesize  three  values  of  1^  from  the  three  antenna  pairs  containing  antenna 
4  and  average  these  three  1^  values  together  to  produce  1^  for  each  0.01-s 
data  frame. 

The  mm-wave  data  contain  trends  that  have  no  correspondence  to  meteorolo¬ 
gical  trends,  nor  do  the  trends  in  one  ram-wave  channel  appear  related  to 
trends  in  other  channels.  These  trends  are  apparently  caused  by  drift  in 
receiver  gains.  Consequently  we  filter  the  four  intensities  and  six  phase 
differences  using  single-pole  digital  low-pass  filters  having  a  roughly  125-s 
time  constant.  The  filter  value  is  subtracted  from  whe  phase  differences  but 
divided  into  the  intensities.  We  retain  the  unfiltered  mean  value  of  the  four 


Intensities.  (The  figures  of  merit  are  obtained  from  unfiltered  data.) 
However,  this  filtering  requires  that  we  discard  about  5  min  of  data  from  the 
beginning  of  each  data  tape  to  initialize  the  filter. 

The  mm-wave  data  contain  the  expected  periods  of  phase-lock  loss.  These 
are  of  variable  duration  but  typically  less  than  20  s.  There  are  a  variety  of 
types  of  these  data  dropouts.  The  mildest  type  appears  to  have  intact  data, 
but  the  intensity  drops  to  half  its  usual  value.  The  most  severe  type  reduces 
all  mm-wave  channels  to  noise.  We  detect  these  phase-lock-loss  periods 
visually  using  computer-generated  microfilm  stripcharts.  We  note  the 
beginning  and  ending  times  of  each  such  loss  on  these  stripcharts  and  enter 
them  in  the  parameters  file  for  use  in  subsequent  processing.  Thus  we  can 
exclude  the  phase-lock  losses  from  subsequent  processing.  This  greatly 
complicates  our  programming  because  we  must  archive  the  beginning  and  ending 
times,  determine  what  averaging  intervals  contain  such  a  loss,  determine 
whether  or  not  the  loss  extends  across  the  boundary  of  an  averaging  interval, 
keep  account  of  the  number  of  samples  averaged  in  each  interval,  avoid  the 
intervals  in  the  high-frequency  FFT,  and  replace  the  data  with  fake  values  for 
both  the  digital  filtering  and  low-frequency  FFT. 

The  mm-wave  data  also  occasionally  contain  spurious  spikes  of  unknown 
origin  in  one  or  several  channels.  Sometimes  these  spikes  come  in  groups  of 
many  randomly  spaced  spikes.  We  identify  these  spikes  from  the  microfilm 
stripcharts  and  note  the  time  Interval  containing  the  spike.  We  then  treat 
the  interval  as  a  phase-lock  loss  and  exclude  all  mm-wave  data  in  this  inter¬ 
val  from  the  data  processing.  This  creates  the  same  complications  for  our 
processing  as  do  the  phase-lock  loss  periods.  A  sophisticated  "despiker" 
algorithm  could  be  used  to  salvage  more  data. 

On  some  data  tapes,  one  or  more  mm-wave  channels  will  exhibit  a  step 
change  occurring  in  a  time  less  than  0.01  s.  Later,  other  such  step  changes 
may  occur,  either  in  the  same  direction  or  the  opposite.  At  present  we  are 
not  processing  the  data  tapes  containing  this  type  of  error. 

Most  of  the  mm-wave  data  contain  noise.  One  source  of  noise  is  harmonic 
vibration  of  the  steel  I-beam  that  supports  the  mm-wave  receivers.  This  noise 
tends  to  come  and  go  and  may  be  generated  by  wind  gusts  stressing  the  receiver 
semitrailer,  which  transmits  the  stress  to  the  ground  and  then  to  the  I-beam 


supports*  The  differential  motion  of  parts  of  the  I-beam  toward  and  away  from 
the  transmitter  causes  an  artificial  phase  difference.  The  most  common  and 
largest  amplitude  of  this  type  of  noise  is  at  roughly  5  Hz.  Higher  frequen¬ 
cies  have  also  been  observed.  A  form  of  noise  in  the  phase  differences  that 
does  not  seem  to  be  caused  by  I-beam  vibration  is  a  roughly  24-Hz  oscillation 
modulated  by  a  roughly  3-Hz  oscillation  giving  a  popcorn-on-a-string  effect. 
This  noise  can  cover  entire  data  runs  and  persist  from  one  run  to  another. 

Even  higher  frequency  noise  sometimes  appears  in  intensities  as  well  as  in 
phase  differences.  The  mm-wave  Fourier  transform  is  used  to  identify  the  fre¬ 
quencies  and  bandwidths  of  these  noise  types.  Time-domain  convolution  filters 
are  used  on  the  time  series  of  the  four  intensities  and  six  phase  differences. 
Low-pass  filtering  is  especially  effective  on  the  phase  differences.  For  any 
one  of  the  10  time  series  there  may  be  one,  none,  or  several  such  filters 
acting  on  different  frequencies  with  different  bandwidths. 

4.4  Corea  for  Micrometeorological  Data  Problems 

For  those  instruments  that  are  not  functioning  properly  during  a  data  run 
or  session,  we  simply  process  their  signal  anyway.  It  is  then  necessary  to 
ignore  the  reduced  data  from  these  instruments  when  the  averaged  results  are 
Interpreted.  In  session  1,  this  is  the  case  with  the  Lyman-a  hygrometer  at 
station  1  as  well  as  with  one  of  the  optical  wind  and  Cn  meters.  The  wet 
bulb  in  the  psychrometer  at  station  1  is  dry  over  most  of  session  1. 

The  platinum  resistance-wire  thermometers  have  numerous  cold  and  hot 
spikes  that  have  the  same  signature  as  those  observed  in  the  marine  surface 
layer  by  Schmitt  et  al.  (1978).  These  are  particularly  evident  in  session  1. 
In  the  marine  case  the  wires  rapidly  become  contaminated  by  salt  spray;  the 
hydroscopic  nature  of  the  sodium  clorlde  causes  random  evaporative  cooling  and 
condensatlve  heating  of  the  wire  because  the  fluctuating  humidity  continually 
drives  the  adhering  salt-water  solution  away  from  equilibrium.  We  hypothesize 
that  our  signals  are  contaminated  by  some  hydroscopic  particles  adhering  to 
our  wires.  This  effect  has  never  before  been  observed  in  over-land  data.  We 
cannot  test  this  hypothesis  because  the  experiment  session  ended  before  the 
effect  was  observed.  However,  considering  only  session  1  data,  it  is  a 
plausible  hypothesis  because  of  the  high  relative  humidities  (above  85Z)  and 
the  close  similarity  of  the  signal  to  that  observed  for  salt-spray- 


contaminated  data.  We  do  not  know  the  source  of  hydroscopic  aerosol,  but  we 
speculate  that  it  may  be  related  to  fertilization  of  the  soil  or  herbicide 
spraying  or  washing  of  the  wires  in  a  weak  solution  of  baking  soda  during 
manufacture.  Field  tests  were  conducted  during  session  4  under  high  humidity 
conditions,  but  the  effect  did  not  reappear  during  the  test.  In  session  2  the 
electronics  were  tested  as  a  source  of  these  spikes  but  found  not  to  be  at 
fault.  However,  the  hot  and  cold  spikes  were  present  in  data  taken  during 
session  2. 

The  cold-  and  hot-spike  phenomenon  has  serious  consequences  for  our  sta¬ 
tistics.  As  shown  by  Schmitt  et  al.  (1978),  temperature  structure  parameters 
and  heat  fluxes  are  seriously  overestimated  in  marine  surface-layer  data.  The 
effects  are  spectrally  broad  so  that  simple  digital  filtering  will  not  help. 

In  effect,  the  resistance-wire  signal  is  responding  to  a  combination  of  the 
air  temperature  variation  and  the  time-derivative  of  the  humidity.  In  the 
future  we  will  examine  whether  or  not  this  effect  can  be  parameterized.  If 
so,  then  it  may  be  possible  to  use  the  temporal  derivative  of  humidity  from 
the  Lyman-a  hygrometer  to  correct  the  temperature  signal  from  the  resistance 
wire  that  is  within  the  Lyman-a  gap.  Such  a  correction  will  be  a  run-by-run, 
trial-and-error  approach  based  on  the  parameterization  having  variable  parame¬ 
ters.  In  the  mean  time  we  calculate  contaminated  statistics  with  a  cautionary 
note  as  to  their  accuracy. 

There  are  occasional  spikes  in  one  or  more  of  the  100-Hz  channels  that  are 
of  electronic  origin  unrelated  to  the  above-mentioned  signal  contamination. 

We  observe  these  on  computer-generated  microfilm  strlpcharts  and  note  times 
immediately  preceding  and  following  these  spikes.  The  data  processing  then 
skips  all  of  the  given  station's  100-Hz  and  25-Hz  data  within  the  Interval 
defined  by  these  noted  times.  This  is  the  same  procedure  as  was  used  for  the 
mm-wave  processing,  and  it  Imposes  the  same  complications  on  the  processing 
sof  tware . 

During  experiment  sessions  1  and  2,  all  three  axes  of  both  sonic  anemome¬ 
ters  produced  numerous  spikes  that  were  triggered  by  power-line  noise  and 
sonic  electronics  that  were  excessively  sensitive  to  such  effects.  These 
spikes  occurred  at  a  rate  of  roughly  one  per  several  seconds  per  sonic  axis. 
Their  occurrences  and  strengths  are  random  and  sometimes  coincident  in  several 
axes.  Fortunately,  in  almost  all  cases  a  spike  consists  of  only  one  erroneous 


datum.  A  despiking  subroutine  has  been  written  that  recognizes  the  spikes  by 
the  sudden  change  of  signal,  first  in  one  direction  then  in  the  other,  that 
exceeds  a  threshold  set  at  1  standard  deviation  of  the  time  series.  This 
threshold  was  chosen  by  trial  and  error.  If  the  standard  deviation  of  the 
despiked  data  differs  significantly  from  that  of  the  original  data,  then  the 
despiking  is  repeated.  This  procedure  removes  the  spikes  and  has  very  little 
effect  on  the  data.  Of  many  thousands  of  spikes  removed,  only  four  double¬ 
datum  spikes  of  relatively  low  amplitude  remain  in  the  data  of  session  1. 
Although  the  despiking  subroutine  can  remove  double-datum  spikes  by  first 
averaging  the  data  in  temporally  adjacent  pairs,  these  few  remaining  spikes 
are  judged  to  be  too  insignificant  for  this  procedure.  Moreover,  the 
remaining  spikes  can  be  excluded  from  the  data  processing  by  previously 
discussed  methods.  The  problem  of  the  sonic  spikes  is  thus  solved. 

The  existence  of  sonics  spikes  prohibited  nulling  of  the  sonic  electronics 
in  the  field  for  sessions  1  and  2.  Normally  this  null  calibration  produces 
zero  signal  from  the  sonics  when  the  axes  are  covered  so  that  no  wind  is  pre¬ 
sent.  Thus  all  the  sonic  axes  have  offset  wind  readings  amounting  to  as  much 
as  ±0.6  m/s  per  axis.  This  is  not  a  problem  for  the  vertical  axes,  because  we 
subtract  <w>  for  all  of  our  calculated  statistics  involving  w,  and  the  ver¬ 
tical  axes  were  maintained  within  0.1°  of  vertical,  so  horizontal  wind  fluc¬ 
tuations  negligibly  contaminate  the  measured  vertical  fluctuations.  The 
offset  is  serious  for  the  horizontal  axes.  Fortunately  there  is  an  adjacent 
prop-vane  anemometer  with  which  to  compare  wind  readings,  so  the  offset  can  be 
measured. 

Take  the  sonic  A  component  as  an  example.  The  true  sonic  A  component,  A£ , 
is  given  in  terms  of  the  measured  and  despiked  A  component,  A  ,  and  the  offset 

°*  by 

A  ■  (A  -  0.)  K_i  ,  (82) 

t  m  A 

where  K  is  the  wake  correction  factor  given  by  Kaimal  (1980b).  The  inverse  of 
K  appears  on  the  right  side  of  (82)  because  Kaimal  defines  K  as  the  ratio  of 
measured  to  true  sonic  components.  Also,  K  is  a  function  of  the  angle  between 
the  wind  vector  and  the  sonic  arm.  We  would  like  to  calculate  the  A  component 
from  the  prop-vane  anemometer,  Ap,  and  require  that  <Ap>  »  <At>.  However,  we 
do  not  a  priori  know  the  sonic-derived  angle  between  the  wind  vector  and  the 


sonic  arm  because  of  the  offset.  Thus  we  cannot  make  the  wake  correction  to 

A  or  eliminate  data  when  the  wind  is  outside  of  the  ±45°  acceptance  angle, 
in 

Instead,  we  average  (82),  eliminating  data  when  the  vane  angle  is  outside  the 
±45®  acceptance  angle,  and  similarly  average  Ap.  Replacing  <At>  by  <Ap>,  we 
have  from  (82), 

<A  >  =*  <A  K_l>  -  0.  <K"l>  ,  (83) 

pm  A 

wherein  K  is  calculated  using  the  vane  angle  to  infer  the  angle  between  the 
wind  vector  and  sonic  arm.  The  offset  is  obtained  from  (83): 

0*  ■  (<A_  K_l>  "  <A  >)/<K_1>  .  (84) 

Am  p 

To  obtain  0^,  we  thus  calculate  the  three  averaged  quantities  in  (84)  using 
the  vane  angle  for  both  K  and  the  ±45®  acceptance-angle  criterion.  Equation 
(84)  is  not  exact,  but  it  need  only  be  a  good  statistical  approximation.  We 
have  had  excellent  results  using  (84).  That  is,  using  (84)  for  0A  in  (82)  and 
the  analogous  calculation  for  the  B  sonic  component,  we  then  obtain  the  angles 
between  the  wind  vector  and  sonic  arms  from  the  sonic  data.  We  use  this 
sonic-derived  angle  to  determine  both  K  and  the  ±45®  acceptance-angle  cri¬ 
terion.  The  resulting  At  gives  <Afc>  *  <Ap>  to  excellent  approximation.  The 
simpler  method  of  setting  K  ■  1  in  (84)  and  not  eliminating  any  data  based  on 
the  acceptance-angle  criterion  gives  poorer  results. 

5.  OVERVIEW  OF  THE  DATA  PROCESSING  SOFTWARE 
5.1  Introduction 

As  mentioned  in  Sec.  1,  the  data  processing  software  is  designed  to  pro¬ 
duce  answer  files  that  do  not  limit  interpretation  of  the  data,  to  cure  afore¬ 
mentioned  problems  within  the  data,  and  to  produce  graphic  visualization  of 
data  validity  at  several  steps  in  the  processing.  The  first  step  is  to  make  a 
copy  of  each  data  tape  using  a  PDP-11/70  computer.  This  protects  the  original 
data  tape  from  being  destroyed  or  altered  in  subsequent  use.  The  remainder  of 
processing  is  done  on  a  Cyber  750  computer.  This  section  describes  the  soft¬ 
ware  used  on  the  Cyber  750.  Figure  2  is  a  flow  chart  of  the  computer  programs 
and  magnetic  tape  and  disk  storage. 


5.2  FLATCOPY 


The  first  step  in  Cyber  processing  is  program  FLATCOPY  shown  on  Fig.  2. 
FLATCOPY  reads  the  1600-bpi  Flatville  data  tape  (the  copy  made  on  the 
PDP-11/70),  which  has  1024  words  per  block,  converts  the  data  to  60-bit  words, 
and  writes  the  results  to  the  6250-bpi  Cyber-version  data  tape.  It  is  this 
translated  Cyber  version  that  is  efficiently  read  by  the  Cyber  750  in  sub¬ 
sequent  processing.  By  using  6250  bpi  and  blocks  of  4096  words,  a  single 
Cyber-version  magnetic  tape  contains  almost  all  of  the  data  from  a  Flatville 
data  tape;  the  remaining  data,  if  any,  are  discarded. 

FLATCOPY  also  replaces  the  missing  sync  bytes  and  data  frames,  and  locates 
the  missing  data  blocks.  The  times  of  occurrence  of  skipped  data  frames  and 
of  the  approximate  first  and  last  times  bracketing  a  block  loss  are  written  to 
a  disk  file  known  as  the  pedigree  file;  this  file  is  shown  on  Fig.  2.  The 
pedigree-file  information  is  used  in  subsequent  processing. 

FLATCOPY  also  extracts  the  0.39  Hz  data  and  averages  the  Lyman-a  voltage 

counts  as  well  as  platinum  resistance-wire  temperatures  into  2.56-8  bins. 

This  information  is  written  to  a  2.56-s  meteorological  data  base  disk  file 

that  is  used  in  subsequent  processing.  This  2.56-s  data  base  disk  file  is 

shown  on  Fig.  2.  The  following  processing  of  0.39-Hz  data  occurs  in  FLATCOPY. 

The  optical  disdrometer  data  are  converted  to  drop-size  distribution  and  the 

rain  rate  in  volts;  both  are  recorded  every  2.56  s  in  the  disk  file.  The 

vaisaia  signals  are  converted  to  temperature  in  Celsius  and  to  both  relative 

humidity  and  absolute  humidity,  and  recorded  every  2.56  s  in  the  disk  file. 

The  barometer  signal  is  converted  to  millibars  and  likewise  sent  to  the  2.56-s 

disk  file.  Every  2.56  s,  the  psychroraeter  wet-  and  dry-bulb  temperatures  are 

derived  in  Celsius;  absolute  humidity  is  also  obtained,  and  all  three  are 

written  to  disk.  Psychrometer  data  are  rejected  if  a  frame  replacement 

occurred  during  the  psychrometer  commutation  cycle;  otherwise,  spurious 

2 

results  would  be  obtained.  The  optical  cross-wind  and  C  meter  data  are 

n-  2/  3 

obtained  every  2.56  s  in  units  of  meters  per  second  and  m  and  likewise 
written  to  the  disk  file.  The  angle  encoder's  angles  in  degrees  as  well  as 
the  vane  angles  averaged  in  2.56-s  bins  and  prop-vane  wind  speeds  averaged  in 
2.56-s  bins  are  written  to  the  2.56-s  data  base. 

For  the  fourth  and  fifth  experiment  sessions  additional  optical  instrumen¬ 
tation  was  deployed  and  their  outputs  are  converted  to  engineering  units  and 
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written  to  the  2.56-s  data  base.  One  such  instrument  was  an  additional  opti¬ 
cal  cross-wind  and  C  2  meter  at  a  lower  height  directly  beneath  the  station  2 

n 

optical  cross-wind  and  C n2  meter.  The  purpose  is  to  measure  shear  in  the 
optical  cross  wind.  Another  pair  of  optical  sensors,  a  Ha-Ne  laser  scin¬ 
tillometer  with  a  large-aperture,  modulated  LED  scintillometer,  are  used 

2 

together  to  obtain  the  inner  scale  of  turbulence.  The  nominal  Cn  values 
obtained  from  both  scintillometers  using  inertial-range  formulas  for  their 
log-amplitude  variances  are  recorded  in  the  2.56-s  data  base  as  are  the 
resulting  optical  inner-scale  values. 

For  rain,  fog,  snow,  or  blowing  snow  the  2.56-s  data  base  contains  all  of 
the  NOAA  meteorological  information.  An  exception  is  experiment  session  5, 
during  which  the  Lyman-a  hygrometer  and  its  resistance-wire  thermometer  are 
shielded  from  precipitation  in  an  attempt  to  measure  the  turbulence  fluc¬ 
tuations  during  precipitation  events.  For  clear-air  conditions,  data  from  the 
sonic  anemometer,  Lyman-a  hygrometer,  and  resistance-wire  thermometers  must 
also  be  analyzed;  this  is  discussed  in  the  following  sections* 

5.3  METAPE 

Program  METAPE  reads  the  100-Hz  and  25-Hz  microraeteorological  data  from 
the  Cyber-version  data  tape,  effects  cures  to  the  sonic  anemometer  problems, 
converts  the  signals  to  engineering  units,  and  writes  the  results  to  a 
meteorological  working  tape.  This  MET  tape  is  shown  in  Fig.  2.  To  this  end 
it  also  reads  the  2.56-s  data  base  disk  file  and  performs  whole-tape  averages 
of  the  temperature,  humidity,  and  Lyman-a  voltage  counts;  these  quantities  are 
used  in  the  conversion-to-units  for  the  sonic  anemometers  and  Lyman-a  hygrome¬ 
ters.  The  output  to  the  MET  working  tape  consists  of  all  of  the  resistance- 
wire  temperatures  and  Lyman-a  humidities  written  at  100  Hz,  as  well  as  the 
solar  flux,  sonic  wind  components,  sonic  horizontal  wind  speeds,  prop-vane 
wind  speeds  and  angles  written  at  25  Hz;  all  these  quantities  are  written  in 
engineering  units.  In  addition,  two  flags  are  written  at  the  25-Hz  rate  to 
the  MET  working  tape;  these  flags  designate  whether  or  not  the  wind  angle  at 
the  sonic  anemometer  lies  within  the  ±45*  acceptance  range  of  the  sonic-axis 
configuration.  These  flags  are  used  in  later  processing  to  reject  sonic  data 
for  wind  angles  beyond  ±45°. 

The  handling  of  the  sonic  data  is  particularly  complex.  Program  METAPE 
first  reads  the  sonic  components  and  prop-vane  signals  from  the  Cyber-version 
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data  tape.  If  necessary,  the  sonic  signals  are  despiked  in  intervals  of  40.96 
s  of  data  at  a  time  and  then  converted  to  units  of  meters  per  second.  The 
prop-vane  signals  are  likewise  converted  to  wind  speed  in  meters  per  second 
and  to  angle  in  degrees.  These  results  are  written  to  a  temporary  disk  file. 
Meanwhile  the  averages  in  (84)  are  computed  so  that  the  sonic  offsets  can  he 
determined  as  discussed  in  Sec.  4.4.  When  the  entire  Cyber-version  data  tape 
has  been  thusly  processed,  the  offsets  for  the  horizontal  sonic  wind  com¬ 
ponents  are  obtained  from  (84).  The  exhausted  Cyber  tape  is  now  rewound  and 
read  from  the  beginning  to  provide  the  resistance-wire  temperatures,  Lyman-a 
humidities,  and  solar  flux.  Simultaneously,  the  temporary  disk  file  is  re¬ 
wound  and  read  to  provide  the  prop-vane  speeds  and  angles  already  converted  to 
units  and  the  despiked  sonic  wind  components.  The  sonic  offsets  are  then 
subtracted  from  their  respective  components.  Note,  however,  that  determining 
the  sonic  offsets  by  comparison  with  the  prop-vane  anemometer  is  an  option 
that  is  not  used  in  processing  every  data  tape;  alternatively,  the  offsets 
from  the  field  log  can  be  used.  At  the  25-Hz  rate  the  angle  between  the  hori¬ 
zontal  wind  vector  and  each  sonic  axis  is  calculated  and  used  to  correct  the 
sonic  components  for  their  wake  averaging  and  to  set  the  flags  designating 
whether  or  not  the  wind  angle  is  within  the  ±45°  acceptance  range.  The  wake¬ 
averaging  correction  is  given  by  Kaimal  (1980b).  Although  the  wake  of  one 
sonic  transducer  crosses  the  second  sonic  axis  when  the  wind  angle  is  beyond 
i45# ,  wind  tunnel  tests  of  a  related  sonic  anemometer  by  Hanafusa  et  al. 

(1980)  show  no  effect  on  the  wind  measured  along  this  second  sonic  axis. 

Using  the  wake  correction  gives  us  confidence  in  the  ±45°  acceptance  angle. 

The  horizontal  sonic  components  are  converted  to  orthogonal  components  as 
described  in  Sec.  3.2.  The  horizontal  wind  speed  is  calculated  from  the 
latter  sonic  components.  Finally  the  data  are  written  to  the  MET  working 
tape. 

5.4  Stripchart  Data  Visualization 

We  now  discuss  the  programs  in  Fig.  2  labeled  GRAPH  2.56-s  DATA  BASE, 
STRIPCHART  MIL-WAVE,  and  STRIPCHART  MICRO  MET.  These  programs  check  the  data 
validity  and  yield  the  time  Intervals  that  we  must  delete  from  further  pro¬ 
cessing  because  of  erroneous  data.  All  graphs  are  produced  on  microfilm  so 
that  the  enormous  amount  of  graphical  Information  can  be  conveniently  stored. 


In  Sec.  5.2  we  discussed  all  of  the  quantities  that  program  FLATCOPY 
archived  in  the  2.56-s  data  base  disk  file.  The  program  labeled  GRAPH  2.56-s 
DATA  BASE  in  Fig.  2  makes  a  graph  of  each  of  these  quantities  versus  time 
during  the  data  run,  and  one  point  is  plotted  for  each  2.56  s.  In  addition  it 
graphs  the  averaged  rain-drop  size  distribution.  It  also  graphs  all  of  the 
temperatures  on  a  single  graph;  these  are  the  four  2.56-s-averaged  resistance- 
wire  temperatures  and  the  two  dry-  and  wet-bulb  psychrometer  temperatures  as 

well  as  the  two  VSis&lS  temperatures.  These  graphs  show  us  the  trends  in  tem- 
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perature,  humidity,  optical  Cn  and  cross  wind,  rain  rate,  solar  flux,  and 
even  barometric  pressure.  At  this  point  the  malfunction  of  instruments  be¬ 
comes  evident.  For  instance,  the  fact  that  the  wet  bulb  in  the  psychrometer 
at  station  2  is  actually  dry  is  evident  in  session  1,  and  program  METAPE  is 
then  commanded  to  use  the  humidity  from  the  psychrometer  at  station  1  to  con¬ 
vert  the  Lyman-a  signals  from  station  2  to  physical  units.  The  failure  of  one 

o 

of  the  optical  cross-wind  and  Cn  meters  is  evident  in  session  1.  The  trends 
in  temperature  reveal  whether  such  trends  will  dominate  temperature  variance 
calculated  over  the  entire  data  tape. 

The  program  labeled  STRIPCHART  MIL-WAVE  in  Fig.  2  produces  a  stripchart  of 
the  six  phase  differences  and  four  intensities  on  one  microfilm  stripchart; 

10  8  of  data  is  displayed  in  each  frame.  A  related  program  can  produce  a 
stripchart  of  the  raw  mm-wave  voltage  squares  if  need  be.  These  strlpcharts 
reveal  the  presence  of  phase-lock  losses,  spikes,  step  changes,  or  other 
problems  in  the  mm-wave  data.  The  time  Intervals  to  be  skipped  in  subsequent 
processing  are  obtained  from  these  strlpcharts.  The  program  also  produces 
printed  summaries  of  the  averages,  rms,  and  extrema  of  the  10  quantities  in 
100  s  intervals.  In  addition,  averages  and  rms  of  the  10  quantities  are  also 
produced  for  20-s  intervals  and  are  graphed  with  the  entire  data  run  displayed 
on  single  graphs  for  each  of  the  quantities.  These  latter  graphs  display  the 
trends  in  the  mm-wave  data. 

The  program  labeled  STRIPCHART  MICRO  MET  in  Fig.  2  is  actually  two  pro¬ 
grams.  One  produces  a  single  microfilm  stripchart  for  data  from  the  four 
resistance-wire  thermometers  and  two  Lyman-a  hygrometers  at  the  rate  of  10  s 
per  frame.  The  second  produces  strlpcharts  of  the  u,  v,  w  components  of  the 
sonic  anemometers  and  overplots  on  these  the  corresponding  u  and  v  components 
calculated  from  the  prop-vane  anemometers;  this  gives  10  curves  per  stripchart 


at  a  rate  of  20  s  per  frame.  These  stripcharts  yield  the  time  intervals  to  be 
skipped  in  later  processing  in  case  of  erroneous  data.  As  with  the  mm-wave 
stripcharts,  summaries  are  printed  of  100-8  averages,  rms,  and  extrema,  and 
20-s  averages  and  rms  are  plotted  for  each  quantity  to  reveal  trends.  The 
extrema  are  used  to  set  the  minimum  and  maximum  bins  for  the  probability 
distribution  function  calculations. 

5.5  PED  2  PAR  and  MODPARM 

Several  parameters  are  needed  for  the  processing  of  the  data.  These 
include  times  bracketing  erroneous  data  that  is  to  be  skipped  during  pro¬ 
cessing,  the  location  of  spikes,  the  integration  period  to  be  used,  etc.  In 
addition,  missing  block  locations  obtained  from  the  pedigree  file  must  be 
provided  to  the  processing  programs.  A  compact  file  called  the  parameters 
file  has  been  defined  to  retain  these  parameters  for  each  tape  to  be  pro¬ 
cessed.  Two  programs  create,  modify,  and/or  display  the  contents  of  the 
parameter  file.  These  are  labeled  PED  2  PAR  and  MODPARM  in  Fig.  2. 

The  PED  2  PAR  program  allows  two  modes  of  operation.  If  a  non-empty  ped¬ 
igree  file  is  provided,  then  PED  2  PAR  reads  the  file  and  extracts  from  it 
the  timing  data  identifying  missing  data  blocks.  PED  2  PAR  searches  the  para¬ 
meters  file  for  the  times  to  be  skipped  during  processing  and  either  replaces 
the  missing  block  times  if  there  is  an  existing  entry  or  makes  a  new  entry  if 
there  is  not.  In  the  second  mode  of  operation,  PED  2  PAR  displays  the  missing- 
block  data  for  all  data  tapes  in  the  parameters  file  in  both  absolute  time 
(Greenwich  Mean  Time)  and  time  relative  to  the  beginning  of  the  tapes. 

The  MODPARM  program  allows  Interactive  access  to  the  parameters  file. 

The  program  first  attempts  to  read  the  file  and  verify  that  the  version  of  the 
file  is  at  least  the  minimum  needed  for  the  current  version  of  the  MODPARM 
program.  MODPARM  then  presents  a  menu  from  which  to  choose  the  data  tape  of 
interest  and  those  parameters  from  that  tape  that  the  user  wishes  to  Inspect 
or  modify.  The  program  then  displays  the  current  values  of  the  parameters, 
and  the  user  is  given  an  opportunity  to  modify  them,  thereby  altering  the 
parameters  file.  An  example  of  modifying  the  parameters  file  is  to  enter  the 
beginning  and  ending  times  of  Intervals  to  be  skipped  in  subsequent  pro¬ 
cessing;  Identification  of  such  time  intervals  using  our  stripcharts  was 
discussed  in  the  previous  section. 


As  shown  in  Fig.  2,  the  parameters  file  is  read  by  the  MILDRV  program  to 
obtain  the  parameters  needed  for  performing  its  processing  steps.  The  output 
file  MILSUMS  contains  a  copy  of  the  values  of  the  parameters  file  so  that  a 
record  is  available  to  confirm  the  validity  of  parameters  used  in  processing. 

5.6  MILDRV  and  METDRV 

The  quantities  calculated  as  averages  from  the  ram-wave  and  mlcrometeorolo- 
gical  data  were  discussed  in  Secs.  2.2  to  2.6  and  in  Sec.  3.3.  These  quan¬ 
tities  are  calculated  as  running  sums  (rather  than  averages)  in  the  programs 
MILDRV  and  METDRV  for  the  mm-wave  and  mlcrometeorological  data  respectively. 
The  advantages  of  outputting  running  suras  over  time  intervals  much  shorter 
than  a  data  run  will  become  clear  in  the  next  two  sections.  As  shown  in  Fig. 
2,  MILDRV  reads  the  data  from  the  Cyber-version  tape  and  outputs  the  running 
sums  to  a  disk  file  called  MILSUMS,  and  METDRV  reads  the  MET  tape  and  outputs 
to  METSUMS  (actually  METSUMS  consists  of  several  distinct  disk  files).  We 
output  the  running  sums  for  specified  averaging  periods;  at  present  we  use 
25.6-s  periods.  The  one  exception  is  the  structure  parameters,  which  we  out¬ 
put  as  averages.  Also,  the  extrema  found  in  each  25.6-s  period  are  recorded 
in  the  disk  files.  The  PDFs  and  JPDFs  are  also  calculated  by  MILDRV  and 
METDRV  and  written  to  the  disk  files  as  unnorraalized  histograms.  The  number 
of  points  summed  is  also  written  to  these  disk  files;  this  number  is  variable 
because  erroneous  data  are  skipped,  and  it  is  different  for  mm-wave  data  ver¬ 
sus  mlcrometeorological  data.  This  number  can  differ  between  mlcrometeorolo- 
glcal  stations  because  only  one  station  may  have  erroneous  data  that  is 
skipped,  and  different  statistical  quantities  from  a  given  station  may  have 
different  numbers  of  points  summed  because  sonic  data  may  be  dropped  due  to 
wind  angles  outside  the  ±45*  acceptance  range  and  because  only  some  statistics 
depend  on  sonic  data.  The  data  reduction  in  MILDRV  and  METDRV  obviously 
depends  on  the  strlpcharts  in  an  essential  way  because  the  strlpcharts  yield 
the  time  Intervals  containing  erroneous  data  that  is  to  be  skipped. 

It  is  MILDRV  that  cures  problems  in  the  nm-wave  data,  as  discussed  in 
Sec.  4.3.  MILDRV  calculates  the  polynomial-fit  correction  to  all  28  mm-wave 
channels  and  synthesizes  1^  for  experiment  session  1.  MILDRV  also  performs 
the  digital  filtering  needed  to  remove  erroneous  trends  in  the  mm-wave  data, 
and  to  remove  high-frequency  noise  and  noise  due  to  vibration  of  the  I-beam 
holding  the  mm-wave  receivers. 
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MILDRV  has  a  second  node  of  processing  in  support  of  the  ton-wave  fast 
Fourier  transform  program,  MILFFT.  In  this  mode  MILDRV  reads  the  Cyber¬ 
version  data  tape,  performs  the  polynomial-fit  correction  to  the  28  ram-wave 
channels,  synthesizes  1^  if  session  1  data  are  used,  and  calculates  the  phase 
differences.  However,  in  this  mode  MILDRV  does  not  perform  any  digital 
filtering.  In  so  doing,  MILDRV  uses  the  parameters  file  to  locate  erroneous 
or  lost  data  and  thereby  Identify  all  sequential  5.12-s  intervals  of  tem¬ 
porally  continuous  data.  The  four  intensities  and  six  phase  differences  are 
then  output  to  a  temporary  disk  file  spanning  the  entire  data  tape  in  blocks 
of  5.12-s  duration.  This  temporary  disk  file  contains  gaps  caused  by  erro¬ 
neous  or  lost  data.  A  new  5.12-s  block  is  attempted  beginning  with  the  first 
good  datum  following  erroneous  or  lost  data.  Thus  the  temporal  extent  of  a 
gap  is  from  the  end  of  one  5.12-s  block  to  the  next  good  datum  that  is 
followed  by  5.12  s  of  temporally-contlnuous  and  error-free  data.  The  tem¬ 
porary  disk  file  is  used  by  MILFFT. 


5.7  Graphing  the  25.6-s  Statistics 

The  program  labeled  GRAPH  25.6-s  STATISTICS  in  Fig.  2  is  actually  two 
programs;  one  graphs  the  results  in  MILSUMS  and  the  other  those  in  METSUMS. 
This  program  takes  the  running  sums  in  the  disk  files,  produces  averages  by 
dividing  by  the  appropriate  number  of  points  summed,  converts  such  quantities 
as  mean  squares  to  more  meaningful  quantities  like  central  variances,  and 
finally  graphs  these  quantities  versus  time  with  one  point  every  25.6  s  per 
statistic. 

At  this  point  these  graphs  give  us  a  look  at  the  trends  in  all  of  our  sta¬ 
tistics,  only  a  few  of  which  are  redundant  with  the  trends  visualized  In  the 
previous  graphing  of  the  2.56-8  data  file  and  the  stripcharting.  The  graphs 
of  the  extrema  are  particularly  sensitive  to  erroneous  data  that  might  have 
eluded  us  or  were  wrongly  included  because  the  time  intervals  skipped  were 
Incorrect  when  the  data  were  input  to  MILDRV  or  METDRV.  If  an  error  Is  found, 
the  runs  of  MILDRV  and/or  METDRV  can  be  repeated. 

MILSUMS  and  WTSUMS  are  repairable  files.  If  an  error  occurs  in  a  25.6-s 
bln,  then  MILDRV  or  FCTDRV  can  recalculate  just  this  one  bln  without  repro¬ 
cessing  the  entire  MET  or  Cyber-version  tapes.  It  is  then  possible  to  rewrite 
the  MILSU^  or  METSUMS  files  including  the  corrected  25.6-s  bin. 


We  check  these  graphs  for  statistical  stationarity  and  specify  times  tj  to 
t2  over  which  the  statistics  are  sufficiently  stationary.  A  given  data  tape 
may  have  more  than  one  period  of  stationarity  without  being  stationary  as  a 
whole.  An  extreme  example  would  be  a  data  tape  that  consists  of  clear-air 
data  for  part  of  the  run  but  rain  data  over  the  remainder  of  the  run.  In  this 
case  we  can  use  the  graphs  to  specify  times  tj  to  t2  for  later  accumulating 
the  clear-air  data  and  times  t^'  to  t2*  for  accumulating  the  rain  data.  By 
studying  these  graphs,  we  can  include  or  exclude  data  at  will  in  subsequent 
processing. 

5.8  MILRDR  and  METRDR 

Given  an  interval  t^  to  t2  of  stationary  statistics  as  described  In  the 
previous  section,  the  programs  labeled  MILRDR  and  METRDR  on  Fig.  2  read  the 
disk  files  MILSUMS  and  METSUMS  and  sum  the  25.6-s  partial  sums  from  tj  to 
t2  while  also  summing  the  numbers  of  points  in  each  25.6-s  partial  sura.  These 
total  sums  from  t^  to  t2  are  converted  to  averages  by  dividing  by  the  total 
number  of  points  summed,  and  quantities  like  mean  squares  are  converted  to 
more  meaningful  quantities  like  central  variance.  The  PDFs  and  JPDFs  stored 
for  each  25.6  s  are  summed  and  then  normalized.  The  extrema  of  the  25.6-s 
data  are  found.  The  resulting  averages  (which  Include  such  quantities  as 
variances  and  cross-correlations),  the  PDFs  and  JPDFs,  and  the  extrema  are 
identical  to  those  that  we  could  calculate  from  tj  to  t2  if  we  skipped  the 
intermediate  step  of  storing  25.6-s  running  sums  and  histograms.  The  excep¬ 
tions  to  this  identity  are  the  structure  parameters  that  are  stored  in  the 
disk  file  METSUMS  as  averages.  METRDR  then  averages  these  25.6-s-averaged 
structure  parameters  from  tj  to  t2.  This  is  actually  an  advantage  for  the 
time-lagged  structure  parameters  because  wind  speed  variability  is  less  on  a 
25.6-s  time  Interval  than  from  tj  to  t2,  and  the  relationship  between  time- 
lagged  and  spatial  structure  parameters  relies  on  wind  speed  stationarity. 

Final  normalizations  of  ram-wave  data  are  performed  by  MILRDR.  For 
Instance,  the  MCF  and  fourth  moments  of  the  field  are  converted  to  their  nor¬ 
malized  values,  as  discussed  in  Secs.  2.5  and  2.6. 

Finally,  MILRDR  and  METRDR  write  the  accumulated  statistics  to  the  disk 
file  labeled  ANSWER  FILES  in  Fig.  2.  The  answer  files  may  contain  only  one 
set  of  statistics  for  a  given  data  tape  if  only  one  set  of  times  t^  and  t2  is 


specified,  but  it  may  contain  additional  sets  if  additional  times  tj'  and  t 2' 
are  specified.  We  have  already  noted  the  advantages  of  data  validation  and 
data  selection  in  creating  the  files  METSUMS  and  MILStJMS  of  25.6-s  running 
sums.  In  addition,  this  partial  sum  procedure  allows  additional  processing  of 
differing  sets  of  times  t^  to  t2  to  be  performed  in  the  distant  future  without 
recalculation  at  the  high  data  rates  of  100  Hz  and  25  Hz. 

5.9  MILFFT  and  METFFT 

Sections  2.9  and  3.9  describe  what  spectra  and  cross-spectra  are  calcu¬ 
lated  and  how  this  is  accomplished  for  mm-wave  and  micrometeorological  data, 
respectively.  The  programs  MILFFT  and  METFFT  in  Fig.  2  accomplish  these 
transforms  and  write  the  averaged  results  to  the  answer  files.  In  addition, 
mm-wave  spectra  that  are  not  averaged  in  frequency  bins  are  saved  and  examined 
so  that  noise  in  the  mm-wave  intensities  and  phase  differences  can  be  iden¬ 
tified  for  subsequent  filtering  done  by  MILDRV  in  producing  the  MILSUMS  file. 

Section  5.6  describes  the  temporary  disk  files  created  by  the  second  mode 
of  MILDRV  processing.  These  files  are  then  read  by  MILFFT  to  produce  the 
Fourier  transforms.  METFFT  creates  its  own  analogous  disk  files  for  Fourier 
transformation  without  a  second  mode  of  processing  that  exists  within  METDRV. 
METFFT  does  not  use  the  parameters  file  and  therefore  cannot  transform  data 
tapes  that  have  block  drops;  data  tapes  preceding  tape  12  of  session  2  are 
therefore  excluded.  However,  times  bracketing  bad  data  observed  on  the  micro- 
meteorological  stripcharts  are  entered  so  that  bad  data  do  not  enter  the  disk 
files  to  be  transformed. 

In  this  mode  METFFT  reads  the  MET  working  tape  and  outputs  5.12-s  blocks 
of  temporal ly-continuous  error-free  data  to  four  temporary  disk  files.  The 
choice  of  these  5.12-s  blocks  and  the  temporal  extent  of  gaps  between  them,  if 
any,  is  the  same  as  for  MILDRV  in  Sec.  5.6  (except  block  drops  are  not  con¬ 
sidered).  There  are  two  disk  files  for  each  micrometeorological  station.  One 
file  holds  the  100-Hz  data,  namely  the  temperatures  from  the  two  resistance- 
wire  thermometers  and  the  Lyman-a  hygrometer.  The  other  file  holds  the  25-Hz 
data,  namely  from  the  three  axes  of  the  sonic  anemometer  and  the  three  100-Hz 
quantities  averaged  by  fours  to  produce  25-Hz  data.  The  temporal  integrity  of 
the  100-Hz  data  derives  from  any  errors  within  the  100-Hz  data,  whereas  that 
of  the  25-Hz  data  derives  from  errors  in  both  the  100-Hz  data  and  the  sonic 


anemometers  and  from  the  wind  angle  being  outside  the  cone  of  acceptance  for 
the  sonic  anemometer. 

Once  the  mm-wave  or  micrometeorological  disk  files  are  created,  M1LFFT  and 
METFFT  use  the  same  versatile  Fourier  transform  program  to  produce  spectra. 

The  data  are  read  from  these  files  in  sets  of  40.96  s  of  continuous  data,  that 
is,  eight  5.12-s  blocks  that  have  no  gaps  beween  them.  The  FFT  is  performed 
on  a  pair  of  quantities  simultaneously  and  written  to  another  temporary  disk 
file,  then  another  pair  is  read  from  the  first  file,  and  so  on.  The  user  can 
select  a  given  number  of  available  40,96-s  intervals,  and  each  is  Fourier 
transformed  and  written  to  disk.  These  transforms  are  read  from  the  disk. 

The  desired  power  spectra  and  cross-spectra  are  formed  for  each  40.96-s  inter¬ 
val  and  are  then  averaged  together  to  increase  statistical  reliability.  These 
spectra  still  retain  their  full  dependence  on  frequency  and  can  be  graphed  to 
diagnose  instrument  problems  and/or  saved  in  a  permanent  disk  file.  Such 
graphs  are  especially  important  for  further  mm-wave  processing  because  the 
frequencies  and  strengths  of  the  various  types  of  noise  in  the  mm-wave  data 
must  be  identified  for  subsequent  digital  filtering  before  the  MILSUMS  file  is 
written.  Finally,  the  spectra  are  averaged  into  logarithmically  spaced  fre¬ 
quency  bins  and  graphed  as  well  as  written  to  the  answer  files. 

To  produce  the  low-frequency  Fourier  transform,  the  data  in  the  temporary 
disk  file  described  in  Sec.  5.6  are  low-pass  filtered  with  a  single-pole 
RC-type  filter  and  then  averaged  into  bins  of  about  0.6  s  such  that  there  are 
4096  temporal  bins  per  data  tape  or  per  given  section  of  a  data  tape.  The 
gaps  beween  the  5.12-s  blocks  are  filled  with  interpolated  values.  These  data 
are  then  Fourier  transformed  in  pairs  of  quantities,  as  for  the  high-frequency 
FFT,  and  written  to  a  temporary  disk  file.  When  all  pairs  have  been  trans¬ 
formed,  the  desired  power  spectra  and  cross-spectra  are  formed  and  corrected 
for  the  effects  of  low-pass-filter  and  bin  averaging.  These  spectra  are  then 
averaged  into  logarithmically  spaced  frequency  bins,  graphed,  and  written  to 
the  answer  files. 


5.10  Using  the  Answer  Files 

The  set  of  all  answer  files  from  all  data  tapes  provides  a  very  versatile 
data  base  for  interpreting  the  data  and  preparing  publications.  This  flexibi¬ 
lity  is  the  desired  end  product  alluded  to  in  Sec.  1.  Computer  programs  can 


be  written  and  modified  to  present  the  reduced  data  in  any  graphical  or  tabu¬ 
lar  manner  desired.  For  instance,  programs  can  be  written  to  plot  all  the 
MCFs  versus  antenna  spacing  for  the  clear-air  data,  with  one  curve  for  each 
desired  refractive-structure  parameter  value,  then  alternatively  plotted  ver¬ 
sus  antenna  spacing  scaled  by  a  theoretical  wave  coherence  length. 
Alternatively,  these  clear-air  MCFs  could  be  plotted  for  various  values  of 
surface-layer  stability  rather  than  refractive-structure-parameter  values. 

The  MCFs  from  turbulence,  rain,  fog,  and  snow  data  could  all  be  plotted  on  a 
single  graph  for  comparison  by  accessing  the  answer  files  for  these  diverse 
meteorological  conditions.  One  could  tabulate  all  of  the  temperature  and 
humidity  structure  parameters  and  cross-structure  parameters  from  the  various 
different  methods  of  calculating  them  (spatial,  time-lagged,  optical  scin¬ 
tillation)  to  search  for  data  consistency  or  errors.  Clearly  these  examples 
are  Insignificant  compared  with  the  myriad  relationships  that  could  be 
investigated  using  the  answer  files. 

5.11  Archiving  of  Disk  Files 

The  Cyber-version  data  tape  and  JET  tape  are,  of  course,  magnetic  tapes 
and  can  therefore  be  stored  for  further  processing  as  long  as  desired  or  prac¬ 
tical.  The  various  disk  files  in  Fig.  2  cannot  reside  on  disk  indefinitely 
because  of  disk-storage  costs.  Consequently,  the  2.56-s  data  base  file  as 
well  as  the  MILSUMS  and  METSUMS  files  are  dumped  to  magnetic  tape  and  deleted 
from  disk  after  the  answer  files  are  made.  Thus  these  files  can  be  returned 
to  disk  if  further  data  processing  is  desired  in  the  future.  Likewise,  the 
answer  files  are  dumped  to  magnetic  tape  and  deleted  from  disk  when  desired. 

It  is  expected  that  the  magnetic  tapes  for  the  answer  files  will  be  retained 
indefinitely. 

6.  REDUCED  MILLIMETER-WAVE  AND 
M1CROMETEOROLOG1CAL  DATA 

6.1  Introduction 

As  an  example  of  our  data  processing  results,  we  present  reduced  data  from 
the  answer  files  of  data  tapes  20,  22  to  25,  and  27  of  experiment  session  1. 

It  should  be  remembered  that  the  Intensity  at  antenna  4  is  synthesised,  as 


discussed  in  Sec.  A. 3.  The  probability  distribution  functions  are  normalized 
to  unit  area  under  a  curve  connecting  the  histogram  peaks.  That  is,  the  sum 
of  the  probability  over  all  bins  is  made  equal  to  the  number  of  bins  divided 
by  the  range  of  the  independent  variable. 

6.2  Tabular  Results 

Tables  4  and  5  are  mm-wave  statistics  averaged  over  the  whole  of  session  1 
for  tapes  24  and  20,  respectively.  In  interpreting  these  results  one  should 
keep  in  mind  that  the  root-mean-square  (rms)  intensity  decreased  by  30%  in  an 
almost  linear  fashion  over  the  duration  of  tape  20  data,  with  a  much  less  dra¬ 
matic  change  in  rms  phase  differences.  On  the  other  hand,  tape  24  had  nearly 
constant  rms  intensities  and  rms  phase  differences. 

As  expected,  the  intensity  variance  is  much  less  than  unity.  Ideally  all 
four  antennas  should  give  equal  variances.  Agreement  is  good,  but  the  synthe¬ 
sized  intensity  at  antenna  4  gives  the  least  agreement  for  tape  24.  The  mini¬ 
mum  and  maximum  Intensities  are  normalized  by  dividing  by  the  mean  intensity. 
Table  4  shows  that  over  the  whole  of  tape  24  the  intensities  swung  from  about 
0.5  to  greater  than  1.5  of  the  mean  intensity. 

The  phase  variances  in  these  tables  differ  from  one  antenna  pair  to  the 
next,  as  they  should  because  of  the  differing  antenna  separations.  The  mini¬ 
mum  and  maximum  phase  differences  likewise  show  corresponding  variability.  As 
seen  in  Table  4,  the  largest  observed  phase  differences  exceeded  3  radians. 

The  mixed  intensity-difference  and  phase-difference  correlation  given  in  Eq. 
(7)  is  very  much  smaller  than  either  the  intensity  variance  or  phase- 
difference  variance.  Thus  amplitude  and  phase  are  relatively  uncorrelated. 

The  phase-type  figures  of  merit  are  much  smaller  than  any  of  the  phase- 
difference  variances;  this  verifies  the  consistency  of  the  phase-difference 
measurements. 

The  p-type  figures  of  merit  in  Tables  4  and  5  suggest  that  the  MCF  should 
not  be  trusted  to  better  than  about  5%  accuracy.  In  Tables  4  and  5  the  p  type 
and  I  type  are  close  to  unity,  as  desired,  even  for  antenna  pairs  3,  5,  and  6, 
which  use  the  synthesized  intensity  at  antenna  4.  The  variances  of  both  p- 
type  and  E-type  figures  of  merit  show  that  there  is  very  little  temporal 
variation  about  the  mean  values  of  these  figures  of  merit. 
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Table  4. — Millimeter-wave  statistics  from  the  answer  files 
obtained  by  averaging  all  of  tape  24  of  session  1 ,  July  1983 


Antenna 

number 

Intensity 

variance 

Intensity  Statistics 

Normalized 

minimum 

intensity 

Normalized 

maximum 

intensity 

1 

0.0210 

0.540 

1.70 

2 

0.0219 

0.523 

1.74 

3 

0.0216 

0.423 

1.75 

4 

0.0287 

0.433 

1.74 

Antenna 

pair 

Phase- 

difference 

variance 

Phase  Difference 

Minimum 

phase 

difference 

Statistics 

Maximum 

phase 

difference 

Eq.  (7) 

1 

0.081 

-1.32 

1.08 

-7.12  (-6) 

2 

0.315 

-2.46 

2.37 

-8.54  (-6) 

3 

0.703 

-3.25 

3.09 

-1.50  (-5) 

4 

0.193 

-1.67 

1.75 

-7.90  (-6) 

5 

0.563 

-3.24 

2.87 

-9.86  (-6) 

6 

0.408 

-2.47 

2.29 

-4.87  (-6) 

Antenna 

pair 

P  type 

Figures  of  Merit 

p-type 

variance 

Z  type 

E-type 

variance 

1 

1.07 

2.85  (-5) 

1.10 

4.66  (-4) 

2 

1.05 

4.48  (-4) 

1.03 

7.41  (-4) 

3 

1.02 

5.85  (-5) 

0.934 

3.13  (-3) 

4 

1.03 

2.57  (-4) 

1.01 

2.06  (-3) 

5 

1.00 

3.61  (-5) 

0.944 

4.18  (-3) 

6 

0.976 

8.43  (-5) 

0.869 

1.75  (-3) 

Antenna 

Phase  type 

numbers 

(a  variance) 

14244 

8.16  (-4) 

14345 

2.51  (-3) 

24346 

1.36  (-3) 

44546 

9.24  (-4) 

1434446 

2.49  (-3) 

Note:  Numbers  in  parentheses  indicate  powers  of  10. 


Table  5. — Millimeter— wave  statistics  from  the  answer  files 
obtained  by  averaging  all  of  cape  20  of  session  l,  July  1983 


Intensity  Statistics 


Antenna 

number 

Intensity 

variance 

Normalized 

minimum 

intensity 

Normalized 

maximum 

intensity 

1 

0.0059 

0.694 

1.34 

2 

0.0046 

0.716 

1.26 

3 

0.0066 

0.638 

1.43 

4 

0.0066 

0.635 

1.34 

Phase 

Difference  Statistics 

Antenna 

Phase- 

Minimum 

Maximum 

Eq.  (7) 

pair 

difference 

phase 

phase 

variance 

difference 

difference 

1 

0.024 

-0.58 

0.68 

-2.06  (-6) 

2 

0.104 

-1.32 

1.19 

-3.59  (-6) 

3 

0.193 

-1.50 

1.74 

-1.50  (-6) 

4 

0.063 

-1.13 

1.02 

-2.18  (-6) 

5 

0.163 

-1.54 

1.41 

-7.09  (-7) 

6 

0.131 

-1.41 

1.38 

-1.45  (-6) 

Figures  of  Merit 


Antenna  p  type  p-type  E  type  E-type 

pair _ variance _ variance 


1 

0.979 

1.27 

( 

-4) 

0.89 

1.28 

(-3) 

2 

1.00 

7.10 

( 

-5) 

0.88 

1.16 

(-3) 

3 

0.964 

1.10 

( 

-4) 

0.81 

1.53 

(-3) 

4 

1.04 

1.71 

( 

-4) 

1.06 

3.02 

(-4) 

5 

1.01 

1.96 

( 

-4) 

1.02 

1.75 

(-3) 

6 

1.02 

4.34 

( 

-4) 

0.99 

6.05 

(-3) 

Antenna  Phase  type 

numbers _ (a  variance) 


1&2&4 

3.00 

(-4) 

14345 

5.83 

(-4) 

24346 

4.99 

(-4) 

44546 

8.57 

(-4) 

1434446 

4.93 

(-4) 

Note: 

Numbers  in 

parentheses  indicate  powers  of  10. 

Table  6  is  a  summary  of  the  micrometeorological  data  from  tape  24,  which 
were  taken  at  about  11:30  a.m,  July  1983.  These  averages  were  obtained  from 
480  s  into  the  run  to  its  end  at  2400  s.  The  values  of  the  structure  parame¬ 
ters  imply  that  the  humidity  structure  parameter  overwhelmingly  dominates  the 
calculated  ram-wave  refractive-index  structure  parameter.  The  values  of  tem¬ 
perature  structure  parameters  from  optical  and  fine-wire  instruments  are  in 
excellent  agreement.  The  small  but  nonzero  value  of  the  mean  of  the  cross¬ 
stream  wind  component  is  a  measure  of  the  inaccuracy  of  determining  the  true 
streamwise  direction. 

6.3  Graphical  Results 

Figure  3  shows  the  intensity  probability  distribution  function  (PDF)  for 
antenna  1  for  session  1,  tape  24.  Gaussian  and  lognormal  PDFs  having  the  same 
mean  and  variance  are  also  plotted.  There  is  little  difference  between  the 
Gaussian  and  lognormal  PDFs  for  such  small  values  of  variance.  However,  the 
data  are  clearly  skewed  and  more  nearly  follow  the  lognormal  PDF,  as  is 
expected  theoretically. 
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Figure  3. — Probability  distribution 
function  of  intensity  for  antenna  1 , 
session  1,  tape  24,  July  1983. 
Gaussian  and  lognormal  distributions 
are  also  shown.  Normalized  Intensity 
is  scaled  to  unit  mean  value.  The 
distribution  is  normalized  to  unit 
area  under  the  histogram. 
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Table  6. — Summary  of  microraeteorological  data  from  tape  24, 
taken  at  11:30  a.m.,  July  1983 


l« 

I. 

k 

i 

t 

L 

I 

S> 


N 

f 


Average  Values 

Humidity 

19  g/m3 

Temperature 

32°C 

Wind  speed 

5.3  m/s 

Wind  angle 

10° 

Pressure 

993  mb 

Solar  flux 

94%  of  full  sun 

Wind  stress 

-0.14  (m/s)2 

Humidity  flux 

0.1  (g/m2)/s 

Temperature  flux 

0.03°C  ra/s 

Stability 

-0.03 

Square  Roots  of  Variances 

Humidity 

0.72  g/m3 

Temperature 

0.35°C 

Wind  speed 

1.2  m/s 

Wind  angle 

11* 

Streamwise  wind  component 

1.1  m/s 

Cross-stream  wind  component 

1.0  m/s 

Vertical  wind  component 

0.54  m/s 

Structure  Parameters 

c  2 

n 

from  optical  scintillometers 

2  x  10" 13  m"2/3 

C  2 

i 

2 

from  optical 

0.03°C2  m"2/3 

CT2 

from  resistance  wires 

0.03°C2  m"2/3 

V 

from  Lyman-a  hygrometers 

0.2  (g/m3)2  m"2/3 

s, 

from  resistance  wires  and 

0.075°C  (g/ra3)  m"2/3 

ah 

Lyman-a  hygrometers 

Si  2 

for  radio  frequencies* 

5.9  x  10"12  m"2/3 

2  2 

♦Obtained  from  Eq.  (79),  using  CT  ,  C  ,  and  C_  from  the  resistance-wire 
thermometer  and  Lyman-a  hygrometer  with  and  A^  obtained  from  the  radio 
refractive-index  equation. 
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Figure  4  shows  Che  phase-difference  PDF  for  antenna  pair  4  of  session  1, 
tape  24.  A  Gaussian  PDF  having  the  same  mean  and  variance  is  also  plotted. 
Clearly  the  phase  difference  is  Gaussian,  as  is  expected  theoretically. 

Figure  5  shows  the  structure  function  of  phase  versus  receiver  separation. 
Tape  20  corresponds  to  the  6:50  p.m.  data  and  tapes  22  to  27  correspond  to 
times  8:40  a.m.  to  3:50  p.m.  of  the  following  day.  The  six  receiver  separa¬ 
tions  correspond  to  the  six  antenna  pairs.  At  a  given  receiver  separation, 
the  value  of  the  structure  function  of  phase  is  just  the  phase-difference 
variance  for  the  corresponding  antenna  pair.  Also  plotted  is  a  theoretical 
prediction  for  the  6:50  p.m.  data  assuming  a  14-m  horizontal  outer  scale. 

Figure  6  shows  the  modulus  of  the  mutual  coherence  function;  the  3:50  p.m. 
data  of  tape  27  are  not  included.  Since  the  intensity  PDF  is  lognormal  and 
the  phase-difference  PDF  is  Gaussian,  it  follows  theoretically  that  the  MCF 
should  equal  exp(-D/2)  where  D  is  the  sum  of  the  structure  functions  of  log- 
amplitude  and  phase;  this  relationship  is  satisfied  to  within  0.5Z  in  the 
ram-wave  data.  Although  not  shown,  the  imaginary  part  of  the  MCF  is  very  much 
smaller  than  the  real  part,  as  is  expected  theoretically. 
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Figure  4. — Probability  distribution 
of  phase  difference  in  radians  for 
antenna  pair  4,  session  1,  tape  24, 
July  1983.  A  Gaussian  distribution 
is  also  shown.  This  distribution  is 
normalized  to  •«*it  area  under  the 
histogram. 
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Plgure  7. — The  normalized  variance  of  Intensity  versus  Its  Inertlal- 
range  formula.  The  straight  line  shows  theoretical  calculations. 
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